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Abstract 


The work presented in this thesis outlines the development of a new 
technique for visualization and for quantitative measurements of velocity, 
temperature and pressure of a gaseous flow. These are achieved by shining 
a single frequency laser beam into a flow which is seeded with an atomic 
species. The laser is tuned through the absorption frequencies of the 
seeded species and the absorption profile is detected by observing fluores- 
cence as the atoms relax back to the ground state. The flow velocity is 
determined by observing the Doppler shift in the absorption frequency. 
Spectroscopic absorption line broadening mechanisms furnish information 
regarding the static temperature and pressure of the moving gas. 

The primary motivation for the present study is the lack of an adequate 
diagnostic tool for high velocity flows. The technique is not limited to 
a particular flow density and is also nonintrusive since it does not intro- 
duce a local disturbance into the flow. Experiments have been conducted 
in the free stream and in the bow shock of a conical model mounted in a 
hypersonic wind tunnel. The results indicate that the experimental uncer- 
tainties in the measurement of average values for the velocity, temperature 
and pressure of the flow are 0.1, 5 and 10 percent respectively. 
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Chapter I 
INTRODUCTION 


The work presented in this thesis outlines the development of a new 
technique for visualization and for quantitative measurements of velocity, 
temperature and pressure of a gaseous flow. These are achieved by shining 
a single frequency laser beam into a flow which is seeded with an atomic 
species. The laser is tuned through the absorption frequencies of the 
seeded species and the absorption profile is detected by observing fluores- 
cence as the atoms relax back to the ground state. The flow velocity is 
determined by observing the Doppler shift in the absorption frequency. 
Spectroscopic absorption line broadening mechanisms furnish information 
regarding the static temperature and pressure of the moving gas. 

Traditional flow diagnostic techniques such as pitot probes, static 
pressure taps, total temperature probes, and hot wire anemometers intro- 
duce a disturbance at the measurement point in the fluid under investiga- 
tion.^ In contrast, the technique presented here does not introduce a 
probe into the flow and is therefore free of such local disturbances. 

Laser induced fluorescence from dyes has been found useful for both 

flow visualization and quantitative studies in incompressible fluid 
('2 3') 

flows. v * J For the gas phase, researchers have relied on electron beam 
fluorescence^ and electric discharge induced fluorescence.^ Of the 
two methods, electrom beam fluorescence has also given quantitative den- 
sity, temperature^ and velocity^ information in low density gases. 

f 81 

Index of refraction techniques' - . such as Schlieren, shadowgraph and in- 
terferometry have been applied to the gas phase for flow visualization and 
density measurements. Research is underway to extend these flow 
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visualization techniques to low density flows by using sodium vapor 

(9") 

seeding to enhance the index of refraction. ' 

Tracer particles in a flowing gas have extensively been used in con- 
junction with the laser Doppler Velocimeter (LDV) in subsonic and tran- 
sonic flows. In the hypersonic regime, however, these particles 

tend to lose their ability to follow the flow. The "particle lag"^^ 
problem accentuates with increasing flow velocities and lower densities. 

To overcome this difficulty, smaller particles must be used. As the size 
of the particles decreases, however, their effectiveness as a light scat- 
terer also decreases, thereby lowering the signal levels. The Resonant 
Doppler Velocimeter (RDV) technique presented in this thesis uses atomic 
size particles and therefore has virtually no particle lag problems. 

The RDV has been proposed by R. B. Miles in 1975 and first demon- 
strated in 1977. 

Figure 1.1 depicts the main features of the technique as applied to 
an idealized two level atom. When the laser frequency is tuned into the 
resonant absorption frequency in the frame of reference of the atoms, 
the atoms undergo an absorption process followed by fluorescence emission. 
The fluorescence signal is detected by a photomultiplier and recorded as 
a function of the laser frequency. The beam of the tunable laser is split 
in two directions. The lower beam in Figure 1.1 perpendicularly inter- 
sects an atomic beam of sodium. As the laser is tuned, maximum absorption 
and therefore maximum fluorescence is observed from the atomic beam when 
the laser frequency equals the atom's rest frame absorption frequency. 

This is due to the fact that the atomic and laser beams are perpendicular 
so the atoms have zero average velocity in the direction of the laser beam. 
The upper beam in Figure I. 1 intersects the seeded gas flow at an angle 
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so that the fluorescence peak occurs at a different frequency due to the 
Doppler shifted absorption line center. A comparison of the two frequen- 
cies furnishes a quantitative value of the Doppler shift from which the 
velocity component of the flow along the laser beam may be determined. 

An ensemble of two level atoms interact with an envelope of frequen' 
cies, an effect which is known as broadening of the spectral line. Var- 
ious factors including temperature and pressure have a bearing on the 
broadening. Hence, the resolution of the different broadening components 
permits the determination of the static temperature and pressure of the 
flow. 

Sodium is used as the seeded atomic species in this experiment for 
the reasons listed below. Sodium atoms possess two strongly absorbing 

( 14 ') 

line manifolds in the reddish-orange part of the spectrum. These 

lines are so strong that number densities of several atoms per cubic 
centimeter have been detected in the laboratory. Although sodium 

is not a two level atom, the structure of these lines is relatively un- 
complicated simplifying the interpretation of the spectrum. The asso- 
ciated absorption frequencies are close to the frequency corresponding 
to the peak of the lasing range of the Rhodamine 6G dye laser pumped by 
the 5145 A line of an argon ion laser. 

The RDV technique has been demonstrated in helium flows since helium 
does not react with sodium and the quenching of sodium by helium is neg- 
ligible. Helium can also be expanded to very high Mach numbers from 

room temperature since it lias a very low condensation temperature. This 
results in a low static temperature and pressure in the flow. Hence, 
the broadening is a small fraction of the Doppler shift, thereby enhanc- 
ing the velocity measurement accuracy. The sodium follows the helium 
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flow well since it equilibrates thermally after only a few collisions. 

Measurements of velocity, temperature, and pressure have been made 
in a free stream flow and in the vicinity of a model. The experimental 
components are described in Chapter II. Theoretical considerations pre- 
sented in Chapter III are used to optimize the different operating para- 
meters associated with the RDV and determine some of the limitations of 
the technique. The experimental results are presented in Chapter IV. 
Conclusions and suggestions for further study are discussed in Chapter V. 
The three sections of the appendix present a short discussion of atomic 
structure and broadening effects, statistical data analysis and the com- 
puter programs employed in the data reduction. 


Chapter II 

EXPERIMENTAL COMPONENTS 


A, Introduction 

The different elements included in the Resonant Doppler Velocimeter 
(RDV) experiment are presented in this chapter. The facility producing 
the hypersonic helium flow is described first (Section II. B). The pitot 
pressure and total temperature probe measurements are included in the 
same section. The measurements resulting from these classical flow diag- 
nostic tools are compared to those obtained from the RDV in Chapter IV. 
The design of the atomic beam device and the factors affecting its per- 
formance are presented in the following section (Section II. C). The 
sodium seeding apparatus is described next (Section II. D) . The optics 
employed in the experiment including the procedure for tuning the dye 
laser frequency to the sodium absorption frequency and monitoring its 
frequency stability are then discussed (Section II. E). The computer 
data acquisition scheme is presented in the final section (Section II. F), 
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B. Test Facility 
1 . Tunnel 

The hypersonic helium flow investigated in this study was generated 
at Princeton University's six inch (0.152 m) diameter hypersonic facility, 
He3, located within the Gas Dynamics Laboratory. A constant stagnation 
temperature was maintained by passing the test gas through a 61 m coil 
of room temperature steel pipe upstream of the stagnation chamber. The 
helium was then expanded through a converging diverging nozzle which was 
contoufed to provide a uniform free stream flow over the core region of 
the axisymmetric tunnel. All the data reported in this study were obtained 
at the set of operating conditions listed below. 

P 1 = 225±1 psi (1550±7 KPa) T* = 295°±1°K 

where P* is the stagnation pressure and T t is the stagnation temperature. 

All measurements were made at one axial location in the test section, 
at a position 0.900±0.003 m downstream of the nozzle throat or 0.216±0.001 m 
from the end of the nozzle. This position will henceforth be referred to 
as the measuring station. The various classical and spectroscopic data 
scans were taken across the wind tunnel diameter. The test section dia- 
meter at that point was 0.151±0.001 m. The free stream conditions, listed 
below, were determined from pitot pressure, stagnation pressure and tem- 
perature and the assumption of the isentropic expansion of a perfect gas. 

M n(l/m 3 ) T(°K) P(torr, Pa) u(m/sec) Re/£(l/m) 

13.8 7.5xl0 23 4.55 0.34, 46 173b 8.5xl0 6 

M is the Mach number, n is the number density, T is the static temperature, 

P is the static pressure, u is the velocity and Re/S- is the Reynolds num- 
ber per unit length. The run time of six minutes was limited by the 
steam ejector pumping system. 
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2. Pressure Measurements 

During a run the stagnation pressure was monitored using a 500 psi 
Heise gauge graduated in 1 psi increments. It was held at 225 psi (1550 
KPa) . The static wall pressure in the plenum chamber Q.61±Q.l m down- 
stream of the measuring station was monitored by a 800 torr Pennwalt 
(Wallace and Tiernan Division) Model FA160 gauge graduated in 5 torr in- 
crements. Both gauges are of the Bourdon tube type. 

Pitot pressure surveys were made at the measuring station with two 
pitot probes. For free stream measurements, a probe with a circular cross 
section was used. The inside diameter of the tube was 0.023 (0.58 mm) and 
the outside diameter was 1/32" (0.79 mm). The tip had the shape of an 
inverted truncated cone with a 10° half angle. Surveys with a model in- 
stalled in the test section on the tunnel axis were made with a pitot probe 
flattened to 1.0±0.2 mm external height and having an opening of 0.1±0.05 mm. 

Two Pace transducers, of the reluctance type, with a range of 5 and 
15 psi were used to monitor the pitot pressure. They were referenced to 
vacuum and calibrated using a mercury filled U tube manometer. The accuracy 
of the pitot pressure reading was 0.04 psi (270 Pa). 

3. Temperature Measurement 

The wall temperature was measured 0.09±0.0Q3 m upstream of the survey 
station using 20 gauge chromed alumel thermocouple wires fastened to the 
tunnel wall through a brass plug. The stagnation temperature was monitored 
with a model CASS-18G-12 Omega thermocouple. Both thermocouples used an 
ice reference. The voltages of the wall and stagnation chamber thermocouples 
were read by a Hewlett Packard model 425 AR and model 413 AR multimeters. 
These voltmeters have a recorder output through which the temperatures 
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were monitored with the computer throughout the run (Section II.F). 

A 0.005 inch (0.127 mm) diameter chromel/alumel thermocouple wire was 

used to measure the total temperature along a line extending from the wind 

tunnel wall into the core region. The spotwelded thermocouple junction 

was ground down to ensure that the wire, positioned perpendicular to the 

flow direction, would be of uniform diameter. This is a prerequisite for 

the unambiguous determination of the total temperature through the recovery 

ny) 

ratio. The recovery ratio, n, is defined by v ■ 



where T is the static temperature, T* is the stagnation temperature and 
T is the adiabatic wall temperature. In the continuum flow regime the 
value of the adiabatic wall temperature falls between T and T*. The tem- 
perature assumed by an insulated wall, past which flows a high speed gas 
stream, is controlled by physical phenomena in the boundary layer. Fluid 
layers far from the wall deliver viscous shear work to layers near the 
wall. The consequent temperature rise of the inner layer is necessarily 
accompanied by heat conduction away from the wall, thus tending to limit 

r is*) 

the temperature rise. '- J 

The temperature at the thermocouple junction is relatively unaffected 
by end losses (heat conduction to the prongs supporting the wire) for 
large enough values of the L/D ratio where L is the wire length and D its 
diameter (in the present work it is equal to 30). The supporting prongs 
were 0.005 inch in diameter which is smaller than is commonly used, reduc- 
ing the end losses even further. Hence it is reasonable to assume that 
the temperature measured by the thermocouple is the adiabatic wall tem- 


perature. 
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The recovery ratio is a function of the Knudsen number and in the 

fl9) 

hypersonic regime is independent of the Mach number. 

The nondimensional Knudsen number is defined by 

Kn = A II. 2 

where A is the mean free path and D is the dimension of interest, the 
thermocouple wire diameter. 

From kinetic theory, the viscosity, y, is dependent on the mean free 


path 


( 20 ) 


y = 3 y puA 


II. 3 


where p is the density, @ y = 0.5^^ and the average atomic speed is^*^ 


“ = / 8kT 
/ Trm 


II. 4 


where k is the Boltzman's constant, T is the static gas temperature and m is 
the mass of one atom. By substituting equations II. 3 and II. 4 into II. 2 
one obtains 


» _ h _ nr m_ 

n / 2 Re 

where y is the ratio of specific heats and M and Re are the Mach and 
nolds numbers respectively. These results are used in Chapter IV to 
pret the total temperature surveys. 


II. 5 

Rey- 

inter- 


- 11 - 


C. Atomic Beam Device 

An atomic beam device is used to create a directed beam of collision 
free particles. In the present experimental work it is used as a frequency 
reference. According to kinetic theory, the motion of an atom in a gas 
may be regarded as purely random and is accompanied by continual collisions 
of the atoms with each other and any surfaces present. Between collisions 
the atoms describe straight paths with uniform velocities which are distri- 
buted according to Maxwell's law. Corresponding to the average velocity 
there exists a mean free path A which is inversely proportional to gas 
pressure, 

If a circular aperture with diameter less than A connects a gas filled 
vessel, with a highly evacuated chamber, the atoms colliding within one 
mean free path from the aperture and having velocity vectors directed to- 
wards the aperture will have a high probability of escaping into the evac- 
uated region. The atoms then continue in straight paths, filling a solid 

( 22 ) 

angle which is determined by the collimating properties of the aperture. 

The atoms which pass through a second aperture downstream of the first are 
constrained to move along nearly parallel paths defined by the two aper- 
tures and the distance between them. This forms an almost collisionless 
beam of atoms. Collisions in the highly evacuated regions occur only in 
those rare occasions when faster atoms overtake slower ones travelling 
along the same path. 

The velocities of the atoms in the beam depend on the temperature of 
the source. Their distribution is not quite Maxwellian since the method 
of beam formation favours faster atoms. This is due to the fact that col- 
lisions taking place in and near the source aperture must be taken into 
account even when the mean free path in the source is considerably greater 



- 13 - 


than the aperture size. The probability of loss of atoms from the beam 

by this mechanism increases rapidly as the atomic velocity becomes less 

than the most probable velocity of the atoms in the source. This accounts 

r 231 

for most of the low velocity deficiency. J The agreement between the 
Maxwellian and the experimental velocity distribution is quite good for 
velocities above the most probable velocity. 

Figure II. 1 is a schematic of the apparatus employed in the ^ jsent 
study. Sodium metal contained in a quartz crucible is placed in a tungsten 
basket. A current applied to the basket heats the sodium to about 180°C. 

The temperature is monitored by a chrome 1/alumel thermocouple in physical 
contact with the crucible. The vaporized sodium atoms leave the surface 
in all directions, some of them reaching the bottom aperture. The top aper- 
ture produces a collimated beam in the test section. This atomic beam is 
intersected perpendicularly by the dye laser beam entering and exiting the 
device through Brewster angle windows. The detection optics are positioned 
perpendicularly to both the laser and atomic beam. A lens images the fluo- 
rescence signal onto a thermo-electrically cooled (Product for Research 
Model TE-104) RCA Model C31Q34 photomultiplier tube as depicted in Figure 
11 , 2 , 

Chambers B and C are evacuated by a Vactronic Econovac model 20 dif- 
fusion pump through ports welded to the walls of the device. Each vacuum 
line is immersed in a liquid nitrogen cold trap before joining at the pump. 
This prevents the circulation of sodium atoms between the chambers and en- 
hances the vacuum by nearly an order of magnitude. The pressure recorded 

-7 -5 

by a Vectronix model LDG-266 cold cathode gauge was 5x10 torr (7x10 Pa) 
at the pump's inlet port. The pressure in Chamber B, monitored with a 
Norton Vacuum Equipment model NRC 831 ion gauge using a NRC 507 tube 
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was 10 _o to 10"^ torr (1C) to 10” Pa). These pressure readings were 

taken with the sodium source at the operating temperature. At 180°C the 

-4 -2 (44) 

vapor pressure of sodium is about 10 torr (10 Pa). This corres- 

ponds to a mean free path A of*' 20 '* 


A = -J- - ■ = — i — = 3,5 cm II. 6 

/2Vd n /2o n 
c 

where a ~ 10~ 13 cm 2 *- 24 '* is the elastic collision cross section for a sodium 
c 

sodium collision and n is the number density. (The helium-sodium collision 
cross section is one order of magnitude smaller.) As the distance from the 
crucible to the aperture is 22 mm and the diameter of the bottom aperture 
is 0.030" (0.76 mm) the majority of the sodium atoms leave chamber A with- 
out undergoing a collision. This insures some directionality to the beam 

_ A 

in B and enhances its intensity. In B, at a pressure of 5xl0~ torr (7xl0~ 
Pa), the mean free path is of the order of 400 cm, which is appreciably 
longer than the upper aperture diameter of 0.020" (0.5 mm) and the distance 
to the observation region of 3.2 cm. These design conditions lower the 
probability of a collision taking place in the probe volume. The atomic 
beam possesses some divergence due to the finite dimension of the apertures 
and their separation. The maximum divergence angle a for a non-colliding 
atom is 


a = tan" 1 [(d 1 +d 2 )/£] II. 7 

where dj, d 2 and Z are the lower and upper aperture diameters and their 
separation respectively. For a separation of 9 = 8.6 cm the angle is 
a = 50' 56". At worst this would correspond to a broadening half width 
half maximum (HWIIM) , Av, of 


Av = v_, — sin a 
21 c 


/SkT 
/ Tim 


v^sin a 
c 


II.S 
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where is the stationary absorption frequency, u and m are the sodium 
atom thermal velocity and mass respectively, c is the speed of light and 
T the temperature of the source. For a source temperature of 1S0°C and 
the sodium transition wavelength of 589 nm one obtains a HWIIM of 16 MHz. 

A HWIIM of 10 MHz is obtained from the spectrum of the sodium line cor- 
responding to the 589 nm wavelength. The natural and laser HWIIM are 5 MHz 
each. Their sum is 10 MHz which equals the experimental width obtained 
from the atomic beam device spectra. To sum the natural and laser line- 
widths is a simplification which is not rigorously correct due to the fact 
that the laser linewidth is not Lorentzian (Section IV. 4, Appendix A). 

Since the sum of the natural and laser linewidth equals the linewidth ob- 
tained from the spectra of the atomic beam device it is clear that the 
broadening contribution due an uncollimated atomic beam with an upper bound 
calculated by using equation II. 8, is negligible. If the laser beam does 
not intersect the atomic beam perpendicularly, a frequency shift and addi- 
tional broadening takes place. This is due to the fact that the atoms 
have a velocity distribution along the atomic beam. The broadening effect 
is negligible for small deviations from 90° but the resulting Doppler shift 
of the central absorption frequency has to be taken into account for velo- 
city measurements (Section IV. B. 3). 



v 


Figure II. 3 


D. Tunnel Seeding Apparatus 


To facilitate seeding the flow through the stagnation chamber at high 
pressure the apparatus shown in Figures II. 3 and II. 4 was devised. A small 
amount of purge helium is used to produce a hot helium sodium mixture which 
is then seeded into the main helium supply. The oven is maintained exter- 
nal to the small stagnation chamber with only a needle penetrating into it. 

The helium purge supply is preheated while passing through a 1/4" stain- 
less steel coiled tube of 0.35" wall thickness (part 23) and then flows into 
the oven. In the oven the hot helium gas is forced to sweep very close to 
the sodium contained inside a quartz crucible (part 37) . It therefore mixes 
with the sodium vapor and the mixture then flows towards the nozzle through 
the inner needle (part 6). This is a stainless steel 1/4" tube of 0.02" wall 
thickness which again is resistively heated. 

The oven constitutes a heat sink as it contains a relatively large 
amount of metal. Due to its surface area it loses heat to the surroundings. 
To compensate for the above two losses an additional, independently powered, 
heating element (part 18) is installed inside the oven around the crucible 
containing the sodium. 

The resistance of all heating elements are calculated so as to utilize 
the maximum power capability of the power supplies. A Mallory type VA-6000 
DC power supply capable of 9600 watts at 48v is used for the coil. A 
Harrison Laboratory model 810-A capable of 450 watts at 60 volts DC and 
a Varian model 927-007 with a power output of 3 KVA at 20 volts AC are used 
to heat the oven and needle respectively. Figure II. 5 is a schematic of 
the power and helium supplies. The voltage to the coil is applied between 
parts 26 and 45. The voltage to the inner needle is applied between parts 
45 and ground. The outside jacket (part 22) is insulated electrically from 
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the resistively heated elements by an electrical insulator (part 27) . 

Thermal insulation is provided by the evacuated chamber between the inner 
and outer needles (parts 6, 5) and between the oven and coil assembly and 
the cylinder (part 22) . The pressure inside the oven is monitored by a 
500 psi Heise gauge of the Bourdon type graduated in one psi increments. 

It is maintained at 25 psi (0.17 MPa) above the stagnation pressure using 
a differential pressure regulator with a range of zero to 90 psi. 

The temperatures of the preheated helium, of the bottom of the sodium 
crucible and of the needle are all monitored by chormel/alumel thermocouples. 
Two junctions are attached to the inner needle, one at about 25 mm from the 
tip on part 1 and the other at 5 cm from the tip on part 6. The oven thermo- 
couple is operated in conjunction with an Omega model 50K on/off temperature 
controller with a 0-600°C range at 5°C increments. The thermocouple junc- 
tion fastened to the 1/4" (0.0635 m) stainless steel needle (part 6) is 
connected to a IVheelco model 401C on/off temperature controller with a 0- 
1000°C range at 10°C increments. The purge flow is obtained from industrial 
grade helium delivered in bottles. 

The tunnel starting procedure consists of several steps. First, the 
helium purge flow through the sodium loaded oven is established while the 
tunnel is evacuated by a 40 CFM Stokes pump. The power supplies to the 
seeding device are turned on until fluorescence is observed in the test 
section. At this stage the tunnel is started by gradually raising the stag- 
nation pressure. The pressure in the oven increases simultaneously, con- 
trolled by the differential pressure regulator. Due to the large amount 
of the main helium flow through the stagnation chamber when the tunnel is 
operating, additional power has to be supplied to the seeding device to com- 
pensate for the convective heat losses. 
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During a run, the needle tip is maintained at approximately 700°C 
while the oven temperature is set to 250°C. It was found that the preheat- 
ing coil does not have to be energized at all since enough heat is avail- 
able from the other heated components. To determine the helium purge mass 
flow rate through the oven, the pressure drop in one of the supply cylinders 
is monitored as a function of time. If internal leaks in the pressure regu- 
lator are neglected then the mass flow rate, M, is 


• AM IVV AP 

At RT At 


II. 9 


where M and IV are the mass and atomic weight of sodium respectively, V and 

T are the cylinder volume and temperature, AP is the pressure drop in time 

At in the cylinder and R is the universal gas constant. For the seeding 

• -5 

device injecting into vacuum, M = 4.3x10 Kgm/sec. When the tunnel is 
operating, M = 3.1x10 Kgm/sec. The difference is due to the fact that 
the needle tip is located inside the converging part of the nozzle and 
senses a pressure which is lower than the stagnation pressure when the 
tunnel is operating. This raises the flow rate through the oven. 

To enhance the signal to noise ratio for the collection optics focused 
into the flow, the sodium density should be as high as possible. The limit, 
imposed by radiation trapping considerations is 8.2x10^ atoms/cc if one 
demands K^L = 0.1 as determined in Chapter III. Clearly 


M. 


Na 


n XT iti x . 
Na Na 


n.. m. , 

M^ e He He 


II .10 


where M is the total mass, n the number density and m is the mass of one 
atom. The helium number density is 


n^ e = 2.69x10 


19 P/P o 
T/T 


atoms , _ ._17 atoms 
= 7.5x10 
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where T q = 273.15°K and P q = 760 torr (1.013xl0 5 Pa). T and P refer to 
the temperature and pressure in the test section. The mass flow in the 
tunnel is 0.0462 Kgm/sec based on the isentropic flow assumption and the 
physical throat area of the nozzle. Hence, the sodium mass flow rate equals 
2.9x10 ** Kgm/sec. The above calculation assumes a uniform helium density 
across the test section and uniform mixing of the sodium atoms in the helium 
flow. Thus one obtains a seeding ratio of 6.3xl0 -4 parts per million on a 
mass basis. 
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E. Optics 

A Spectra Physics model 580 single frequency tunable dye laser is 
ppmped by a Coherent Radiation model 53A Argon ion laser. The pump laser 
is rated at 2.5 watts output at 5145 A. At 1 watt input, the dye laser 
output ranges from 30 to 70 milliwatts, dependent on the quality of the 
alignment. A schematic of the setup is depicted in Figure I I. 6. 

The output beam is split in several directions. The Spex 1700 Czerny- 
TUrner monochromator is used for coarse wavelength adjustments. The obtain- 
able resolution according to the manufacturer is 0.1 A or 10 GHz at the 
sodium line corresponding to the 589 nm wavelength. A sealed pyrex ampule 
containing sodium metal under vacuum is used for finer wavelength adjust- 
ments. The cell is heated to raise the atomic vapor pressure. The resolu- 
tion is determined by the linewidth of the sodium absorption profile in the 
cell. In this case the temperature broadening predominates and at 150°C 
one obtains a HIVHM of 782 MHz. The sum of twice the HIVHM plus the ground 
state hyperfine splitting provides a resolution of about 3.5 GHz. As the 
dye laser is capable of a 4 GHz linear frequency scan, it is relatively 
easy to find the sodium lines in the atomic beam device once fluorescence 
in the sodium cell has been observed. The longitudinal modes and frequency 
stability of the dye laser are monitored by a Jodon model SA-1500 confocal 
spectrum analyzer with a free spectral range of 1500 MHz and a linewidth 
of 7 MHz, A photodetector measures the laser intensity to record laser 
intensity fluctuations during a scan so that these may be removed during 
the data processing stage. 

One of the beams is directed into the sodium seeded tunnel. The fluo- 
rescence signal from the seeded atoms is picked up by an EMI type 95745 
photomultiplier tube. Fgure II. 7 is a schematic of the tunnel collection 
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optics. Lj is a 113 mm focal length, 53.3 mm diameter lens. The distances 
from the probe volume at the measuring station to and from to the 
0.0135" (0.34 mm) diameter aperture Aj are 22.2 and 23.2 cm respectively. 
The detection acceptance solid angle is defined by 

0 2tt 

0, = f sin 6d0 f d<f> = 2tt(1-cos 0) 11.12 

o o 


where 

diameter of L. 

tan ® ~ 2x(distance from probe volume to L^) " tan ^ ^0 ^ ) 11.13 

_3 

The collection efficiency, £2/ 4tt, equals 3.6x10 . The probe voluri at the 

-4 3 

measuring station is 1.8x10 cm for a laser beam perpendicular to the 
collection optics axis. This volume is defined by the laser beam radius at 
the e 1 points and the collecting aperture diameter, A^. The two lenses 
between the aperture and the photomultiplier tube were matched so as 
not to degrade the collection efficiency of L^. These lenses were employed 
to enable the use of an interference filter as depicted in Figure II. 7. 

The collection optics are mounted on an xyz translator. The laser 
beam position relative to that of the probe volume inside the tunnel defined 
by the aperture A^ is kept constant by passing it through a mirror sequence 
as depicted in Figure II. 8. Since the beam enters each translation stage 
parallel to its direction of movement, the following stage does not perceive 
any change in beam orientation. All the data presented in this experimental 
work was obtained using the xyz translator. Another possible configuration 
would be to expand the beam into a sheet of light and scan the position of 
the detection optics only. The disadvantage of this method is that while 
the signal is collected from the probe volume only, scattered light is con- 
tributed by the whole width of the beam. 
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F, Data Acquisition Scheme 

A Hewlett-Packard 1000 16 bit minicomputer is incorporated into the 
experiment for data acquisition as shown in Figure II. 9. Two A/D conver- 
ters are used since all the data are in analogue form. A four channel 12 
bit Preston model GMAD/1 with a variable sampling rate of 350 to 500,000 
Hz for all channels is used to measure four inputs: the dye laser frequency 
scanning voltage ramp; the signals out of the two photomultipliers used to 
monitor the fluorescence intensities from the atomic beam device and the 
tunnel; and the output of the photodetector used to probe the laser inten- 
sity. The data from the set of four inputs are read simultaneously within 
a window of 10 nsec. A thousand points per channel at a rate of 530/4 Hz 
are obtained per laser frequency scan. Due to nonlinearities in the elec- 
tronics producing the high voltage ramp driving the dye laser frequency 
scan, one obtains only about 600 useful points per channel. Hence the 
scan time is 7.3 sec. 

Three additional channels are connected to the 14 bit Preston model 
GMAD/4 A/D converter which has a fixed total sampling rate of 1750 Hz and 
can read up to 64 channels sequentially. The inputs are the tunnel stag- 
nation temperature, the test section wall temperature and the position 
along one axis of the xyz collection optics translator. The axis chosen 
determines the position of the probe volume along the tunnel diameter. 

The input range of both A/D converters is ±10 v. Hence, to take ad- 
vantage of the full resolution of these instruments, the weaker signals are 
passed through D.C. amplifiers. A Princeton Applied Research model 115 
preamplifier is used for the laser intensity and two Dymec. amplifiers 
model 2460A are used for the two photomultiplier signals. 
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A Tektronics R564B storage oscilloscope with a dual trace vertical am- 
plifier is used to monitor the two fluorescence signals from the photo- 
multipliers during the laser scan. The horizontal amplifier is driven by 
the laser frequency scanning voltage ramp. 


[ Chapter III 

(■j, 

RESONANT DOPPLER VELOCIMETER OPERATING CONDITIONS 
A. Introduction 

A calculation for the determination of the optimal operating conditions 
for the RDV is presented in this chapter. The calculation is representative 
only since it applies to the He3 facility in the Gas Dynamics Laboratory 
of Princeton University as described in Section II. B. 

The procedure for deriving the velocity from the Doppler shift is dis- 
cussed in the beginning of the following section. The effect of the number 
density of the seeded sodium on laser light absorption and radiation trap- 
ping is presented next. A discussion of the limit which has to be imposed 
on the laser power level to reduce saturated absorption and power broaden- 
ing follows. Finally, the expected signal level based on the optimal sodium 
number density and laser power level is calculated. The phenomena of the 
sodium atoms lagging behind the helium flow, optical pumping and radiation 
quenching are also touched upon briefly. 
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B. Velocity Calculations 


The velocity of the sodium atoms can be determined from the Doppler 
shifted absorption frequency. From equation A. 29 

IH. 1 


-► 

u % 
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A 

where k is a unit vector in the direction of the propagation vector k, u 
is the atomic velocity vector and c is the speed of light. The equation 
may be rewritten in terms of 
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Expressed in terms of a Doppler shift, Ato, the calculated velocity 


component, u, parallel to the propagation vector, k, is 


u = Ato 


“ki 


Ato , 

2tt A ki 


III. 4 


where is the absorption wavelength and Ato equals to-to^. 

The Doppler shift ^s measured by simultaneously scanning the dye laser 
frequency through the stationary absorption line in the reference cell and 
through the shifted absorption line in the wind tunnel while monitoring the 
fluorescence intensity from both. Since the ground state splitting is of 
the order of the velocity shift, one may use it to obtain a relative laser 
frequency calibration as detailed in Chapter IV. For accurate measurements 
of the Doppler frequency shifts, the reference frequency has to be deter- 
mined as accurately as possible. The collimated atomic beam in the refer- 
ence cell is collision free, hence the absorption linewidth is limited only 
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by the natural and instrumental broadening as discussed in Section II.C 
and Appendix h. The number of collisions a sodium atom undergoes in the 
tunnel during the time interval between injection and detection is an in- 
dication of how closely the sodium atom follows the helium flow. The 
smallest collision frequency and largest convective velocity should be used 
to obtain a lower bound on the number of collisions. These conditions are 
realized in the test section of the wind tunnel. The collision frequency 
per sodium atom, z, is 


z = n,. a 


/ 8kT 
He u c / fry 


A. 58 


where n u is the helium number density, a is the sodium helium collision 
cross section, y is the reduced mass, k is the Boltzman constant and T is 
the temperature. All the experimental investigations performed to determine 
the collision cross section for the sodium lines were carried out at tem- 
peratures of 400°K and higher. A cross section for the temperature range 
of only a few degrees Kelvin is measured in this experimental investigation 
(Section IV.D.l). For T = 4.55°K and n^ e = 7.5x10“'’ 1/m^ determined from 
a pilot survey (Section IV. B. 2) then z = 142 MHz. The cross section, 
c,= 126x10 ^ cm“, obtained by McCartan and Farr^^ at T = 415°K was used 
for a first iteration. 

The velocity in the test section based on isentropic flow calculations 

is 1736 m/sec and the distance from the injection point to the probe volume 

-4 

is 0.95 m. A sodium atom therefore travels 5.4x10 sec before being 

detected. At the above calculated collision frequency, it undergoes at 
4 

least 7.7x10 collisions. One may thereby assume that the sodium atoms 
seeded in the tunnel follow the helium flow closely since they undergo many 
collisions between the point of their injection into the flow and the probe 


volume. 
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The two strong sodium lines with transition frequencies in the fre- 
quency range of the Rhodamine 6G dye laser are the at 5896 A 0 
(3^S 1/2 "► s2p 1 / 2 ) and the D 2 at 5890 ^ C3 2 Si/ 2 ^ 3 2 P 3 ^ 2 ) transitions. A 
manifold of four and six hyperfine lines is associated with the Dj and D 2 

2 2 

transitions respectively (Appendix A) . The D 2 3 S l/2* F=2 3 P 3/2 * F=3 

hyperfine line is used for a frequency reference since it is almost three 
times more intense than the next strongest line. 
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C. Laser Power and Sodium Number Density Limits for RDV in Helium 

The temperature and pressure may be determined by fitting a Voigt 
profile to the experimental points (Appendix A) . The broadened lineshape 
may be distorted compared to the theoretical one by strong laser absorp- 
tion, radiation trapping and power broadening in addition to optical pump- 
ing (Appendix A). These effects may be reduced by controlling the amount 
of seeded sodium and laser irradiance. 

From the pitot survey and the stagnation conditions of 295+1 C K and 
2?5±1 psi (1550±7 KPa),one obtains T = 4.55+0.09% P = 6.7±0.2xl0 -3 psi 
(45±2 Pa) and M = 13.8±0.08 for the temperature pressure and Mach number 
existing in the test section. From equation A. 31 the Doppler HWHM, Av^, 
is expected to be 81.8 MHz. From equation A. 72 the collision HWHM, Av c> 

is expected to be 12.8 MHz where the cross section at a temperature of 
.(25) 


415 K 


has been used for a first iteration. The natural, Av , and 

n 


laser, Av^, are 5 MHz each. The total Lorentzian HWHM, Av^, is 


Av, = Av + Av 
L c n 


III. 5 


Assuming that the laser HWHM can be simply added to the Lorentzian component 
then the total Lorentzian HWHM 


Av, = Av + Av + Av„ 

T n r\ C 


22.8 MHz 


c n l 

The ratio of Lorentzian to Gaussian linewidths, a, is 
Av, 
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a = 


Av G //£n2 


0.252 
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where Av^ is the Gaussian HWHM and Av^ equals Av^. From tables based on 


numerical calculations 


(26) 


one obtains the field free Voigt HWHM of 


Av y (I=0) =94.7 MHz 


III. 7 
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If a 1% error in the HWHM of the Voigt profile due to power broadening 
is acceptable, then 

AUyfl t 0) = 1.01 Av y (I = 0) = 95.6 MHz III. 8 

and the linewidth ratio, a', is 0.248, The power broadened Lorenzian 
HWHM is 


Av 

Av'(I + 0) = a' — 

L /Sn2 


= 24.4 MHz 


III. 9 


The saturation intensity, is 


Srr hcAv 


h^21> = 


L = 104.5 W/m 2 


A. 98 


1 + 


21 


where g is the degeneracy, h is Planck's constant, c is the speed of light, 
and is the transition wavelength. The subscripts 2 and 1 refer to the 

upper and lower states respectively. The laser intensity, I, is related 
to the saturation intensity by equation A. 100 and may be calculated for 
the above conditions 


A^t ( I 5^ 0) ™ a' fl ^ 01 2 9 

1 = I s (w 21 ){[— ^ Tfr r- l} = I s ^ aJ 2i :>{[ ~ a 3 -1} = 15 w/m 2 III 
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The intensity distribution for a cylindrical beam with a Gaussian profile 

. (27) 
is 


I' = I e 
o 


-2 (— ) 2 

4 V dP 

" dA 


where w is the beam waist at the e~ point and r is the radial distance 
from the laser beam axis, Therefore the power is 

2 


Hi w 


P = / I ' 27Trdr = 


-5 2 

The power equals 2.4x10 w assuming I Q equals 15 W/m at a beam waist of 1 mm. 
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This then is the maximum total laser power that can be used to illuminate 

a 2 mm diameter cross section area in the flow without causing distortion 

of the data due to saturation and power broadening. The laser beam may 

alternatively be expanded into a sheet of light since enough power is 

available at the output of the dye laser to illuminate a large portion of 

the flow. The calculation was based on the values &2 = ^ anc * = 3 which 
2 2 

hold for the 3 F=1 3 P 3/2’ F=2 trans i t: *- on * This transition pro- 

vides the lowest saturation intensity. Hence, if the laser power is kept 

below P = 2.4x10*5 IV (for a 1 mm waist), all other lines connecting levels 
2 2 

3 Sj /2 anc * 3 w iH possess an even smaller saturation broadening com- 

ponent . 

The sodium D 2 transitions have different strengths (Figure A. 2). The 

intensity of the laser beam at the probe volume is therefore dependent on 

the hyperfine line being pumped due to absorption along the beam's path as 

depicted in Figure 111,1. A stronger line would absorb more of the laser 

beam power through its passage in the sodium seeded region towards the 

probe volume. Also, significant absorption of the fluorescence signal may 

take place between the probe volume and the detection optics. The last 

effect known as radiation trapping is also dependent on the line strength. 

If major differences in absorption are allowed for different transitions 

the interpretation of the spectra would become difficult. In order to 

avoid this problem the density of the sodium must be limited. 

Due to broadening effects in the tunnel , the structure of the upper 
2 

state 3 P is not resolved. The relative strengths of the two transi- 
tions 3 2 S jy 2 , F=l,2 3 2 P^2 are 12/48 and 20/48 respectively. From 

equation A. 82 the ratio of the laser irradiances, after passing a distance 
L through the absorbing medium for an optical depth, K^L, equal to 0.1 is 
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given by 
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One may calculate the maximum number density of sodium for an optical 
depth of 0.1 by using equations A. 106 and A. 111. 
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where is the absorption frequency, c is the speed of light, is 
the transition frequency, g is the degeneracy, A is the Einstein coeffi- 
cient for spontaneous emission and n° is the zero field number density, a 
was defined in equation A. 108. The detuning, x, is 


v-v 
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21 
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The term in square brackets is the Voigt integral (Appendix A) . On 
line center x = 0. Then 
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Use has also been made of equation A. 15 and A. 16 to obtain 
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where r = 2.82xl0 - ^ m is the classical electron radius, f, „ and f„, are 
o 12 21 

the absorption and emission oscillator strengths and T is the natural HWHM. 


From equation III. 15 for L on the order of 0.1 m and K^L = 0.1 
o o 

n i n ? 13 3 

— = 8.19X10 1 atom/in 

Si 
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At equilibrium, for X = 589x10 m and a temperature of 4.55 K J 
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The thermal upper state number density may be neglected. Hence 
0 

n l 13 3 

— = 8.19x10 atoms/m 

*1 
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The laser power available at the probe volume varies less than 1.6 ? o for 

2 

the two ground states of 3 5^/2 ^ soc ^ um density is kept below the cal- 
culated value. The ratio of the detected fluorescence emission intensity 
of the stronger transition manifold to the weaker transition manifold, if 
trapping is included for the same optical depth of 0.1 is from equation 
A. 82 
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where the first bracket relates to the absorption of the laser beam and 
the second bracket to radiation trapping. The ratio obtained from an in- 
finitely short absorption length equals the ratio of the relative strengths. 


20/12. The difference between the two ratios is 2.5%. For k L = 1 and 


w 
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0.01 the ratios are 1.304 and 1.662 which differ from the ideal value by 
22% and 0.25% respectively. In the present experimental investigation 
the optical depth is of the order of 0.1. 

The power absorbed by the sodium atoms is given by equation A. 87. 

It also equals the energy of a single photon multiplied by the net photon 
absorption rate. The number of photons emitted equals the number absorbed, 
times a radiation quenching factor defined by 

O - Intensity of fluorescence with foreign gas 

' ” Intensity of fluorescence without foreign gas 

Quenching is a result of nonradiative decay of the upper state due to in- 
elastic collisions with foreign atoms or molecules. Quenching occurs when 
the atomic excitation energy is converted to translational, rotational or 
vibrational energy. The rate per unit volume of spontaneously emitted 
photons, n, is 


• I 

n = fiw 



19 

4.5x10 


1/sec-m 
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The values used for this calculation are the sodium D 2 transition frequency, 
2 

I = 15 IV/m , K^L = 0.1, and L = 0.1 m. Sodium quenching by helium is neg- 
ligible^ J and therefore Q * 1. 

The detection optics are discussed in Chapter II. The rate of col- 
lected photons, N, if the probe volume and the collector solid angle are 
taken into account, is 

N = n • V • ~ III. 23 

From the above calculations, the sodium density has to be kept below 
13 3 

8.2x10 atoms/m to eliminate the distortion of the absorption lineshape. 
The laser power has to be kept below 2.4x10 ^ W so as not to add an appre- 
ciable saturation component to the Lorentzian linewidth. Under these 
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7 

conditions one obtains 2.9x10 photons/sec for a probe volume of 1.8x10 

3 -3 

m and a collection efficiency of 3.6x10 (Section II. D). 


Chapter IV 
RESULTS 


A. Introduction 

The experimental results are presented in this chapter. Pitot pres- 
sure and total temperature surveys were taken to determine the effect of 
sodium seeding on the helium flow in the wind tunnel. The characteristics 
of the atomic beam and the laser spectral purity are then analysed since 
they affect the quality of the data obtained from the wind tunnel. A 
discussion of the data analysis procedure is presented next and is followed 
by the results obtained from three different experimental configurations. 

The atomic beam device is used for velocity measurements in the free 
stream flow of the wind tunnel in the first experiment. A laser double pass 
configuration in the free stream flow of the wind tunnel is employed ‘next 
and dispenses with the use of the atomic beam device. The first config- 
uration is applied in the third experiment to perform measurements in 
the vicinity of a shock wave surrounding a conical model placed in the 


wind tunnel. 
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B. System Calibration 
1. Total Temperature Surveys 

Two total temperature surveys were taken at the measuring station 
(Section II.B.l), without and with sodium seeding through the needle, as 
depicted in Figure IV. 1 and IV. 2 respectively. The ordinate represents 
the adiabatic wall temperature and the abscissa is the distance along 
the wind tunnel diameter from the wall into the hypersonic core. 

The adiabatic wall temperature measured in the free stream without 
sodium injection was 278±2°K. Hence the recovery ratio, ri (equation II. 1) 
is 0.95±0.008. This value is based on a measured stagnation chamber tem- 
perature of 295±1°K and a static temperature of 4.55±0.09°K calculated 
from isentropic relations based on the stagnation temperature and the Mach 
number . 

The Knudsen number, Kn (equation II. 2), obtained in this experimental 
work is 0.02 where the Reynolds number, Re (equation II. 5), is based on the 
wire diameter. The recovery ratio obtained by P. Materna^^ for the same 
Knudsen number was 0.96. 

The injection needle tip was maintained at 700°C while the survey of 
Figure IV. 2 was taken. The injection system is responsible for some heat- 
ing of the center of the flow. The heating manifests itself through a 
peak in the core region as depicted in Figure IV. 2. The total temperature 
at the peak is 305±3°K (taking into account the recovery ratio). It is 
also clear that some asymmetry is present in the flow, as discussed in 
the following section. 
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2. Pitot Surveys 

The measuring station, as discussed previously (Section II.B.l), is 
maintained at 0.216 m downstream of the nozzle exit because of experimen- 
tal constraints. This position facilitates scanning the optics along the 
tunnel diameter as far as possible towards the tunnel wall. Most of the 
previous experiments performed by other workers in the He3 facility were 
conducted closer to the nozzle exit than the measuring station. Hence, 
two sets of pitot surveys were taken to compare the flow conditions at the 
measuring station (Figure IV. 3 to Figure IV. 5) to those existing at 
0.049±0.001 m downstream of the nozzle exit (Figure IV. 6 to Figure IV. 8). 
Figure IV. 3 and IV. 6 are the surveys obtained with the injection needle 
removed, corresponding to the normal mode of operation of the wind tunnel 
without sodium seeding. The injection needle was positioned 0.1 inches up- 
stream of the nozzle in Figure IV. 4 and IV. 7 with no injection. Sodium 
injection was taking place in Figure IV. 5 and IV. 8. The abscissa is the 
distance from the tunnel wall and the ordinate is the pitot pressure. 

There is no apparent difference between Figures IV. 7 and IV. 8 and between 
Figures IV. 4 and IV. 5 apart from the vicinity of the compression wave. 

The concave part of the expanding side of the nozzle produces compression 
waves as depicted in Figure IV. 9. The pitot probe which is scanned across 
the tunnel diameter, 0.216 m downstream of the nozzle exit, senses these 
compression waves and registers them as a pressure peak on both sides of 
the inviscid core. The compression wave really surrounds the inviscid core 
since the tunnel is axisymmetric. It appears that the presence of the needle 
in the stagnation chamber increased the diameter of the inviscid core in 
the test section and also reduced the compression wave strength near the 
core boundary. This effect was particularly noticeable at the measuring 







inches Hg)l(KPa) 



r igurc 



X(cm) 








-55- 


station. The pitot surveys also indicated the presence of flow asymmetry 

as previously observed with the temperature survey. The compression wave 

surrounding the inviscid core is not of uniform magnitude. 

To determine the Mach number from the Pitot pressure, p!j> one has to 
( 311 

calculate the ratio v J 


P 2 P 1 P 2 P ? 

4* <4 (/) c4) 

pL p L *1 r 2 

*1 1 


IV. 1 


P with the superscript t is a stagnation pressure and without it a static 
pressure. The subscripts 1 and 2 refer to the pressure upstream and down- 
stream of the shock standing in front of the pitot probe. The expansion 
from the stagnation chamber towards this shock is not isentropic when 
sodium seeding takes place as is evident from Figure IV. 2 since the flow 
is not adiabatic. Then dlearly P* is not the stagnation chamber pressure 
but an unknown quantity. Therefore when the flow is seeded the Mach number 
may not be found from equation IV. 1 by substituting the stagnation chamber 
pressure for P*. The pitot pressures measured in the hypersonic core for 
all the six different surveys are the same to within the experimental error 
and are equal to 3.9±0.08 inches Hg (13200±270 Pa). The four surveys per- 
formed with no seeding taking place may be used for Mach number determin- 
ation. The Mach number obtained in the inviscid core, based on the pitot 
and stagnation pressure is M = 13.8+0.1 assuming an isentropic expansion 
of a perfect gas. The area ratio of the nozzle is 2.7x10 ^ which corres- 
ponds to a Mach number of 17.95. The nozzle is designed for a Mach number 
of 16. The difference between the value of 17.95 and 16 is due to the fact 
that the boundary layer thickness is taken into account which limits the 
expansion. The measured Mach number of 13.8 is lower than the design 
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value of M = 16. To determine whether this was due to contaminated helium, 
the same gas was used in the He2 facility and in the He3 facility with its 
original nozzle and with the He4 nozzle. All these nozzles were designed 
to operate at Mach 16. The resulting pitot pressures in the core are 
summarized in Table IV. 1 


Table IV. 1 


configuration 


He2 tunnel 
He2 nozzle 


He 3 tunnel 
He3 nozzle 


He 3 tunnel 
He4 nozzle 


pitot pressure (inches Hg) 2.4±0.08 


3.9+0.08 2.4+0.08 


Mach number 


16.2+0.1 13. 8±0. 1 16.2+0.1 


The He3 nozzle is 6 inches longer than the He2 and He4 nozzles. Hence the 
displacement thickness at its exit is larger which would tend to decrease 
the expansion. It is thus the nozzle's design which is responsible for the 
lower than expected Mach number. In the presence of seeding, the Mach num- 
ber may be found spectroscopically as discussed in Section IV. C. It is 
reduced from the above measured value due to the injection process (Section 
IV.B.l). 


3. Atomic Beam Device and Frequency Calibration 

The atomic beam spectra of the D 2 line of sodium is depicted in Figure 

IV. 10 for three different laser scan starting frequencies. Each of the 

larger and smaller peaks contain three hyperfine transitions associated 
2 2 

with the 3 F =2 a ' nc * 3 s ^ 2 ’ F=1 8 roun d state respectively. 

2 2 

The 3 Sjy 2 > F=2 ■+ 3 hyperfine transition is mainly responsi- 

ble for the larger peak in the spectra, since it is almost three times 
stronger than the next strongest hyperfine transition. The hyperfine lines 
in each peak are closely spaced in relation to the spectral purity of the 
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laser. The natural width of the transitions causes some overlap between 
them. These reasons in addition to the fact that only 10 sampling points 
were used to map the spike and the fact that a large amount of scattering 
existed in the atomic beam probe volume, are responsible for the unresolved 
structure of the line. 

The experimentally observed HWHM of the large peak is 10 MHz. The 


natural HWHM is 5 MHz (Appendix A) . If the large peak represents the 


2 2 

5 Sjy 2 > F=2 **■ 3 Pj/ 2 * P =3 hyperfine transition only, then the experimental 


contribution to the HWHM is 5 MHz. This implies that the laser HWHM is 
5 MHz if it is assumed that the laser linewidth may be simply added to 
the Lorentzian natural linewidth. 

2 

The three lines emanating from the 3 F=1 state are not resolved 

for the same reasons. Since the line strengths of the two transitions 


3“Sjy 2 , F=1 3 r 3 y 2 » F=1 > 2 (Figure A. 2) are comparable, one may not 

assign the peak to any transition in particular. Therefore, for velocity 


2 2 

measurement purposes it is the A S l/2* F 2 * A P 3/2 , F=3 transition which 


provides the best accuracy as a frequency reference. 

The dye laser frequency is tuned by applying a high voltage ramp to 
a piezoelectric crystal which translates one of the end mirrors of the 
cavity. The intracavity etalon installed to facilitate single frequency 
operation has to be tuned so as to track the cavity mode. Since the fre- 
quency-voltage relation is nonlinear, it has to be calibrated. Figure A. 10 
depicts the calibration procedure which is performed in conjunction with 


the atomic beam device. First, the minimum laser frequency is adjusted so 


2 2 

that the strongest transition (3 Sjy ? , F=2 •+ 3 P^ 9 ) appears just as the 


laser starts scanning. The laser is then tuned to a lower initial frequency 
so that the absorption line corresponding to the strongest transition falls 
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2 2 

on top of that of the second transition (3 F=1 + 3 P of the 

first scan. This is repeated a third time for the last frequency range, 
providing a piecewise calibration. Linear interpolation was used for 
frequencies falling between the four calibrated points. The determina- 
tion of these four points is again subject to the error sources mentioned 
in relation to the atomic beam device spectrum. For frequency calibration 
purposes the computer collects 2000 points instead of a mere 1000 for a 
regular run to enhance the resolution. 

A more accurate calibration scheme, which averages over the laser fre- 
quency fluctuations is presented in Section IV. D. 2. 


4. Laser Spectral Purity 

The jitter obtained with the Spectra Physics model 580 dye laser is 
within the resolution of the Jodon 7 MHz linewidth spectrum analyser. In 
the absence of amplitude noise and frequency jitter, a scanning interfero- 
meter would display a perfect bell-shaped curve. When the frequency jitter 
of the dye laser has a single component at a frequency which is lower than 
the sweep rate of the interferometer, several of these bell-shaped curves 
will appear on the oscilloscope screen, one for each spectrum analyser 
scan, as depicted in Figure IV. 11. a. When the dye laser has a jitter at 
a frequency which is higher than the interferometer scan rate, the inter- 
ferometer transmission varies rapidly with time as the dye laser frequency 
jitters back and forth over various portions of the bell-shaped response 
curve as shown in Figure IV.ll.b. 

The utilization of the spectrum analyser to measure dye laser jitter 
as outlined is not accurate. It was used however to obtain an upper limit 
on the laser frequency jitter. The maximum frequency excursions are about 
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V * • 


25 MHz and occur at a rate of about 5 Hz. The frequency excursions due 
to frequency jitter introduce an offset in the velocity measurement since 
the center of the absorption lineshape varies from scan to scan. A 25 MHz 
frequency excursion corresponds to a velocity measurement error of 15 m/sec 
at the sodium transition wavelength. More reliable information regard- 
ing the error is obtained in Section IV.D.l and IV. D. 2. 
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C. Data Analysis 

The data consisting of somewhat more than 600 points per laser scan 


were reduced by fitting with a Voigt profile, V(x). x is the detuning 
normalized by the Voigt HWHM, Av . 


(33) 


x = 


v-v 

c 

Av.. 


IV. 2 


where V q is the line center frequency. 

The fitting function was generated using an approximation devised by 
(33) 

Kielkopf v J which is accurate to within 0.0001 of the peak value of the 
function. The function is defined so that 

V (x = 0) = 1 IV. 3 

V (x = 1) = 1/2 IV. 4 

Each of the six hyperfine transitions (Section II. B, Appendix A) com- 
prising the D 2 line was frequency shifted and intensity weighted (Figure A. 2) 

to generate the complete theoretical spectrum. The best fit was determined 

(34) 

by a computer program incorporating a least squares routine ZXSSQ^ 
supplied by I.M.S.L.^^ and described in Appendix B. 

A minimum of five parameters, needed to fit a complete D ? spectrum are 
listed below. 

1) The ratio of Lorentzian over Gaussian HWHM, a (equation A. 108), deter- 
mines the shape of the curve. 

2) The Voigt HWHM determines the frequency scale. 

3) The frequency, F^, of the 3^S^ ? » F=2 -> F=1 transition deter- 

mines the position of the Doppler shifted spectrum along the frequency 
axis . 

4) The background intensity, I q , determines the base line. 

5) The intensity normalization factor, Ij, determines the intensity scale. 


V 


-63- 


Once an optimal value for the parameters has been established one 

may determine the Lorentzian (LHWHM) and Gaussian HJVHM (GHWHM) from the 

(331 

Voigt HWHM and a by using the relations 
£ = 0.099 

2 


l = 


IV. 5 


1 + e £n2 + J { l-ejln2) 2 + 


4Zn2 


a 


g = 


1 - (l+eto.2) ft * £Jln2 ft 


IV. 6 


ftn2 

LHWHM = ft • HWHM 

GHWHM = g HWHM /ftn2 IV. 8 

The collisional and natural broadening mechanisms are Lorentzian. Their 
combined linewidth is therefore a sum of the collisional (CHWHM) and natural 
(NHWHM) HWHM (Appendix A. ) . For simplicity it is assumed that the laser 
HWHM (IHWHM) may be added to the Lorentzian broadening component. This 
could be justified based on the analysis of the atomic beam device spectra 
(Section IV. B. 3). The result is 

CHWHM = LHWHM - NHWHM - IHWHM IV. 9 

The temperature, T, is given by 


. GHWHM, 2 
o'- 628 j 


where T is 273°K and GHWHM is in MHz. 
o 


A. 35 


If the large spike in the atomic beam spectrum of frequency FREFL is 
2 2 

due to the F = 2 + 3 transition, equation A. 29 and the data 

from Figure A. 2 may be used to find the velocity component, u, in the 
laser beam direction. 


u = [FI + (60+34)xi0 - FREFL] • X 


21 


IV. 10 


. '• v 
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where th e transition wavelength. The Mach number can be deter- 


mined from the relation 


(36) 


M = — = — — 


IV. 11 


a /YRT 
/ W 

where a is the velocity of sound, y is the ratio of specific heats for 

helium, R is the universal gas constant and, W is the atomic weight of 

(37) 

helium. The static pressure, P, is given by 

1 

y-1 


P = P t [( I ^-)M 2 ] y ~ 1 [ ] 


IV. 12 


2yM - (y-1) 

where is the pitot pressure* Since M is the result of a spectroscopic 

t 

measurement the value of P is based on local quantities (P , M) and not 
on the stagnation chamber pressure (Section IV. B. 2). Finally, a value 
for the pressure broadening cross section, o c , can be calculated. From 
equation A. 67 

•1 

-) [P. /rrskr] 

vf/T 


P/P -1 5 — 

■ c - CH "™ hd=’ t p o Jwn 


IV. 13 


where P and T are the static pressure and temperature respective, T q is 

273°K and P Q is 760 torr (l.OlxlO 5 Pa), k is the Boltzmann constant and 

y is the sodium helium reduced mass. 

The least squares fitting method, described in Appendix B, involves 

2 

the calculation of the sum of squares, y 

2 N yir f cv 2 

X “ = 2 — 3 —r 

j=i j 


B . 20 


where N is the number of sample points, y^j is the fluorescence intensity 

from the tunnel at a frequency x^ , f(x^) corresponds to the theoretical 

model and o. is the standard deviation. An estimate of a., the conse- 
J J 

quence of counting statistics and instrumental uncertainties, has to be 
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found. The fluctuations in the observations as a result of the counting 
statistics are not due to any imprecision in measuring the time interval 
or counting the number of events. Rather they are a result of the fact 
that a random sample of events distributed randomly in time contains num- 
bers of events which fluctuate from sample to sample. Physically this is 
due to the discrete nature of the photons in the laser and fluorescence 
electromagnetic field, the sodium atoms seeded in the flow and the photo- 
electrons produced in the photomultiplier tube. These fluctuations follow 

( 38 *) 

the Poisson distribution 1 J as discussed in Appendix B. 

The instrumental uncertainties are due mainly to the laser frequency 
jitter and to noise induced in the shielded cables, 25m long, connecting 

the experiment to the A/D converter (Section II. F). The sum of squares, 

2 ' 2 
X , may also be expressed in terms of the variance of the fit and 

~2 

the variance of the data a (Appendix B) . 



2 


B. 22 


where n is the number of parameters, a is a characteristic of the dis- 
persion of the data y^ in the parent population j and is not descriptive 
' 2 

of the fit. is a characteristic of the spread of the data about the 

2 2 

fit but is also an estimate of a . x is the ratio of the estimated var- 
iance to the parent variance, times N-n, which makes it a convenient mea- 
sure of the goodness of the fit. If the fitting function is a good esti- 

'2 ~~2 

mate of the parent function, should agree with a and 


X^ - N-n IV. 14 

The parent distribution standard deviations, o_. , have to be used in 
equation B.20 to be able to take advantage of the "Goodness of Fit" test. 
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The fitting programs are executed for two standard deviations. Using 

2 

0^ = 1 one clearly does not obtain X = N-n. Assuming that the model 

describes the physics and that the data are consistent with the model, 

2 2 
then x = N-n, if the value v = a = x /(N-n) is chosen. The para- 

meter uncertainties are then^°) 

a k = [(a) _1 ]J( 2 • a fk-l.,..,n) IV. IE 

where relates to parameter 3^. a is defined in Appendix B. For small 
residuals a may be approximated by 
N 3£ (x ) 8f Cx ) 

■ A “V - 

where the fitting function f is evaluated at the solution to the least 
squares problem as designated by the subscript o. This approximation is 
incorporated into the computer program of Appendix C to calculate the un- 
certainties in equation IV. 15. 

No counting was performed in the present experiment, though the photo- 
multiplier used to monitor the tunnel fluorescence is inherently a digital 
device. The output was passed through an integrator, and a voltage was 

measured. Using CL = /y „ , where y„ is the observed voltage, one again 
2 

cannot use the x "Goodness of Fit" test since the ratio of the number of 
counts to the output voltage is unknown. Making the same assumptions as 


for the choice ck = 1, then 

* t(5) _1 ]J, /2 • o.xr- 


(k = l,,..,n) 


where 


x ^- 

" N - n 
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and 


2 N y,,-f(x.) 2 

^ j=l ^ 

The residues are defined as 


IV. 18 


R. 

J 


. y id ( ?-£ 

0 . 

J 


IV. 19 


Each y. . belongs to a different parent population with mean f(x.). Once 
r 3 j 

[y_-f(Xj)] is divided by Oy all the residuals are encompassed into one 
parent population. The residuals are plotted versus for every fit. 

They would be randomly distributed around zero for a good fit. 

A histogram of the residuals was drawn for each fit. The range of the 
residuals was divided into a hundred slots. The histogram was fitted with 
a Gaussian using three parameters. 

1) mean of histogram 

2) standard deviation of histogram 

3) probability normalization factor. 

A good fit should have a histogram which follows the Gaussian envelope 
closely. ^0) 

An additional criteria for the consistency of the data is that the 
observations be randomly distributed within the histogram. The data points 
are plotted in sequential order, versus their histogram frequency of occur- 
rence. 

Apart from using two different standard deviations, the fits were 
performed for the following cases: 

1) Using all the generated data points. 

2) Neglecting the points surrounding the smaller peak corresponding to 

2 

transitions involving the 3 F=1 state. 
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3) With one set of intensity normalization factor, Voigt HWHM and ratio 
of Lorentizian over Gaussian HWHM. 

4) With two sets as in 3 , one for each peak in the line. 

The errors in the quantities calculated through the parameters are 


obtained by using 


(41) 


°x = * ( lr )2 A3 k 

X k=1 33 k k 


IV. 20 


2 

where x is the quantity of interest, 3 k is the parameter, A(3 k is the var- 
iance of the parameter and n is an integer signifying the number of para- 
meters. Equation IV. 20 holds if the parameters are uncorrelated. For 
errors including a systematic component equation IV. 20 yields too large a 
standard deviation. Systematic errors will occur if optical pumping, sat- 
uration effects, absorption effects or excessive laser frequency jitter are 


present. 
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D. Experiment 

1. Free Stream Flow Investigation 

The Doppler shifted spectrum obtained by shining the laser into the 
seeded flow and collecting the fluorescence signal from the probe volume 
is depicted in Figure IV. 12. There are somewhat more than 600 data points 
represented by the crosses. The continuous line is a fitted Voigt function 
discussed in section IV. C/ and Appendix A. The two noisy spectra on the 
left are due to the unshifted sodium line obtained from the atomic beam 
device. The three vertical lines inside each peak corresponding to the 
tunnel spectrum represent the theoretical line strength and location along 
the frequency axis of each hyperfine transition as determined by the fitting 
routine. 

The Doppler shift furnishes the velocity component along the laser beam 
direction. For Mach number calculations, the total velocity vector has to 
be determined. A velocity vector in a plane parallel to the tunnel axis 
was measured by using two different laser beam angles. 

The bottom of figure IV. 12 includes a plot of the residuals. Figures 
IV. 13 and IV. 14 depict the histogram and the frequency of occurrence of 
the observation points within the histogram respectively. It was established 
that using = v^yT is preferred over ck = 1 (Section IV. C), since then 
the histogram follows the Gaussian envelope more closely and the distribu- 
tion of data points within the histogram is more random. These two tests 
are not very sensitive since they do not indicate the presence of systematic 
errors very well. By fitting one peak out of the two in Figure IV. 12 and 
comparing the other with that predicted by the fit one can determine whether 
optical pumping, absorption effects, saturation or laser frequency jitter 
have affected the measurement. 



Figure IV. 12 
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The larger peak has been chosen for the fit since it contains more 
points, particularly in the low frequency tail. The tail of the spectrum 
is affected mostly by the Lorentzian component of the HWHM as discussed in 
Appendix A. Absorption effects tend to raise the intensity of the experi- 
mental observations of the smaller peak. Optical pumping tends to lower 
the intensity of the experimental observations of the smaller peak, a 
similar effect is observed due to the laser frequency jitter. Laser fre- 
quency jitter also tends to distort the lineshape. 

The quality of the data does not have much of an affect on velocity 
and Voigt HWHM measurement, but has a strong affect on the ratio of Loren- 
tzian to Gaussian linewidths. These three quantities, the velocity, the 
Voigt HIVHM and the ratio of Lorentzian to Gaussian linewidths, remain con- 
stant within the experimental accuracy throughout the hypersonic core in 
the tunnel. By applying the single peak fit criterion of "Goodness of Fit' 
one obtains the values based on four fits summarized in Table IV. 2. 

The velocity uncertainty is based on the laser frequency jitter as 
determined from the atomic beam device spectrum. The uncertainty of ±16 
m/sec is an absolute error and does not depend on the value of velocity 
measured. 

Time of flight broadening is neglected as it contributes less than 
0.5 MHz to the HWHM for a Gaussian beam waist radius of 0.001 m. (equation 
A. 55). The scaling of the collision cross section, a , with temperature 
may be determined by the expression 

Aw c = a c n G A. 66 

where Aw £ is the collision HWHM, n is the number density and u is the rel- 
ative atomic velocity. The collision HWHM is clearly a function of two 


if 
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Table IV. 2 

Voigt HIVHM (Half Width at Half Maximum) 

Ratio of Lorentzian over Gaussian HIVHM 

Collisional HWHM 

Gaussian HWHM 

Temperature 

Pressure 

Mach Number 

Collisional Cross Section 
Velocity Component Along Laser Beam 
Velocity Vector 


141+2 MHz 
0.43±0. 03 
48± 3 MHz 
107+3 MHz 
8+0. 4°K 

0.5610.03 torr (75+4 Pa) 
11+0.3 
(187114) 10“ 6 cm 2 
1265116 m/sec 
1749128 m/sec 
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thermodynamic variables; the density and the temperature, T, since the 
relative velocity, u, is a function of temperature 



A. 64 


where k is the Boltzmann constant and y is the reduced mass. Equation A. 66 
may be expressed in terms of the collision WVHM over the density. 

Aw 

— — = a u = KT a IV. 21 

n c 

where K and a are constants. This equation was fitted to two data sets. 

One obtained during this experimental investigation at 8°K (Table IV. 2) 

f 251 

and the other obtained by McCartan and Farr v 1 at 415°K. The value ob- 
tained for a is 0.38±0.03 which is lower than a = 0.41 as calculated by 
Lewis and McNamara ' (Appendix A). 


2. Double Beam Experiment 

The spectrum obtained by shining the laser into the tunnel with a mir- 
ror reflecting the beam back on itself is shown in Figure IV. 15. The first 
and second peaks are the Doppler down shifted manifold. The third and 
fourth peaks are the Doppler upshifted D£ manifold. The downshifted peaks 
are smaller in intensity since they are due to the reflected laser beam 
which has undergone more transmission and reflection losses. 

The same computer fitting routine is used as that for the single beam 

experiment with two additional parameters. One is the frequency, F2, of the 
2 2 

3 $ y / 2 ’ " P^/ 2 > F = 1 transition and the other is the intensity nor- 

malization factor, 12, corresponding to the downshifted laser beam. 

The ratio of the Lorentzian to Gaussian width and the Voigt HWHM are 
assumed equal for all transitions. The different variables are obtained 
in the same manner as before. The velocity may also be calculated by 



• i gurc 
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halving the Doppler shift between any two transitions connecting the same 
hyperfine levels. The values obtained by deconvoluting the spectrum are 
essentially the same as for the single beam experiment of the last section 
apart from the velocity sensitivity which is appreciably improved. The 
improvement in the velocity measurement accuracy is due to a more accurate 
frequency calibration scheme. Calibration is established in this experi- 
ment by assuming a frequency (v) versus ramp voltage (v) relation of the 
form 


V = CjV 2 + c 2 v IV. 22 

The constant c^ and c 2 are additional parameters determined by the least 
squares routine. Figure IV. 16 depicts the calibration method performed 
in conjunction with the fitting process. For this presentation only one 
transition per peak is assumed. The top graph is the spectrum intensity 
versus ramp voltage. The bottom graph is the laser frequency versus the 
ramp voltage. The frequency separation in the manifold, Av, is known 
from the literature. One may therefore write 


Av = C 1 (v 2 -Vj) + c 2 C v 2 -v 1 ) 


IV. 23 


2 2 

Av = c 1 (v 4 -v 3 ) + c 2 (v 4 -v 3 ) 

C l* C 2* v i» V 2’ v 3 anc * v 4 are a *l determined by the least squares routine. 
In the present experimental work all the hyperfine transitions in the D 0 
manifold are utilized for the calibration. Nine instead of just seven 
parameters are used the two additional ones being c^ and c 7 . The effect 
of laser frequency drift is also minimized by the calibration scheme. This 
calibration is employed in the interpretation of all the spectra. 
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Since the laser frequency jitter is associated with a time which is 
much shorter than the one second scan time through one of the broadened 
peaks, this frequency jitter is averaged. The averaging process is an 
additional reason for a better velocity measurement accuracy. The uncer- 
tainty obtained is ±2 m/sec, an absolute error which does not depend on 
the value of the measured velocity. The total temperature, T*", of the 
flow is calculated using the expression^ 

2 

T r = + T IV. 24 

P 

where u is the flow velocity, c^ is the specific heat at constant pressure 
and T is the static temperature. The spectroscopically calculated values 
are depicted in Figure IV. 17, by the crosses which signify error bars. The 

t 

ordinate represents the total temperature, T , and the abscissa is the dis- 
tance from the center line. The errors are mainly due to an uncertainty 
in the measurement of the angle between the laser beam and the wind tun- 
nel axis. The five crosses with a black dot in their centers are measure- 
ments taken on the center line during five, six minute duration, wind tun- 
nel runs. Each corresponds t:o an average taken over three to seven laser 
scans. The five are separated along the abscissa for graphical presenta- 
tion only. The solid line represents the thermocouple total temperature 
measurement. The shaded area is the error associated with that measure- 
ment. It is due to an uncertainty of ±2°K and ±0.008 in the value of the 
adiabatic wall temperature and the recovery ratio respectively. The total 
temperature and its uncertainty based on the thermocouple measurement are 
not sensitive to the value of the static temperature. The static tempera- 
ture obtained from isentropic calculations (Section II. B) was used. The 
discrepancy between the crosses and the solid line are due to a systematic 


INTENSITY 



VOLTAGE 






error in the measurement of the angle between the laser beam and the wind 
tunnel axis of symmetry. 

3, Cone Flow Field Investigation 

A 10° half angle cone was mounted in the test section, with the apex 
0.067±0.001 m upstream of the measuring station. The pitot probe survey at 
the measuring station is shown in Figure IV. 18. The abscissa is the dis- 
tance from the wall and the last point is the cone surface. The probe vol- 
ume was scanned across the tunnel diameter from the inviscid region into the 
shock and boundary layer of the cone. The temperature and velocity of the 
flow were calculated from the spectroscopic data as in Section IV.D.l. 

The sodium helium collision cross section obtained from the free stream 
flow conditions was used to spectroscopically determine the pressure in the 
vicinity of the shock and inside the boundary layer. This is a good approx- 
imation since the main temperature dependence of the pressure broadening 
mechanism is due to relative velocity and density effects and not to the 
cross section as discussed in Appendix A. The results are presented in 
Table IV. D. 3. 

The uncertainties in the relative position of the points are 0.003 
inches. The position of the cone surface is known to within 0.01 inch. 

The static pressure in the free stream based on the pitot survey is calcu- 
lated using the spectroscopically determined Mach number. The velocity 
vector in the free stream is inclined at an angle, a, of 43°40' from the 
laser beam. The velocity vector downstream of the shock is calculated 
assuming a (maximum) stream line turning angle of 10°. This changes a to 
44°30* (maximum). The pressure and temperature on the cone surface from 
compressible flow relations are 4.1 torr (547 Pa) and 27. 7° K respectively. 




Table IV. 2 
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The calculation is based on the free stream Mach number obtained from the 
pitot surveys. 

The pressure, P, temperature T, velocity component, u^, velocity vector, 
u and Mach number, M, are depicted in Figures IV. 19 and IV. 20 along the 
ordinate. The solid line is the pitot pressure survey drawn for position 
comparison while the abscissa is the distance from the cone's surface. 

The data near the cone's surface is less reliable due to strong light scat- 
tering and a weakening of the fluorescence signal which may be due to a 
depletion in the number of sodium atoms. The source of this drop in signal 
intensity is unknown, but may be caused by loss of sodium to the cone's sur- 
face . 

For flow visualization purposes the laser beam is expanded into a 
sheet of light by using two cylindrical lenses as shown in Figure II. 6 
The normal to the light sheet lies in the plane of the figure. A 10° 
half angle cone was mounted in the test section. The resulting picture is 
shown in figure IV. 21. Since the shock forms a cone around the model, its 
intersection with the sheet of light produces an ellipse. The long axis of 
this ellipse is hidden by the model. The fluorescence is more intense 
downstream of the shock since the density is larger and the laser is tuned 
to a frequency which highlights the particular velocity existing behind the 
shock. The shock angle measured by drawing a line from the cone tip tan- 
gent to the fluorescing ellipse is 13°50'±10'. This corresponds to a 
f311 

Mach number of 7.34±0.34. The Mach number obtained from the spectro- 
scopic calculation is 11±0.3. The discrepancy is related to hypersonic 
interaction. Si^ce the shock angle in hypersonic flows is very shallow, 
the shock interacts with the boundary layer near the tip of the cone. The 
cone appears to the flow with a virtual half angle which is larger than the 
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real half angle. The interaction weakens further downstream since the 
shock moves into the inviscid flow. The shock is therefore curved es- 
pecially near the tip and one may not calculate its angle as outlined. 


Chapter V 

CONCLUSIONS AND SUGGESTIONS FOR FURTHER STUDY 


A. Conclusions 

1. Introduction 

As a result of the experimental work described in this thesis, the 
Resonant Doppler Velocimeter (RDV) technique has been established as a 
means of measuring average components of temperature and velocity in com- 
pressible helium flows. Furthermore, the sodium helium collision cross 
section at low temperature was determined using the spectroscopically 
measured temperature and collision line broadening coupled with a pitot 
tube measurement of the pressure in the free stream. 

Using the cross section data, the RDV technique is capable of simul- 
taneously measuring velocity, temperature and pressure. This was demon- 
strated by mapping the flow field in the vinicity of a conical model. 

2. Doppler Shift Uncertainty 

The Doppler shift and the resulting velocity uncertainties are sum- 
marized in Table V.l. 

Table V.l 

single beam double beam 

Doppler shift uncertaintiy (MHz) 28 3 

Velocity uncertainty (m/sec) 16 2 

The experimental uncertainties are due to counting statistics, instru- 
mental uncertainties and systematic errors. An inspection of the Doppler 
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shift data indicates that the position of the spectra (from the atomic 
beam device and the wind tunnel) along the frequency axis and not their 
shape is responsible for the uncertainty in measurement. The limiting 
factor influencing the position of the spectra is the lasfci’ frequency 
jitter and drift. The error in velocity measurements is thus absolute. 
Hence, the larger the measured velocity component, the smaller the rela- 
tive error. 

No averaging of the frequency jitter was performed with the atomic 
beam spectra since the jitter was larger than twice the natural HWHM of 
sodium. Frequency jitter averaging was possible with the tunnel spectra 
(through the use of the data reduction scheme with the Least Square 
routine) since twice the Voigt HWHM is an order of magnitude larger than 
the laser frequency jitter. The large velocity uncertainty associated with 
the single beam experiment reflects the lack of averaging with the atomic 
beam device spectra. 

Were the laser frequency jitter and drift and systematic errors to be 
eliminated, then the sensitivity for average velocity measurements would be 
limited by counting statistics and turbulence. The effect of turbulence 
on the average measurements depends on the relative value of the charac- 
teristic turbulence time and the time it takes the laser to scan through 
the spectrum. If these times are of the same order, then the turbulence 
would be manifested through an erroneous linewidth measurement. For a 
scanning time shorter than the turbulence time, scan to scan inconsisten- 
cies would occur in the velocity measurement though the linewidths would 
not be appreciable effected. Random noise is added to the counting statis- 
tics if the scanning time is longer than the turbulence time. The scanning 
time has to be adjusted for the latter possibility since only then can aver- 
aging be performed. Longer scans also improve counting statistics. 


With an actively stabilized laser incorporating a feedback loop one 
can easily obtain jitter which is smaller than 1 MHz. This is smaller than 
the natural HWHM of sodium. Since the natural lineshape is a Lorentzian 
it may be fitted using a Least Squares routine to find the position of all 
the hyperfine transition frequencies within the sodium D 9 manifold. A 
lower than 0.59 m/sec velocity uncertainty can be expected using a laser 
stabilized to 1 MHz. This is due to the fact that frequency jitter averag- 
ing is then facilitated not only with the wind tunnel spectra but also with 
the atomic beam spectra through the fitting process. 

3. Linewidth Uncertainties 

The temperature, cross section and pressure sensitivities due to line- 
width measurement uncertainties are summarized in Table V.2. 

Table V.2 

Temperature Uncertainty 5°o 

Cross Section Uncertainty 9% 

Pressure Uncertainty 10 ? o 

The quantities in Table V.2 are calculated from the Lorentzian and Gaussian 
linewidths (Section IV. C). These in turn are derived from the Voigt HWHM 
and the parameter which is defined (Appendix A) as the ratio of the Lorent- 
zian over the Gaussian widths. An inspection of the Lorentzian and Gaus- 
sian linewidth data indicates that their measurement sensitivity was limited 
by the uncertainty in the parameter a and not the Voigt HWHM, The value of 
a is influenced appreciably by the presence of systematic errors in the 
spectrum (Section IV. C) and by a low signal to noise ratio (due to counting 
statistics) . 


T~ V 
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4. Scan Limitations 

The tunnel operating time of six minutes was too short to allow for 
the optimization of the power supply settings (to maximize the signal to 
noise ratio or to suppress absorption effects) . It takes the temperatures 
of the different components in the seeding device time to equilibrate. 

Also, the power supply settings vary frpm inn to run since the quality of 
helium deteriorates with increasing numbers of purification cycles. Hence 
only a few laser frequency scans per run were taken. 

The tunnel fluorescence intensity was fed to the GMAD/1 which samples 
during a time smaller than 10 nsec at 330 Hz. The small value of the duty 
cycle raises severe signal to noise ratio problems. An integrator with a 
time constant of the order of the reciprocal of the sampling frequency was 
used in order to alleviate this difficulty. Too large a time constant in 
some of the laser scans gave rise to spectrum distortions which also caused 
scan to scan inconsistencies in the value of a. 

An improvement can be obtained by scanning the laser in discrete steps 
and using a photon counting device to collect the signal for a preset time 
at each frequency. The information stored in the photon counter would be 
transferred to the computer before the laser shifts to a new frequency. 

This would eliminate the integrator and the systematic errors associated 
with its time constant. Since counting is performed, the "Goodness of Fit" 
criterion (Appendix B and Section IV.C) may be applied to the data. 

Photon counting enhances the signal to noise ratio allowing a reduc- 
tion in the seeded sodium number densities and laser power. A lower sodium 
density and laser power would suppress the systematic errors due to absorp- 
tion effects, saturation and power broadening. 
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5. Flow Perturbation 

About 850 watts are injected into the main helium flow in the stagna- 
tion chamber with the purge helium flowing through the sodium seeding device. 
The mass ratio of the two helium flows is less than 0.01. The effect of 
mass and energy injection into the flow is discussed in this section. 

No disturbance to the flow due to sodium injection manifests itself 
in the pitot surveys (Section IV. B. 2). On the contrary, the diameter of 
the core region seems to enlarge as a result of the presence of the injec- 
tion needle. The needle tip was positioned 0.1 inch upstream of the throat. 
The improvement in the core area may be due to the main helium flow under- 
going two expansions while flowing past the needle tip. The first expan- 
sion occurs at the base and the second at the tip of the truncated cone 
(part 2, Figure II.C.l). The resulting expansion waves could interact 
favorably with the compression waves downstream of the throat thus enlarg- 
ing the core region. 

Mach number information cannot be derived from the pitot pressure 
scans obtained while sodium injection was taking place since the flow then 
is not isentropic (Section IV. B. 2). An additional cause resulting in stag- 
nation pressure loss is the presence of compression waves. The expansion 
waves emanating from the needle tip are reflected off the nozzle wall. 

These reflected waves could turn into compression waves if the curvature 
of the nozzle wall is concave enough. 

The Mach numbers obtained from the spectroscopic measurements (Section 
IV.D.l) and from the interpretation of the cone flow picture (Section IV. D. 2) 


are summarized in Table V.3. 


Table V.3 


Spectroscopic measurement 11±0.3 

Shock angle measurement 7.34±0.34 

The shock angle measurement as outlined in Section IV. D. 3 is incorrect 
due to the hypersonic interaction effect. An accurate measurement may be 
obtained by illuminating the complete flow with an expanded collimated 
laser beam instead of a sheet of light (Section IV. D. 3). The shock angle 
should then be measured far enough downstream from the cone tip where the 
hypersonic interaction is negligible. This picture was not taken because 
of helium supply limitations. A picture of a shock surrounding a pitot 
probe is presented in Figure V.l . No quantitative results can be ob- 
tained from that picture since the shape of the pitot probe is not a sim- 
ple one. 

Figure IV. 17 indicates good agreement between the total temperature 
thermocouple measurement and the total temperature calculated from the 
spectroscopic measurements. 

6. The Utility of the Technique 

The technique has been demonstrated in a hypersonic wind tunnel since 
the experimental uncertainty in the velocity measurement, limited by laser 
frequency jitter (Section V.A.2) is absolute. Hence the higher the mea- 
sured velocity the lower the relative error. 

A meaningful deconvolution of the Voigt profile to obtain both the 
temperature and the pressure can be performed only for a parameter a on the 
order of one (Appendix A) . This requirement dictates the preferred range 
of values for the temperature and pressure in the test section of the wind 


tunnel. 



Figure 




Helium was used as the flowing gas since it neither reacts with nor 
quenches sodium. 

Q- 

As discussed in Section V.A.2, accuracies of the order of centimeters 
per second can be expected with a 1 MHz linewidth laser. Hence the mea- 
surement of velocities much lower than those encountered in this experi- 
mental investigation would be possible while still keeping the relative 
error down. The experimental facility is currently undergoing modification 
in an effort to utilize the technique in a Mach 2.9 nitrogen flow, 
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B . S uggestions for Further Study 

1 . Simultaneous Mult !<joint Data Collection 

The spatial resolution in the flow field measurements are due to the 
lens aperture combination mounted in front of the photomultiplier tube and 
the focusing of the laser beam, The flow field properties in a plane can 
be measured by illuminating a plane with an expanded laser beam and imaging 
into a photodiode array. The size of the plane is determined by the magni- 
fication of the collection optics and the area of the array. The spatial 
•resolution will be limited by the matrix density of the array and the number 
of channels which the data acquisition system is able to handle. 

2 . Multiple Velocity Components 0 

The Doppler shift facilitates the measurement of the velocity component 
along the laser beam. Three laser beams intersecting at different angles in 
the flow will be needed to obtain the velocity vector. A chopper can be 
used to block all but one beam at a time. The fluorescence intensity 
picked up by the photomultiplier tube may be fed into the computer and gated 
by the chopper to separate the components. Using a chopper wheel with two 
holes along a diameter, one may split the laser beam into three components 
designated in Figure V.2 by the black dots. The dimensions of the apertures 
should be such that only one beam at a time intersects the seeded flow. 

G 

3. Elimination of the Doppler Effect 

Two photon spectroscopy by two counter propagating laser beams has been 

{' 43 ') 

used to eliminate Doppler broadening. Circularly polarized light has 

to be employed to suppress the Doppler background from single beam two 
photon absorption if this type of experiment is performed in a cell. 
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By applying this technique to spectroscopy in the tunnel, one may elim- 
inate not only the Doppler broadening component but also the Doppler shift 
due to the convective velocity of the atoms. The use of linearly polarized 
light causes two Doppler shifted and broadened spectra to appear on both 
sides of the Doppler free absorption peak in the wind tunnel. Each is due 
to two photon absorption by a single linearly polarized beam, and, as com- 
pared to the spectrum in a simple cell, both of these spectra are well 
removed from the Doppler free peak. 

The velocity may be calculated using as a reference frequency the un- 
shifted absorption peak as in the single beam experiment (Section IV.D.l). 

The velocity may also be calculated by dividing the shift between the two 
Doppler broadened spectra into two, to obtain the Doppler shift as in the 
double beam experiment. The Lorentzian linewidth component is provided by 
the Doppler free absorption peak. It can be compared to the deconvoluted 
value obtained from the two side spectra from which the temperature can be 
determined. 

The technique could be useful in turbulence measurements since the pres- 
sure and temperature broadening mechanisms can be decoupled. 

An additional advantage of two photon spectroscopy is that the excita- 
tion and fluorescence frequency are different. Laser scattering can thus 
be strongly filtered facilitating measurements near surfaces. 

The sodium two photon transition with the largest cross section which 

is accessible with the R6G dye laser is 3S ■+• 4D. The total cross section 

-19 "> 

integrated over the frequency response of the atom is = 2.8x10 cm“ 

MHz (W/cm^)"*.^^ For an irradiance of 15 W/m^ used in the calculations 

-22 2 

of Chapter III one obtains Op = 4.5x10 cm MHz. The total cross section 

-7 2 

for the sodium D 2 line is 1.1x10 cm MHz which is fifteen orders of 
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magnitude larger. 

To improve the signal strength one could seed the flow with more atoms. 
This is possible due to the fact that trapping is less of a problem for the 
3s ■> 4p transition (3302 A) than for the D lines since the absorption 
oscillator strength, f, is only^ 4 5 "* 0.0152 for the 3S -*■ 4P transition. 

The laser scanning time in this experiment was about 7.3 seconds. Longer 
scanning times would improve the statistical averaging procedure. The laser 
power used for a 1 mm Gaussian beam radius was 2.4x10”^ IV. The power avail- 
able from the dye laser used in this experimental investigation for a 1 V« 
argon laser pump beam is 50 mW. With this power and with tighter focusing 
the light intensity could easily be raised much higher, strongly increasing 
the two photon absorption cross section. 

4 . Accurate Frequency Calibration 

A hot low pressure sodium vapor cell may be used in place of the atomic 
beam device for a frequency reference. Doppler broadening would predominate 
in such a cell. A temperature of 373°K results in a Gaussian HWHM of 734 
MHz (equation A. 31). The laser frequency jitter would be averaged (Section 
V.A.2) since it is an order of magnitude smaller (Section IV. B. 4) than 
twice the Doppler linewidth in the cell. The calibration of the frequency 
scale would become possible in a manner similar to the one employed in the 
double beam experiment (Section IV. D. 2) by using the Least Square procedure. 
A velocity measurement sensitivity of the order of 2 m/sec would be at- 
tainable using the suggested experimental setup. Saturation spectroscopy 
in a sodium cell may also be used as a frequency reference. The advantage 
to be gained over ordinary absorption spectroscopy in a cell is that the 
linewidths are adjustable over a wider range. The disadvantage is that 
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the spectrum is complicated by cross over terms. Saturation spectroscopy 
of the D2 line furnishes twelve instead of just six peaks. 

A frequency reference cell and its related detection optics is not 
needed in a double beam experiment. In this case the laser beam has to be 
reflected back on itself, an experimental requirement which it is not always 
possible to fulfill. This constraint is removed by employing a hot sodium 
cell as a frequency reference. 
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1. Introduction 

Sodium has an atomic number of 11, an atomic weight of 22.9898, a 
melting point of 97.81°C and a boiling point of 882. 9°C. Sodium is a 
soft bright, silvery metal which floats on water with a specific gravity 
of 0.97. It decomposes water with the evolution of hydrogen and the for- 
mation of NaOH. The flame, characteristic of this reaction in air, is 
due to the reaction of hydrogen with oxygen to form l^O. Sodium may or 
may not ignite spontaneously on water depending on the relative amount of 
oxide to metal surface area exposed to the water. It normally does not 

ignite spontaneously in dry air at temperatures below 115°C. Seven iso- 

23 

topes are recognized but only j^Na exists naturally. 

The electronic energy levels of sodium and optical pumping are dis- 
cussed in the next section. Broadening mechanisms are presented in 
Section VI. A. 3. The complete absorption line shape including natural 
temperature and pressure broadening is derived in Section VI. A. 4. The 
derivation also incorporates saturation and power broadening. This is 
followed by a short discussion of radiation trapping. 


2. Energy Levels of Sodium 

2 2 6 1 

The electronic configuration of sodium is Is 2s 2p 3s . Only 
the outer shell electron is involved in the transitions corresponding 
to the visible region of the spectrum. By treating sodium as a hydrogen- 
like atom with a coulomb potential one obtains energy levels which are 
dependent on the principal quantum number n = 1, 2, ... only. The orbi- 
tal angular momentum quantum number (A) degeneracy present in hydrogen is 
lifted by taking the screening effect of the outer shell electrons into 
account using the central field approximation. For a definite value of n, 
then ft = 0, l, ...n-1. Spectroscopic notation associates ft = 0, 1, 2, 

3, 4, 5, ... with the letters ft = s, p, d, f, g, h, ... . 

Corresponding to each angular momentum there is an associated 

magnetic moment y generated by the motion of the electrons around the 

ft 

nucleus. The electron itself has an intrinsic magnetic moment y^ by 
virtue of its charge and spin angular momentum sfi. As a result of the 
electron's motion in the central electronic field it experiences a mag- 
netic field $ proportional to its orbital angular momentum 111. The 

X / 

perturbation Hamiltonian is given by 

■ -5. • V ■» • * 

where a is a constant which may be calculated if the central field 
potential is known. 

Capital letters are used to denote total angular momenta for atoms. 
With this notation 

E ls = a§ . t A. 2 

The angular momenta Lfi and are coupled by this interaction. To eval- 
uate the energy shifts, linear combinations of the wave functions |y, L, 
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S, M^, Mg> are taken to form a new zeroth order wavefunction |y, L, S, 
J, M,>. The total angular momentum quantum number, J, corresponding to 
the total angular momentum operator 


-V -V 

J = L + S 


A. 3. 


and the quantum number Mj corresponding to projection of the operator 
J are both good quantum numbers for the new set of wavefunctions. y 
represents all other quantum numbers needed to specify the configuration. 
The values of the quantum number J are restricted to 

|L-S|£J<L + S A. 4. 

by the usual angular momenta sum rule. As S = 1/2 for a single electron 
atom, one clearly obtains two values of J for each L / 0. This inter- 
action together with the relativistic correction to the kinetic energy 
*2 2 

operator -jjj V is responsible for the fine structure of the alkali 

spectra. Specifically each orbital angular momentum quantum number, apart 

from the ground state has associated with it two energy levels. For 

example if L = 1 then J = 1/2, 3/2. The usual spectroscopic notation for 

the atomic fine structure energy levels is Lj. Lfi, Sfi and Jfi represent 

the resultant of all electronic orbital angular momentum spin angular 

momentum s .fi and total angular momentum i of the individual electrons. 

1 1 

The values of |S| , j L | , j j| are obtained by the usual angular momenta sum 
rules. The expression 2|s|+l is known as the multiplicity and represents 
the number of fine structure energy level splittings. Figure A.l. depicts 
the energy levels of Na. 

Although the spin orbit interaction accounts for the primary splitting 


of the energy levels, additional effects exist which modify the energy 
levels still further. In the discussion so far the only property of the 
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nucleus used was its charge Ze where e is the elementary charge and Z 
is the atomic number, There are three additional properties of the 
nucleus which affect the electronic energy levels. First, the energy 
levels depend on the mass of the nucleus and hence on the specific isotope 
being used. Second, the nucleus possesses an intrinsic angular momentum 
th. The magnetic moment associated with f interacts with ^ the magnetic 
field produced by the electrons at the nucleus. Third, the nuclear electric 
quadrapole moment interacts with the electronic electric field gradient 
produced at the nucleus. The last two effects lead to the hyperfine struc- 
ture and their corresponding perturbation Hamiltonian is given by 

B 

H HFS = A J^* J + 21 (21-1) J(2J-1) ^C?"* J) + - I(I+1)J(J+1)} 

A. 5. 


where the magnetic hyperfine structure constant, Aj, is determined experi- 
mentally. The electric quatrupole interaction constant, Bj, is defined 
by 

*2 
r) V 

A. 6, 


*2 

B , = eQ<— ■=■> 
J Sz 2 


A. 


where Q is the nuclear electric quadrupole moment and < — =■> is the average 

dz 

of the field gradient at the nucleus over the valence electron. Bj 

vanishes for 1=0, 1/2 as the nuclear charge distribution is spherically 

symmetric in these cases. Similarly B T is also zero if J = 0, 1/2 for 

J 2 

then the electronic charge density is spherically symmetric and < — ~> 

(49 ) 32 

vanishes. Schneider et. al. foung Aj = 18.9 ± 0,3 MHz and 

Bj = 2.4 ± 0.3 MHz. The angular momenta Ifi and Jfi are coupled by this 

interaction. To evaluate the energy shift linear combinations of the 

functions | y, J, I, Mj, Mj> are taken to form a new zeroth order wave- 

functions |y, J, I, F, Mp>. The total angular momentum quantum number 
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corresponding to the total angular momentum operator 


->-v-v 

F = J + I 


A. 7. 


and the quantum number M p corresponding to the projection of the operator 
F are both good quantum numbers for the new set of wavefunctions. The 
values of the quantum number F are restricted to 

l J - I ll F l J+I A. 8. 

2 2 2 

The hyperfine structure of the sodium 3 Sjy 2 > 3 ^\/2 anc * 3 ^3/2 states 
are shown in Figure A.2.^ 3 ^. The numbers on the arrows representing the 
allowed transitions correspond to the theoretical line strengths. 

The total angular momentum quantum number, F, dipole selection rules 
are AF = 0,±1. These give rise to a phenomenon known as Optical Pumping. 
Exciting the transition F = 1 ■+ 3“P F = !» allows a decay 

channel into the 3 2 S 1 ^ 2 F = 2 level, for then AF = 1. This process removes 

2 

the atom from the 3 F = 1 state so that it will not be available 

2 

for another pumping cycle from 3 F = 1. 

2 2 

Tuning the laser frequency to the 3 S l/2* F = 2 * 3 P 3/2 * F = 3 
transition eliminates this problem, as then the upper level can only decay 
to the lower level by AF = -1. The magnetic quantum number, M p , dipole 
selection rules for spontaneous emission are AM p = -1, 0, +1. For absorp- 
tion and stimulated emission the selection rules are AM p = 0 if AF ^ 0 
for linearly polarized light, AM p = -1 for right circularly polarized 

light and AM p = +1 for left circularly polarized light. The transition 
2 ? 

3 S l/2’ F = 2 -*■ F = 1 will optically pump atoms into the energy 

2 

levels 3 S, F. = 2, M c = ±2 if linearly polarized light is used. 

1/ 1 1 

This is depicted in Figure A. 3. v^ere the heavy and thin lines correspond 
to absorption and spontaneous emission respectively. From the figure 
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it is clear that once the atom has undergone one excitation cycle and has 
ended up in the Fj = 2, Mp = ±2 levels, it is no longer available for another 
cycle since the selection rule for stimulated absorption for linearly 
polarized light is AMp = 0. 

If the transit time of an atom through the excitation volume defined 
by the narrow linewidth laser beam is much longer than the upper state 
lifetime, the fluorescence signal will be considerably reduce. This is 
due to the fact that the atom while in the laser beam will undergo several 
excitation cycles with a high probability that it will finally end up in 
a trap level. 

2 

If a linearly polarized laser is tuned to the 3 F = 2 ■> 

2 2 
3 P » F = 3 transition, levels 3 Pj^* P = 3, Mp = ±3 will not be 

excited, but no trapping of atoms in the ground state will take place. 

A simple analysis will show that apart from this one transition, all 

others in the two sodium D lines trap atoms due to optical pumping. 
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3. Broadening Effects 

a. Natural Broadening 

One may show^ 50 ^ that the expectation value of the dipole moment 
<er> of a two level atomic system behaves very much like a classical oscil- 
lator. With this in mind consider a classical oscillator consisting of an 
electron on a spring under the influence of an applied electromagnetic 
field. Since the electron performs a time dependent motion it acts like an 
antenna radiating an electromagnetic field and causing the dipole to lose 

energy, that is, the dipole is damped. The equation of motion is 

•? eE 

x(t) + 2T x(t) + or x(t) * — — cos tot A. 9. 


to Q corresponds to the transition frequency, to is the forcing frequency, 

E^is the electromagnetic field amplitude and T is the damping coefficient , 

i A. 10. 

3 c 


r = T 


2 2 
e to 

o 


3 

4tte 0 c m 

where m is the electron's mass and r 


.-15. 


2.82x10 m is the classical electron 

U 

radius . 

The average rate of work done by the field on the atom is given by 

2 


the power expressions 

, 27r/to 


TT(eEj‘ 


p * S75 / dt * eE o cos ‘ — 25T 


t/tt 


(to-to o ) 2 + (r) 2 


> o 


A. 11, 

This shows that the steady state forced oscillator absorbs radiation 
dissipating it in the damping process. For radiative damping this involves 
(following the classical treatment) scattering the incident radiation in 
a doughnut pattern characteristic of a dipole. The absorbed power is 
proportional to a Lorenzian lineshape which is a normalized function 
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defined by 


l(w-w r) ^ — 2 

0 (w-w ) 2 +(r) z 


A. 12. 


The half width at half maximum (HIVHM) for the sodium D2 line using 

r 


equation A. 10. is found to be Av = - 5.1 MHz. 

A quantum mechanical treatment was first obtained by IVeisskopf and 
f5' ? 1 

IVigner . It shows that the frequency distribution of the intensity, 
I^ifu), of a transition taking place between an upper level k and a low«r 
level i is 

T^/tt 

A. 13. 


- 'o 




where I Q is the total intensity of the line, is the transition frequency 
and T. . is the HIVHM of the Lorentzian lineshape. If A., is the quantum 

K 1 1 J 

mechanical transition rate from level i to j , then 

A. 14. 


r ki ■ ? \j *1 A ih k > j - 1 > h 


For a transition from an upper level k to a lower level i one defines the 


emission oscillator strength f^ as the ratio between the quantum mechani- 

(45) 


cal and classical rates 

A, • 

ki 

sr = " 3hg 


f ki 


2mw, . 
ki 


Z |<km. |r|im.> 


A. 15. 


k n. i m k 


where nij, is the quantum number associated with the degeneracy of level k 
and g^is the degeneracy of level k. For absorption one defines 

2i f ik * -*k f ki A - 16 ’ 


For the sodium D lines, i is the ground state and k the first excited 


state. Thus 

r 


12. = = |£ f =5 mHz 

2tt 2tt 2tt gj 01 

where F is the damping coefficient calculated in equation A. 10. and 


A. 17. 


,(45) 


f Q1 = 0.982^ , The quantum mechanical rate is close to the classical 


4 
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result since the oscillator strength, f Q j, is close to unity for the 
sodium D lines. The time uncertainty relation may not be derived using 
Heisenberg's rule only since time is not an operator. The equation of 
motion for the expectation value of the operator A is 

ifi ft <A> = < t A » H ^ > + ifi < !l > A * 18, 

Heisenberg's uncertainty relation is 

AA A' > <i [A , H] > | A. 19. 

where the root mean square deviation of A is defined by 

(AA) 2 = <i|>| (A-<A>) 2 |iJj> = <A 2 >-<A> 2 A. 20. 

If A does not depend explicitly on time one obtains from A. 18 and A. 19 
• AH >1 j A. 21. 

provided^ 535 


AA 

A = d<A>/dt 


A. 22 


appears as a time characteristic of the evolution of the statistical 
distribution of A. It is the time required for the center <A> of this 
distribution to be displaced by an amount equal to its "HWHM" AA. 

Kibble et.al.^ 54 ^ measured the lifetime of the upper state 3 2 P,,/ 2 of 
sodium to be x = (16.0 ± 0.5)10 sec. From the uncertainty relation one 


may obtain the HWHM 

AM AE 

Av * _ v ^ 
n 2irf» — 


4tit 


= 5 ± 0.16 MHz 


A. 25 


b. Doppler Broadening 

At any time each atom has some instantaneous velocity. This causes 
the atomic absorption frequency to be Doppler shifted to a new frequency. 

Consider an inertial coordinate frame K and a moving one K' . The phase 
of a planewave is an invariant quantity. This is because the elapsed phase 
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of a wave is proportional to the number of wave crests that have passed 
the observer. Since this is merely a counting operation, it must be 
independent of the coordinate frame. A planewave of frequency w and 
wave vector £ in the inertial frame K will in general have a different 
frequency w' and wave vector in the moving frame K* . The phase of the 
wave <j> is invariant. Hence 

4> = wt - k • x = w't' - lc* • x' A. 24. 

After substituting the Lorentz transformation into equation A. 24. one 

f51) 

obtains the Doppler shift formulas for light waves 
tan e> = sin 6 


y(cos0 - 3) 
a)' = yw(l - $cos0) 


A. 25. 


A. 26. 


where 0 is the direction of k relative to the velocity vector u and 

6 = A. 27. 

c 

y = (1 - 3 2 )' 1/2 A. 28. 

For nonrelativistic gas dynamics 3 « 1. Then 
0 ' =0 

w' = w(l - 3cos0) A. 29. 

The condition for the largest probability of absorption is that w' 
correspond to the transition frequency in the atom's rest frame. 

Assuming a Maxwellian velocity distribution function one obtains 
the absorbed irradiance, J, 

J = J 


(27) 


c / m r -m , c .2 

o vT / 2ikT exp ^fkT ( v7T- ) ( ik- 1 ^ 


A. 30. 


- -m . r . 7 .7. 

s • i r ^nj\i 2kl v.. 

ik ik 

where m is the weight of one atom, T is the temperature of the sample and 


J is the incident irradiance. The Doppler HIVHM is 

° v,, 


AV r 


ik / 2kT ln2 
c / m 


A. 31. 


A. 32 


Rewriting equation A. 30 in terms of the HWHM, one obtains 
J = J Q G(v-v. k ,Av D ) 

where the Gaussian is a normalized function defined by 


G(v-v ik ,Av D ) 


1 / ln2 
AVp / it 


e 


v-v 


-ln2 (- 


Av 


ik-j 2 

D 


A. 33 


Equation A. 31 may be used to determine the temperature as a function of 
the Doppler linewidth. 


T ■ »W 2 2FW 


A. 34 


Normalizing by the broadening HWHM, Av Do , at T q 
V >2 line 

AV D 2 AV D 2 

T ■ Vet -5 ■ V<air> 


= 273°K, then for the sodium 


A. 35 


Do 


for Av q in MHz. 


c. Time of Flight Broadening ^^ 

At low pressures an absorbing atom may interact coherently, without 
undergoing collisions, during its entire transit through the laser exci- 
tation region. As in the discussion of natural broadening, one obtains 
a time of flight limited HWHM which cannot be less than (4TTAt)" 1 where 
At is the duration of the perturbation. Consider a Gaussian beam with a 
waist w and for simplicity assume a plane wave behavior. The field may 
be represented by 

E(x,y,z,t) = £E e "^ r//w ^ cos(kz-wt) A. 36 

where £ is the polarization direction, r is the radial distance from the 
beam axis and z is the distance along the beam. 
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The plane wave assumption holds for particles under the influence of 
approximately uniform phase fronts. It is thus suitable for particle 
with a lifetime which is short compared to the transit time through the 
laser beam. 

Since each particle interacts with the laser beam independently of all 
the others one may translate the coordinates of each in time and rotate 
them in the x-y plane. Choosing u y = 0> the particle enters the laser 
excitation region with velocity u x at an impact parameter of y = b as 
illustrated in Figure A. 4. The time axis is translated so that at t = 0 
the particle lies on the y axis. The transverse dependence of the spatial 
field is transform to 


2 2 
b +(u x tj 


E[r(t)] 


= E e - ( - r/,w - ) = E e 


-cV -t ^) 2 

c Q w' „ W ' 

= E e e 

o 


This is the field the particle experiences in the zero pressure limit. 

The time dependence is the same for all particles and does not depend on 
the impact parameter b. In other words, each particle of a given velocity 
subgroup will have the same time of flight residual HWHM. 

In the absorber’s rest frame 

-cV -c ^ + *) 2 

E(b,z,u ,u ,t) - § E e W e W cos[kz-w(l + — )t] A. 38 

X Z O C 

The phase angle <j> has been introduced to represent the phase of the alter- 
nating optical field upon entry of a given particle characterized by b, u x 

and u . The number of optical cycles during a transit time is given by 
z 

^=^=3.10 8 A. 39 

where u is the convective gas velocity. As this number is large, one may 


assume the effect of <J> to be negligible. The Schrodinger equation has the 
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form 

» H - (« 0 +H ')il> A. 40 

where H q and H' are the field free and perturbing Hamiltonian respectively. 
Consider a two level atom with unperturbed eigenfunctions ^(r) and ^ 2 (r), 
Then 

+ ' iw it + -iw 2 t ^ 

^(r,t) = aj(t)e ^(r) + a 2 (t)e \p 2 ( r) A. 41 

where a^(t) is slowly varying compared to oscillating term. In the dipole 
approximation 

H' = er*E A. 42 

H' is antisymmetric and thus has only nondiagonal matrix elements 

aj - - i <l|H» |2> e ia)t a 2 A. 43 


The density matrix is defined by the products 





/l 1 2 


11 

P 12^ 


fl a i 1 

a l a 2*' 

P 21 

P 22 j 


i a 2 a I* 

U 2 | 2 

/ 


Then^ 

Pjl - + jr <1 |H* | 2> e lWt p 21 + CC 
P 22 = - I <1|H'|2> e iWt p 21 ♦ CC 
P 21 = +|-<1|H'|2> e iwt (P 22 -P n ) 


A. 44 


A. 45 


A. 46 
A. 47 
A. 48 


The wavelength of the electric field is long compared to atomic dimensions 
in the optical frequency range. Inasmuch as the electric field does not 
depend on the atomic coordinates (electric dipole approximation) , the 
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matrix element reduces to 

_(k ) 2 .(Hil ) 2 u 

<l|H'|2> = eE o <l|x|2>e w e w cos[kz-w(l+ ^-)t] 


A. 49 


assuming E q to be polarized along the x direction. The density matrix re- 
lates to a definite velocity class of absorbers. The absorbed power is 
given by 

P = fito Fp 22 A. 50 

As p 22 varies during transit through the laser beam it has to be averaged 
over time. Hence, 

00 dp. 


P ■ lioJF / 


22 


dt 


dt 


A. 51 


F is the rate at which the atoms are injected into the volume in velocity 
space. An inspection of equations A. 47, A. 48 and A. 49 shows that for small 
saturation (Using a lower laser intensity, the upper state population is 
negligible and p^-p^ - -1.) the absorbed power is roughly proportional to 

U y^ 2 2 

[exp-(-^j— ) ] . Instead of performing the density matrix integration, one 

could calculate the distribution of frequencies in the atomic rest frame, 

u t „ 
x 2 

resulting from the action of exp-( 0 on the sharp laser frequency in 

w 

the laboratory rest frame. 


V 2 

00 -2 (— ) . „ 

C(oj) - / e " e -lwt dt 


A. 52 


The above expression should be averaged over the velocity distribution. 
The average random velocity in the wind tunnel for the present experiment is 


/8kT , 

u = / = 156 m/sec 

/ TTm 

This value is less than 10% of the convective velocity of 1736 m/sec. 


A. 53 
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Hen ce , no averaging is performed and for a first approximation one may 
use 


1 f W60. 2 

a " 2 u x 
C(to) = — e 

J 2 u v 


A. 54 


The HWHM is^ 55) 


Av T “ 2tt ^ ln2 w 

For a beam waist of 0.1 mm and a velocity of 1736 m/sec then Av^, = 


A. 55 
3 MHz. 


d. Pressure Broadening and Shift 

The classical model represents an atomic radiator as an elastically 
bound oscillating electron. The problem of collisional broadening is to 
calculate the absorption of such an oscillator when it is disturbed by 
colliding with the other atoms. The equation of motion is 

x + 2Fx + w 2 x = — cos tot A. 9 

o m 

where T is related to the natural HWHM, to is the atomic resonance fre- 

o 

quency, e is the elementary charge, m is the electron mass and E and 
to are the driving field amplitude and frequency respectively. 

One assumes that the collisions do not change the amplitude of the 
oscillator, only its phase. The atom continues its driven motion after 

a collision with a completely random phase between 0 and 2 tt. The instances 
at which collisions take place are also randomly distributed with an aver- 
age time between collisions equal to T ? . By solving equation A. 9, one 
may then average over the random phases and collision times. The time 

averaged power absorbed by the oscillator from the driving electromagnetic 
(57) 


field is then 
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P = 


2tt 

vs l (et cos wt >* dt ■ ir * r) 


A. 56 


The contribution to the Lorentzian HWHM due to dephasing is given by 
the reciprocal of the average time between collisions T 2 . From the kinetic 
theory of gases, the collision frequency per unit volume, z between atom 
a and b is^^ 


/8kT 

z ab ' Wc 


A. 57 


where a and y . are the collision cross section and reduced mass of the 
c ab 

two atoms. The collision frequency per atom a is 
z 


, _ 'ab _ „ „ / 8kT 

2 ab"^- Vc/iy^ 


Thus the HWHM due to collision is 


A. 58 


n, o 


Av c = 2TTT, 


be /8kT 


2tt / -rry 


A. 59 


ab 


Experimentally it was found that the central frequency of a line 
shifts with increasing foreign gas density. J To incorporate this into 
the theory one again assumes a classical oscillator with natural frequency 
o . As the collisions induce phase shifts, the oscillation frequency is 
a function of time. The oscillation displacement, x(t), is 


i / w’ (t')dt' 
x(t) = x e ° 

where x q is the oscillation amplitude. From the Larmor formula 


A. 60 


(51) 


for 


the power radiated by an accelerating charge and taking a Fourier component 
of the displacement one obtains 
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P(U3) 


2 4 i oo 
e or 2 
x 


4ire n _ 2 o 

o 3 lie 


i / to* (t')dt'-iwt 


/ dt e 


A. 61 


This expression has to be averaged over the phase shifts Aw' defined by 
io'(t) = to + Aw'(t) A. 62 

and over the number of collisions. The power radiated by an accelerating 


charge may then be written as 


(56) 


P(<o) = 


2w 4 (ex o ) 2 
3tt c 3 


nua 


2 - 2 

[a)- (w Q +nuS) ] +(nua £ ) 


A. 63 


where n is the number density of the perturbing gas and u is the relative 
velocity defined by 


/8kT 

n - j 

/ iry 


u = 


A. 64 


ab 


The Lorentzian in equation A. 63 is shifted from resonance by 

Av = 
s 2 tt 


A. 65 


where S is a constant having the dimension of an area which can be deter- 
mined experimentally. The Lorentzian has a HIVHM of 
nuo 


Av = -~ 

c 2tt 


A. 66 


which is the same expression as in equation A. 59. From the perfect gas law 
N P 


n “ V " kT 
one obtains 

Av_ = — P 


A. 66 


/ 8 


P/P 


c ' 2 " ° J ^Vab c JfJT 


A. 67 


For helium and sodium at the reference condition of 760 torr and 273°K 


-123- 


Av c = 5.57*10 


23 


P/P 


sr/r f 


A. 68 


for o in cm and Av in Hz. 
c c 


.(25) 


Experiments by McCartan and Farr' 1 * 1 for sodium in helium indicate a 
HIVHM per number density n of 

Aw 


(2.03±0.1)10- 9 ■ - r . a . d - /5 ? - % 
n atoms/ cm 


A. 69 


9 18 
at 415°K for the state up to a density of 11.5x10 atoms/cc. 

Equation A. 59 may be rewritten in terms of a , 


„ Aw , /lTy ab n ,,,, nn -16 2 

a c = ~ / W = ( 126±6 ) 10 cm 


A. 70 


At a temperature of 460° K they obtained 


Aw 


(2. 18*0. 2) 10 ' 9 


The collision cross section using equation A. 70 would then be (128±12)10 
2 

cm . The value for a from equation A. 70 may be substituted back into 
equation A. 68 


A. 71 
-16 


P/P 

Av = (7.04±0.34)10 9 — Hz 

C /F7T” 


A. 72 


For the shift at T = 460°K they obtained 

— -= (-0.07+0. 03) 10~ 19 r ^ d/s ^ c 
n v atoms/cc 


(25) 


A. 73 


By exchanging S with in equation A. 70 one obtains 
S = (-2. 1±0.9) 10~ 16 cm 2 


A. 74 


and the shift is 


P/P 

Av =-(l.l±0.5)10 8 — Hz. 

Sr7T~ 


A. 75 
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The cross section is a function of the temperature. This dependence 
is amplified at low temperature since the colliding particle is then brought 
under the influence of the long range part of the potential. The theoreti- 
cal investigation by Lewis and McNamara, based on the interaction po- 

f 581 

tential calculated by Bayliss v 1 for the ground state and resonance states 

of sodium perturbed by helium, indicates that the broadening rates per unit 

number density (HWHM/number density) scale as T ' . From the spectroscopic 

measurements performed in this experimental investigation (Section IV.D.l) 

0 38 

the rates per unit number density scale as T * . 


e. Instrumental Broadening 

The laser linewidth is influenced by the resonator in which lasing 
takes place. To find the resonator linewidth one must consider a colli- 
mated beam of radiation of intensity I launched inside a laser parallel to 
the resonator axis. After one round trip the intensity of the beam falls 
down to I-IS r where S r represents diffraction, absorption, scattering and 
output transmission losses. For a resonator of length L 


dl 



dt = ‘ 2L/c 


A. 76 


and 


2L/c6 


I = I e 
o 


A. 77 


2L 


The resonator decay time is given by — — and the expression obtained for the 

X* 

passive resonator H1VHM is^^ 


1 C ^r 
Av r = 2? 2L~ 


A.7S 


If it is assumed that the reflection at the optical surfaces and the output 
coupling represent the main loss mechanism in the resonator, then for the 
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Spectra Physics model 580 dye laser 6 r is approximately 0.04 which results 
in a HIVHM of 3.2 MHz. The laser lllnewidth is narrower then the passive 
resonator linewidth since the gain is nonlinearly related to the frequency 
dependent resonator losses. The resonator losses not only determine the 
conditions for the oscillation threshold, but also the theoretical spectral 
purity of the laser output. By considering the noise output produced by 
spontaneous emission into a resonator mode, which cannot be separated from 
the power derived by stimulated emission the following laser HIVHM, Av^, 
is obtained*^^ 



(Av r r 

2P 


A. 79 


where P is the laser output power. The laser HIVHM obtained by using 

_3 

equation A. 79 is 2.3x10 Hz for the dye laser tuned to the sodium 
wavelength with a power output of 30 mW. The narrowest linewidth achieved 
to date is 2kHz R.M.S. The reasons for the discrepancy between the 

experimental and theoretical linewidths are due to external influences. 

The main ones are the thermal and acoustical disturbances to the cavity 
length and the fluctuations in the thickness, position and index of refrac- 
tion of the dye jet. The single mode frequency v is given by 
me 



A. 80 


where m is an integer, ft is the cavity length and n is the index of refrac- 
tion of the medium assumed to be filling the complete cavity. Then 

1^1 - l^fl ■ 1^1 A.Sl 

To keep a 0.3 m long laser stable to within Av = 1 MHz the laser cavity 
length has to remain constant to within Aft * 6x10” 10 m for frequencies 
in the visible range of the spectrum. A model 580 single frequency Spectra 


....... a- . 
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Physics dye laser was used in this experimental investigation. A short 
term linewidth of Av c = ±7 MHz should be obtainable according to the 
manufacturer's specifications. To achieve this narrow linewidth an accum- 
ulator to damp mechanical, pump-induced, dye pressure fluctuations and 
a vibration isolated optical table was used. In the commercial laser, the 
dye was contained between two windows in the laser cavity. The laser was 
redesigned to incorporate a free jet dye stream to eliminate window burning. 
This also reduced the occurrence of mode hops but degraded thS frequency 
stability of the laser, 

The quality of the absorption spectrum depends on the ratio of the 
absorption width to the laser frequency jitter and on the rate at which the 
laser frequency is made to scan across the absorption profile. 
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4, Absorption, Saturation, Power Broadening and Radiation Trapping 
The complete absorption line shape including natural, Doppler and 
pressure broadening, is derived in this section. The derivation also in 
corporates saturation and power broadening. This is followed by a short 
discussion of radiation trapping. 

The solution to the steady state equation of radiation transfer may 
be written as^ 8 -* 


= l u (0)e <wL ♦ ^ (l-e Ka)L ) A. 82 

co 

provided and are independent of position x. yo) an d I^CL) are 
the impinging and emitted intensity from a column of gas of length L. 

and are the volume absorption and emission coefficients respectively, 
given by 


ftco 


e o. ' 4? A 21 "2 


A. 83 


and 


K «i ■ 37 bI 12 b i ( 17 - 77 ) «C“) 


A. 84 


A£^ and are the Einstein coefficients, n is the atomic number 
density, g is the degeneracy, g(co) is a linewidth factor and the subscripts 
1 and 2 refer to the lower and upper levels respectively. k^L is known as 
the optical thickness of the absorber. It determines the amount by which 


a collimated beam of intensity 1^(0) is attenuated as it passes through the 


medium. For k, ,L « 1 
w 

yto - vojo-ty.) * y 


A. 85 


The net absorbed intensity is 


I , = 1 (0)K L - e L 

w,abs oj ^ J to u) 


A. 86 


Neglecting the intensity injected into the beam by spontaneous emission 
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which is emitted in all directions 


Vabs ■ V°> V 


To evaluate the above expression, the population difference has to be 
calculated. The rate equation for the upper state of a two level system 
is given by 

— — — = S 2 (u) - n 2 (u)A 21 - [n 2 (u)B 21 - A - 88 

where S 2 (u) describes the rate at which atoms in the velocity class u, 
u+du are created. The number density and linewidth factor are velocity 
dependent since an electromagnetic wave of frequency w, propagating through 
an ensemble of atoms, interacts with a restricted class of atoms having a 
velocity component u along the propagation vector k given by equation A. 29 


where w 2 j is the transition frequency. The steady state number density of 
the upp<?r level is given by (equation A. 88) 


n 9 (u) 


o /•--■» . n pi I (w) g (u , w) 
S 2 (u) + n T (u)B 12 4 * 

A i f p 1 i p 1 ^ I(a3) 

A 21 + (B 21 +B 12 ) 4i 


where n^,(u) is the total number density of sodium atoms with velocity u. 
The steady state inversion density is given by 


n 2 (u) nj(u) 


n T B 2 

( g7 " g'x a 21 } p(u) 


The total number of atoms in a certain velocity group, n^,(u) , may be written 
as a function of the total number of atoms, n^,, as 


n T (u) = Hip P (u) 
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using the velocity distribution function P(u^. S 2 is defined in a similar 

n 2 n l 

manner. The total inversion density, - — , existing under zero field 

conditions (I = 0) , may be obtained by integrating over the total range of 
velocities. 

o o 

^2 _ : 

g 2 


g 


00 n_(u) n. (u) 

/ [ 


n„ 


1 


S' 


Si 


•]du = -( 


Si 


g l A 21 


•) 


A. 93 


If the thermal population of the upper level is neglected, S 7 vanishes and 


n° equals n^. Equation A. 91 may now be presented in different form. 
n 2 (u) njfu) 


o o 
n n 

c if ' £ p(u) 


*2 


1 + 


IQ) 


I s (u,w) 


I (u,w) is the saturation intensity which is the intensity at which the 

n 2 (u) n x (u) 

population difference, — , drops to one half its zero field 

value. Sodium radiation quenching is not included in equation A. 88 since 
the quenching of sodium by helium is negligible and sodium sodium collisions 
are unlikely for the densities encountered in the experiment. 

The explicit expression for the saturation intensity using the Einstein 
A and B coefficients is given by 
1 8TPflC 1 


I S (U,W) = 


g 2 A 

1 + — 21 
g l 


g(u,w) 


A. 94 


where g is velocity dependent and is the transition wavelength. The 
travelling electromagnetic wave interacts with the homogeneous part of 
the lineshape which is Lorentzian. The homogeneous lineshape is 
, Av, 


g(u,w) = 


T (V“21 )2 *( 4 V 2 


A. 95 
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where Av^ is the HIVHM (i.e. natural plus collision broadening) and 


= 0) - ku 


A. 96 


The preceding results may be substituted into equation A. 84 to obtain the 
absorption coefficient k (u) , 


n 


J, 

. , ftw D I g l 
K M (u ’ - 4i B 12 Si 


P(U) 


1 + 


g(u,w) 


A. 97 


8irhc 


1 


1 + — 
g l 


g 2 


g(u,u) 


21 


By defining 


V“21> 


1 


g 2 

1 + — 


§1 


Sit ficAv. 


21 


A. 98 


one obtains after some manipulations 


o o 

„ , , fito n I , n l n 2, 

4tt 12 '•gj g/ 


A. 99 


/, + IM 

/ i KJ 


7T 


Av. J 1 ♦ J&L- 
L/ X s U 21 


(w-a) 71 -ku) 2 +Av;[l+ - I -^- 


P(u) 


Ll W 


s^2V 

Therefore the absorption coefficient is proportional to a Lorentzian, with 
a HIVHM of 

I' TTTTo 

A. 100 


Av; = Av /l + v 

L L/ I (w 01 ) 


■s'- J 2r 

To find the contribution of all velocities, expression A. 99 has to be inte- 
grated over all u. Even though sodium sodium collisions are neglected, the 
sodium is assumed to be in translational equilibrium with the helium atoms. 
In this case P(u) is a Gaussian and since only the velocity component, u, 
in the direction of propagation k is of interest, the velocity distribution 


function is given by the one dimensional Maxwell Boltzman distribution. 


( 20 ) 
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. JL. u 2 

P(l0 = /JL- e 2kT 

1 J / 2iTkT 


A. 101 


A new variable is defined to simplify the notation 


A 1 WU 

A = ku = — 
c 


A. 102 


P(u)du = — — e 

/ir <5 


-(£) 2 


A. 103 


where 


u / 2kT B ^21 /2kT _ AV P 
c / m c / m 


A. 104 


AVp is the Doppler HWHM and is defined by equation A. 31. The last approxi- 


mation holds for 


21 -5 

— = 10 3 « 1 


as the maximum frequency scan of the dye laser is 5 GHz. Finally 


2 2 t n° „° 

k = ~~ — a_. C— - — ) 

0) 2 g, 21 g, go 


g l 21 g l g 2 / 1(0)) 

/ 1 + t — rrr 




where 


V(x) = / 


Av t\ + IM 
L/ V“2l3 u 

, ..2 . 2 ri I(w) , } 

(u.-a. 21 -A) + AV L [1* 


■ / toil l 
^ IT Av r 


-*.n2(^) 2 
D • 


The integral, V(x), in equation A. 107 is a convolution of a Lorentzian and 

. . 1 -.MnnUnMA m.m n M 4-1.A \frti ft+ MMAp. 1 A C^l) T 4* ■! A n A .. 11 . nl A AAfl 


a Gaussian lineshape known as the Voigt profile. 


It is normalized, 


since both the Lorentzian and the Gaussian are normalized. The following 

/* S *\ 

nondimensional quantities are introduced to simplify the notation^ J 
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a = 


x = 


y = 


Av G //to2 

v ' v 21 

Av G //£nT 

A 

Av G /An2 


A. 108 

A. 109 
A. 110 


After some manipulation one obtains 


V (x) = 


An 2 


Av; 

L 


(tt) 


3/2 


00 

/ 


2 2 


Av_ - 00 (x-y) +a 
b 


dy = 


(tt) 


3/2 


OO 

/ 


,-y 


Av,' , 

L - 00 (x-y) +a 


2 — j dy A. Ill 


where Av' is defined by equation A. 100, 

The shape of V(x) is determined by the parameter a only. The linewidths 
Av g and Av£ determine the scale of the frequency and intensity axis. The 
Voigt profile cannot be expressed in analytical form but has been evaluated 
numerically. 

The data analysis procedure employed in this experimental investigation 
requires the deconvolution of the natural and collisional lineshapes which 
are Lorentzian, The special case of a convolution of two Lorentzians provides 
a Lorentzian with a linewidth equal to the sum of the individual widths. 

The convolution of two Gaussians provides a Gaussian with a linewidth 
squared equal to the sum of the squares of the individual widths. 

Values of the Lorentzian and Gaussian profiles versus detuning off 
line denter in HWHM units, i Av, where i is an integer, are presented in 


Table A.l. 
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Table A. 1 


i 

0.1 

0.2 

0.4 

1 

2 

4 

10 

TrA'Vj^LCiAv^Av^) 

0.99 

0.96 

0.86 

0.5 

0.2 

0.06 

0.01 

/« 2 Av G S(iAv G .Av G ) 

0.99 

0.97 

0.89 

0.5 

0.065 

1.5xl0" 5 

8xl0" 31 


The Gaussian clearly falls off much faster off line center. For 
^ V L <<: ^ V G t ^' e ^ ne P r °fil e would be Gaussian to about three Doppler widths 
before changing over to the more slowly decreasing Lorentzian shape. For 
the opposite case, Av£ >> Av^, the profile would be approximately Lorentzian 
over the entire lineshape. When Av^ =* A\> G it is possible to determine 
both Av£ and Av^ from a detailed comparison of the experimental and theoret- 
ical line profiles.^®) 

Due to the presence of the electromagnetic field, the value of drops 
and its frequency response broadens. The broadening may be related to the 
shortening of the lifetime of the upper state by stimulated emission and 
of the lower state by absorption. The smaller value of is due to 
absorption followed by the reduction in the number density of the ground 
state. These effects, induced by the laser electromagnetic field, are 
known as saturation of the absorption line and power broadening. 

In the present experimental work the absorption line is obtained by 
scanning the laser frequency and monitoring the spontaneous emission inten- 
sity. A hyperfine structure is associated with each D line of sodium. 

Since the hyperfine lines are of different strengths, the radiation emitted 
by the focal volume might be reabsorbed differently by each transition. 

This would distort the spectrum. The effect is known as radiation trapping 
and is discussed more fully in Chapter III. 


Appendix B 

DATA REDUCTION BY LEAST SQUARES 


1. Introduction 

Errors are the difference between a calculated or observed value and 
the true value of some quantity. Since the true value is unknown one has 
to compare the results to those of other experiments or to a theoretical 
calculation, One may identify two types of errors. 

1) Systematic errors which are reproducible and can result from faulty 
calibration, background light scattering, etc. Statistical methods 
are useless in this case. The systematic errors have to be estimated 
by analyzing the experimental conditions. The systematic errors are 
then added to the statistical ones. 

2) Random errors are due to instrumental uncertainties or counting statis 
tics. These errors may be reduced by using more reliable instruments 
or counting more events. 

It is usually of interest to extract the most reasonable estimate of 
the theoretical parameters from the experimental observations, compute the 
uncertainties in the calculations based on the random errors and evaluate 
the confidence in the final result. The accuracy Of an experiment is a 
measure of how close the result of the experiment comes to the true value. 
Therefore it is a measure of the correctness of the result. The precision 
of an experiment is a measure of how exactly the result is determined, with 
out reference to what the result means. Precision is also a measure of how 
reproducible a result is, Clearly the uncertainties in the estimates of 
the theoretical parameters are related to precision while the confidence in 
the final result pertains to the accuracy of the experiment. 
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2. Distributions 

One defines a parent distribution which determines the probability of 
obtaining any particular observation in one measurement. The parent popu- 
lation is a collection of an infinite number of measurements from which the 
parent distribution is obtained. In practical experiments the finite num- 
ber of measurements is a sample from the parent population. Even if one 
cannot determine the true value, one would like to describe the parent 
distribution as well as possible. 

The experimental parent distribution is assumed to equal the theoreti- 
cal one in the limit of an infinite number of observations. Denoting by 
y^ the i observation of y and by N the number of observations, one de- 
fines the mean y by 


1 N 

y = lim rr 2 y . 
N-*° w 1=1 1 

the median y^y 2 by 


B.l 


probability^ <_ y^) = P r °bability(yj ■> y^ 2 ) = 


the most probable value y max by 


i = 1 , , . . ,N 


B.2 


B . 5 


probability (y ) > probability (y.) 

IJlciX X 

If the probability distribution is symmetric then y, y^ 9 and y max are 
all equal to the true value of the quantity of interest apart from systema- 
tic errors. For a nonsymmetric distribution, additional information per- 
taining to the causes for the deviations has to be supplied to determine 

which of u, u. or u if any represents the true value. By definition, 
1/2 max 

y , y. and y are just parameters specifying the probability distribu- 
i / /. infix 

tion. It is y which by convention is considered to be the best estimate of 


the true value. 
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2 

The variance a is defined by 


N 

a 2 = lim i 2(y.-y) 2 = lim Z y 2 ) - y 2 

N-*» " 1 N*» i=l 1 


B.4 


The sample mean includes only a finite number of observations and 
- 1 r 

y = n 1 y i 

IN i = 1 1 


B.5 


For large N 

y - y 


B.6 


The sum of the deviations from y Q , squared, is 

N , N 

^ = z Cy r yJ J 

i=l i=l 1 


B. 7 


This sum is minimized for y Q = y. The sample variance is defined by 

2 1 ^ 2 
s - £ 2 CXi-w) 2 
" i=l 


B.8 


However y is unknown. By substituting y for y one underestimates the sample 
variance. This is evident for N equals one since the sample variance van- 
ishes. Hence one defines a corrected sample variance. 


S 


2 

v 


1_ 

V 


N 

Z 

i=l 


(y^y) 2 


B.9 


v is the number of observations in excess of those needed to determine the 

(39) — 

parameters appearing in the equations. If y is substituted for y, at 

least one observation is needed to determine y. Thus v = N-l and 
N 

s l = ITT 1 (y i '^ 2 “ ° 2 B - 10 

i=l 

The last expression grows without bounds for a single measurement which is 
more in keeping with an idea of a variance. The most useful distributions 


are defined below. 
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a. Binomial Distribution 

If the probability for the success of a single event is p then the 

probability Pg(£,n,p) for the success of l events out of n is 

„ \ n! H,, .^n-#, 

P B (M » p) = &i(n-ii . )T P (1_P) 

and 


y * np 

o 2 = np(l-p) 

The maximum variance is obtained for p = 1/2. Then 


P B (A,n,l/2) 


n! 1 
JU (n-£) f 2 n 


y = n/2 

a 2 = n/4 = y/2 

and the probability distribution is symmetric. The probability distribu- 
tion is asymmetric for n i 1/2 and the variance is smaller then y/2. The 
ratio 


y 

a 


_np_ 


= 

v 1-p 


/np(l-p) 

improves as the size of the sample, n, is enlarged. 


B. 11 


b. Poisson Distribution 

This represents a limiting case of the binomial distribution where 
the average number, y, of successful events, l, is much smaller than the 
possible number, n, since the probability for a single successful event, 
p, is vanishingly small. Then 

p p O,y) = jr e" y b. 12 


B. 13 


-138- 


The advantage of this distribution is that it is specified by the mean 
only. It is useful in counting experiments under the above mentioned con- 
ditions. The ratio 

“ = /p B. 14 

improves as the square root of the mean. 


c. Gaussian Distribution 

This represents a limiting case of the binomial distribution for an 
infinite number of possible events with a finite probability for the sue 
cess of a single event. Then 


P G (y,y,co = 


i 


B. 15 


/2tt a 

The distribution is relatively simple and is accepted by convention and 
experimentation to be the most likely distribution for most experimental 
conditions where the Poisson distribution is inappropriate. The advantage 
of this distribution is that the most probable estimate of the mean, y, 
from a random sample of observations, y, is the sample mean, y. Generally 
one assumes either a Gaussian or a Poisson distribution for the measure- 
ment around the mean. Since these distributions are indistinguishable 
for most physical situations (y » 1) , one assumes that the Gaussian dis- 
tribution is obeyed. This is done due to the fact that the Gaussian dis- 
tribution is more amenable to mathematical manipulations, y is a parameter 
describing the parent distribution. By convention ir is considered to be 


the best estimate of the true value. 


3. Method of Maximum Likelihood 


It is often of interest to determine one characteristic of an experi- 
ment y as a function of some other quantity x. That is to find the function 
f such that 

y = f(x) B. 16 

Instead of making a number of measurements of the quantity y for one par- 
ticular value of x, one makes a series of N measurements y. . one for each 
value of the quantity x = Xj , where j is an index that runs from 1 to N 
and i in this case is one. The procedure is depicted in Figure B.l for a 


measurement y„ at a single point x^ . y^ is one value from the parent 
population associated with the parent distribution, P.. 

J 

If one assumes a Gaussian distribution, the probability for obtaining 
a data point y^ around the true value y^ with uncertainty ck is 


P. = P(y. . ,y . ,a.) 
3 3 3 3 


1 


/Trio. 

3 


4^> 2 

a J 


B. 17 


It is assumed that the mean of y^ , for all i in the parent population. 


is given by the fitting function f at x^ . That is 

"j ■ ft V 


B . 18 


CTj has to be determined from an analysis of the experimental conditions. 

The probability for making the observed set of measurements of the N 
values of y^.. is the product of these probabilities, which is known as the 
likelihood function 


N 

P = [ JI 


j = l /Tiro. 


] e 


i N 9 

2.\ i 3 - 3 r 

J" 1 3 


B. 19 


♦ 
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The method of maximum likelihood assumes that the observed measurements are 

more likely to have resulting from fitting function representing the 

parent distribution then from any other fitting function of the same form 

but with different coefficients. Hence, the fitting function representing 

the parent distribution is obtained by maximizing the probability or mini- 

2 

mizing the sum of squares, x » 


2 N y-.-fCx.) 2 

2 = e r— i-» 2—] 2 


2 [■ 


a. 

1 


B. 20 


The same procedure can be applied to a Poisson distribution. Since 
the resulting expressions are more difficult to manipulate, one generally 
assumes that the shape of the individual Poisson distributions governing 
the distribution of each observation y^ around its mean f(x^) are nearly 
Gaussian. The uncertainties are still assumed to be those associated with 


a Poisson distribution, o. = y... 

3 i J 
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4. Goodness of the Fit 

• 2 

The sum of squares is characterized by the variance of the fit, , 

~ 

which is an estimate of the variance of the data a . By definition 


\> " N-n 


I y JL j = 1 J 3 

N 2 3 

J-l a. 


The square bracket on the right hand side is the deviation of point i out 
of the parent population, y.., from the mean of t,:.- j th parent popula- 

^ J 

2 

tion, f(x^). 1/Oj is a weighting factor which reduces the N different dis- 
tribution into one to enable the sum over j to be taken. The sum in the 

denominator is a normalization factor. The expression for the reduced sum 
2 . 

of squares, X v > iS 

2 y 2 1 N 2 S V 2 

X 2 = 3 — = 77~ — l [-2J 3 —y = — B. 22 

A v N-n N-n . , 1 a. J — r 

' J =1 1 a 

where 


1 N 1 

N 2 
1=1 a. 

which for equal uncertainties ck = a reduces to 

a 2 = a 2 B. 24 

~~2 

a is a characteristic of the dispersion of the data y^ about the 

'2 

parent distribution j and is not descriptive of the fit. is a charac- 

teristic of the spread of the data about the fit but is also an estimate 
2 2 

of o . x v is the ratio of the estimated variance to the parent variance 

which makes it a convenient measure of the goodness of the fit. 

If the fitting function is a good estimate of the parent function, 

1 2 ~~2 2 

S y should agree with o and X v ~ !• If the fitting function is not a 
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good approximation, the deviations will be larger, yielding too large a 
2 2 2 

value for X V (X V > !)• X v < l» does not indicate an improvement of the 
fit. It is simply a consequence of the uncertainty in the determination 

of a. . 

3 


by 


The probability distribution for x with v degrees of freedom is given 
(64) 


±1 . xl 

P (x 2 ,v) ^ i 2 ^ 2 — — 

rcv/ 2 ) 


B . 25 


More important is the integral probability distribution for x defined by 

B . 26 


P v (X 2 ,v) = / P(x 2 ,v)dx 
X X 


This function describes the probability that a random set of N points from 

. 2 

the parent population would yield a value of x as large or larger. The 

2 2 
probability P^ of obtaining x v 51 1 is about 0.5. A larger value of x 

would yield a smaller value of P . 


v 


There are three sources of error contributing to the size of \ . 

t h 

1) The observations i from the j parent population, y^ , are a random 
sample from the parent population. Each y .. is measured at a certain 
Xj and has associated with it an expected value > given by the 
population distribution and assumed to be equatl to f(Xj). The fluc- 
tuations of the y.^-s about the expected value may be greater or less 

9 

than the uncertainties ck due to statistics. The value of x” may be 
improved by repeating the experiment. 

2) The parameters associated with the function have to be estimated accur 

2 

ately. This may be done using the least square method to minimize x • 

3) The choice of the functional behavior of the analytical fitting 
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function f(x) as an approximation to the true function <f(x)> will 

2 2 

influence the value of x • The resultant value of x for several 
different functions f(x) can be compared to determine the most proba- 
ble form for f(x). 
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5. Nonlinear Least Squares 

The least squares solution is straightforward if the fitting function 
is linear in the coefficient. For an arbitrary function one may extrapo- 
late to the method of nonlinear least squares. 

X is considered a continuous function of the n parameters 0^ (k = 

l,...,n) describing a hypersurface in n dimensional space. The space must 

2 

be searched for the appropriate minimum of X • More than one minimum for 
2 

X may exist within a reasonable range of the parameters 0^. It is then 

advantageous to conduct a coarse grid mapping of the parameter space to 

locate the main minima and identify the desired range of parameters over 

which to refine the search. Once this has been done, one may follow three 

2 

methods to obtain the minimum of x • 


a. Gradient Search 

All parameters 0^ are incremented simultaneously with their relative 

magnitudes adjusted so that the resultant direction of travel is along the 

2 

local gradient of x • 

2 


Vx = l fg- A0 
k=l 3(3 k k 


B. 27 


where A0j, characterizes the variation of x^ with 0^. Then define a dimen- 
sionless gradient 

^0 

30 k P k 




\ = - 


B. 29 


£ j ^ 

L '•30 J 
k-1 p k 

The procedure involves assuming an initial set of parameters and com- 


puting X and Then incrementing the parameters simultaneously by 
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«e k ■ -Y k Ae k 


B . 30 


2 2 
and recalculating x and Y k * These steps are continued until x is mini- 
mized. 


b. Linearization of the Sum of Squares 

2 

Instead of searching the x hypersurface for the minimum one may look 
for an analytical approximation to the hypersurface and use it to locate 
the minimum directly. 


Expanding x to first order around some starting point (B. ) . 

K O 


* x ° k.i *h 


[e k" (e kV * x o * E Ifr 

o k=l k 


iB k 


B. 31 


where 


56 k - % - 'Vo 


Equation has to be differentiated to locate the minimum. At the minimum 

2 

2 9y 

of x > the derivatives vanish 

dp(5 


iZ _ sxi 

33, " 8B £ 


+ 2 -ix 


o k=l 9 W*lc 


63 k = 0 


(Z = l,...,n) 


B. 32 


One may introduce the matrix elements. 


( 6) . - I &1 

2 3B X 


B. 33 


m k = 6B k 
(a) 


. .2 2 
1 3 X 


kZ 2 9B k 33 £ 


B. 34 


B. 55 


Then 


b = 


a 


B. 56 


The procedure involves assuming an initial set of parameters and calculating 
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a and 6. Then by solving the matrix equation one can obtain the new values 
of the parameters. These steps are continued until convergence is obtained. 
The gradient method suffers near the minimum as the evaluation of the 
derivatives by finite differences consists of the subtraction of almost 
equal numbers. The linearization scheme works best near the minimum when 
only first order terms may be kept. Marquardt ^ ^ shows that the path 
direction for the gradient and linearization schemes are nearly perpendi- 
cular. Thus the optimum search direction is somewhere in between. A pro- 

("35 33 cl 

gram provided by I.M.S.L. * ' has been used in this experimental in- 

vestigation. 


c. Linearization by Expanding the Fitting Function 

.2 


Instead of expanding x one could linearize the fitting function f(x). 
The function, f(x), expanded to second order is 


n 


n 3f (x) , 

f(x) = f (x) + Z -gf 63,. + 4 ^ 

0 k=l dp k 


n 3~f (x) 


k ' 2 A k=i 9(3 k 913 * 63k63£ 


B. 37 


If the above equation is substituted into equation B.20 and differentiated 
with respect to the parameter increments 6$ g (s is an integer) for mini- 
mization purposes one obtains 



. n 3 2 f (x.) 

1 V O v 3 so 

2 ... 33,33 c P k 

K=1 k s 


n 3 f„(x.) 
Z 


N 


1 


2 A 


0 -C- SB,] = -2 £ i* (x.)- 


n 3f fx.) 


j=l O 2 . 
J 3 


13 


Z 

k=l 


33,, 


1 1 B . 38 


. n n 

j Z Z 

1=1 k=l 


» W 


3f 0 Cx.) 


n 

+ ' Z 
k=l 


33-33 

3 s 


<53J 
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At the minimum l^r- vanishes. Then if terms of order higher than one in 
363 s 

the parameter increments are neglected one obtains 

N 3f (x.) N n 3f 0 (x.) 

3 3 

n 3 2 f (x.) 


8f .t»i> 
36. 


The above result may be put into the matrix form of equation B.36 after 
defining 


N . 3f (x.) i N t a 

( ^k = T ^ y ij -f o (x j^ 36. 3 = ' ~ Z ~ 

i=i o: J J k 

J 3 


O " n CsD [y^ 4 ( X -J ) 3 

2 j=l o 2 33 k 1J 0 3 
J 3 


1 A 

2 33,. 


B. 49 


N . 9f (x.) 3f 0 (x .) ^ 2 f (x . ) . 3V 

n _ y 1_ { ° 3 2—2 J— _ r v _f (v 11 } = i 0 

a k£ . 2 x 3(3, aa •v-S'S ^ J aft a« J ? 

j=i a. 


93 0 


ij 33 k 33 £ 2 33 k 33 £ 

B. 50 

These are the same definitions as in equations B.33 to B.35 
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' ** 


6. Uncertainties in the Parameters 
2 

The x criterion can be rigorously applied to evaluate the validity 
of the model only if the data are uncorrelated and represent a statistical 
distribution. Conversely, if one accepts that the model fully describes 
the physics, x way be used to measure the consistency of the data. Un- 
fortunately in most physical applications both of these conditions break 
down . 

2 

Choosing the data to be correct, x way be directly related to the 
model parameters. Define the error in the model parameter 3 k as 

A3 k = [(a _1 ) k]< ] 1/2 (k - 1 , . . . ,n) B.42 


where a is the matrix in equation VI.B.35. Changing 3 k by A3 k while 

2 

optimizing all the other parameters 3^ f k for minimum x > then the new 
value of x will be 1 greater than the old value. Denoting the likelihood 
function by P, then*' 40 ') 


3 k + 3 k + A3 k , X 2 + X 2 + 1 > P - “P B .43 

/e 

Hence A3 k corresponds to a standard deviation since by definition, the 
Gaussian probability distribution at y±a drops in value by /e ^ compared 
to its value at y. 

To use the method of least squares one has to assume that the uncer- 
tainties are in the dependent variable, y. 

a a 



B . 44 


where the subscripts and j ? correspond to two representative measurements. 
The above assumption does not always hold. The fitting will still be correct 
if the uncertainties in the dependent and independent variable are combined 
and assigned to the dependent variable. If the uncertainties in the 
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independent variable is larger than in the dependent variable it would be 
beneficial to switch their roles. 


Appendix C 
COMPUTER PROGRAMS 


1 . Data Acquisition 

The program reads NCH1 channels, NPPC1 samples per channel from the 
GMAD/1 A/D converter and NCH4 channels, NPPC4 samples per channel from the 
GMAD/4 A/D converter. The results are stored in binary form on tape. The 
records on tape are arranged in the following order. 

ITITLE (19) - contains the date and time, the run number, the sampling 
rate for the GMAD/1 and the total sampling time for the GMAD/1. 

HEAD (10) - contains the calibration of the linear position transducer and 
the laser frequency calibration. 

IPT (180) - Includes the position, stagnation temperature and wall tempera- 
ture in sequential order. The first half of the array contains data 
preceding the GMAD/1 scan and the second half contains data taken 
after the GMAD/1 scan. 

IVOLT (4000) - Includes the laser voltage ramp, the atomic beam device 
photomultiplier intensity, the tunnel photomultiplier intensity and 
the laser intensity in sequential order. 
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2. Data Reducat ion for a One Beam Experiment 

This step had to be divided into three programs due to memory limita- 
tions. Program DC12 converts from binary into decimal and rearranges the 
sequential arrays into separate ones for each device. It eliminates the 
points taken by the A/D converter before and after the laser scan and cal- 
ibrates the position and temperature. Program FT12 fits the experimental 
data with a theoretically derived fit with a least squares routine ZXSSQ 
supplied by IMSL. Program PL12 plots the theoretical curve, the experimen- 
tal tunnel and atomic beam data, the theoretical line strengths for the six 


D 2 transitions and residuals. 


-159- 


© 

CO 

cn 


O' 


to 


E 

Q. 

<M 

© 


jE 

U. 


© 

© 

© 

ill 

ts 

£ 


£ 

►-i 

CJ 

LfJI 

£ 

z 


> 

a 

<E 

Z 

m 


SB 

s 

00 

111 

z 

o 


o 

CL 

© 

© 

© 


a 

o 

lj- 

co 

£ 

© 

a 

LS 

O 

a 

CL 

CJ 


in 

w 

LD 

CL 

O 

-1 

to 


© 

CO 


© 

U3 


© a 

© Ll 

UJ 

w n/ 

I- 

D_ an 

CO 


© Li J 

S Ff 

3 P 

H *-* 


o 

-J 
CL 

Cl as 

Z CO 

© ©i 


(J ^ 


H O * 
% > 

~ Z 

© ' ^ © 
LD ^ h U 
/-* ^ a LD 
© cm in u ~ 

n n w 3 % 

an © LD a an 
Q □ J < UD 


© n * 

^ AS AN 

see 

d ‘ ‘ 

I— r * v as 

rr © 

»-« w w 

_j a a 

SiJd 

LU ^ ^ 


\ a* 

a © 

Li- O 
LU ‘ 
CL 


% at 


CL 
, CJ 

u 



cvj rr 
ro ro 

w w I 
LU 111 
-J -J 
H- H 

*— < I—* 

s- h : 



As 

m 


• 


LU 


CO 

AN 

AS 

m 

3 = 

cn 

CL • 

ro 

a a 

v 


o 

O 3 

UJ 

CL 

E 

~ HI 

i 

o 

S 

a 

CL O 

a 

UJ 


h- — • 

o ^ 

Z 1 

z * 

LU 

© • 

CL 

U! 

• UJ 

Ck CO 

Z h- 

uj ~ 

o z 

r 3 

*— • til 

H- a 

h- 

»-f LU 

© • 


u. f 

LU O 
CD 


© © 
« UJ 

CL h- 
f— -J 
to ~ 
•-» a 


1- 


a 

1 

a 




f-H 

CO 

a 

z 

© 





z 

H 

© 

♦— i 




N 

© 

z 

a 

o 




as as 

CL 

a 

CO 

z 




CJ — 



o 





H iL 

a 

• 

H 

a 




% * 

o 

f— 

© 

Z 


II 


AN A V 


o 

•— « 

© 



n- 

© 'C 

a CL 

a 

X 

a 


a 


tt in 

^ a 

a 




u 

o 

w w 

* © 


a 

a 


© 


w a 

to r 

o 

o 

o 


a 

H 

— i -J 

— ■ © 

K 

a 

a 



© 

H H 

a z 




a 

AS AS 

CD 

*■* ♦ 

© AS 

a 

a 

a 

a 

cn © 


H- H 

*-• ^ a 

cj 

cj 

CJ 

o 

r 


a* *■« 

a • a 




a 

a • 

P 

w w \ 

AN — < CM 1 — 

H 

r- »- 

z 

© a 

cn 

© 

a * — z 


*—» 

•— • 

o 

a a 

*— i 

\ 

•— > •. a a 

X 

X 

X 

CJ 

^ • 

o 


CO © 





a 

i 


in 


c*- 

CL 


LD 


£ - 










PAGE 0004 DC 12 11:02 PM MON., 3 MSP.. 1?S0 


i • ■ 


- 162 - 


7 ' 

t— * 


© 

ro 


o 


o 

CJ 





+ 

CO 
















2 

2“ 

© 











x 




*— < 

Cl 

© 











CE 




2 

uj 

© 




fVJ 




r*s 



2 




CO 

F 




£ 

ro 




*-4 



3 






% 



G 


© 



% 

> 





2 

S— ' 

(j 




a: 

o 

«— • 



F 

F 


• 



<t 

F 

2 

+ 




F 

1 



+ 

I-* 


F 



(_) 

-J 

M 

2 



2 


2 



© 

CO 


-J 


CO 

CO 

O 

71 

<Z 



♦—i 

O 



© 

* 

2 


2 * 

2 

> 

►—» 


/. 

CJ 

© 



CJ 

CJ 

/•“V 

CM 

| 

lu 


<r 

<r 

o: 

* 

if 

G 

CO 

It 


LU 


CO 



CM 

F 


cj c: 

u 

G t 


ni 

CJ 

*— « 

r\ 


CO 


*— 1 

* 

•— < 

w \ 

2 


CO M' 

to 

01 

^7 It 

F LU 

co 

^ It 

F 

UJ 

i— * 

•— 1 


o 

It 

F UJ > 

*— « 

* 

•— « «, 

*~-t 

c 

•> *—• 

«■ 3 

2 

* *- 1 

* 

111 

o 

1 

«*-l 

*—* 

F 

-J 3 M 



% *— * 




•— 2 

O 

»■* 


<d-m 

2 



1 


0 2^ 

01 


— w 

/-V 

U.l 

H CO 

w *-• UJ 

2 

n in 

>—• 

Lsj 


• 

II 


© 

> -* LU 

LU 

•-t 

II F 

•— « 

CU 

onjhl-D 


F — • 

F 

F 3 

F 

c 

►— t 

\*- * 


+ F •> 3 

CO 

w 

~ 3 

* : 

3 


JZ2 

UJ 


3 

2 2 

3 

LU 


F 


> 2 ~ 2 

a 

F 

O 

CD 

F 

in o 

O O ~ 

F 

ro o 

O 

O 

O 


— * 

-J 

o 

II o w ~ LU 

3 

3 

© > 

w 

u 

OJ « 

-> U F 

<x 


“> 

U F 



•— « 

O Q > (J F F 3 


CD 

ro w 


3 


2 

2 



2 

G 

CO 


> 


3 2 2 

U- 

> 


£ 

p 

CD 

o 

O 

U 

r 

DO 


O 

U 

P 

w 

o 

© 

> 



o 

ft 

X 

O Ll 

•— * 

CO 



M 




3 





2 

X 

£ 



LU 



3 




i— i 

UL 




O 

£ 

r 


3 

CH 



UJ 




tu 





(J 

2 

3 





co in 




in ro 






© — cm 





« 

• 


CM CM 

• 

• 




• 

• 





-« ro 

f. 

• 




o !Zi 

X ft 

a x 

52 


(J U U CJ 


u u cj a 


uuuu 


© 

ro 


u a o u 

rw 00 © © ~* cm ro in ld rs- © © © ~* cm ro ^ in co is. co cn © eg ro in cd rwco©©— ‘CMro<ffinc0N.CD 

©©©©©©©©©©©©©©©©©©©©©©SGD©©©©©©©©©©©©©©©©©© 


, l: x 


CD 

CO 

Cn 


- 163 - 


a 

CL 

oc 

a 


<r 


o: 


UJ 

2: 

2: 


Li) 

Dl 

? 

£ 

UJ 


Z 

O 


a 

u 


is 


I 

1 


ro 


• 




0 

0 

0 

O 






w 

CJ 




CL 

h- X 

0 

*— < 






\ 

<r 




LL 


tv 

h* 







Z 


X 



O -J 

* 

a 









0: 


> 

UJ O 

«■— » 

z 






•-* 


>- 

z 


H 

> 

II 

LD 






w 

ro 

r* 

-j 



\ 


a 

© 





X 


*-t 

\ 


cn 

© /■> 

% 

K 

ro 





1 


cn 



z 

ro 

/-“s 

CO 

•f 





cn 

% 

z 

'C 


UJ 

• * 

ro 


►— < 





0 

• 

uj 



p 

C 3 >-» 

*. 

Ci 

w 





a. 

z 

P 



z 

Z UJ w 

*«-* 

z 

CL 





> 

0 

z 

z y 


♦— * 

G • H 

w 

g 

+ 





w 

r 

M 

G h- 



CJ -1 

H 





/ s 


* 



CJ -J 


CL 

cn 2 p 

.J 

z 

© — ’ 

r\ 


© 


/-N 


CL 

cn 0 


LU 

~ rf > 


0 

ro ^ 

© 


ro 



z 

l'l 

►»5 £> 


CO 

\ % 

> 

►—« 

* 0 - 

ro 

>ws 

-L 


CM 

CL 

CO 

•. II 


a 

— 1 *-* n 


H- 


+ 


t—t 


w 


© 

•— t 


-j 

ll W 


M 

« CO ~ 

M 

W 

w 

© 

X 

Oi 

J 

u 

UJ 


— K fO 

L'l 

cn 

•-» 0 w 

V_' 

3 

■?, 

CD 

I 

© 


»— • «, 

© 

Hr 

— J * 

a cm 

0 

CL h- 

H 

K 

h- 

\ 


• • 

UJ 


z 

© 

— P ~ 

z * 

a. 

ro > + 

+ 

+ 

+ 

CO 


— • 

cn 

*w/ 

*— • 

O 

— • > w 

» k > « 

*—• CD 

© 

«-• n —« 

OJ 

ro 

LU 

0 © © 

© © ^ 


(V 

o 

Li- 


en 

UJ 

z 


<x 

LU 

Q- 


UJ 

a: 


<r 

CL 

UJ 


*-h ro H 


CM 


cj 


Di 

O 

Z 


OOO 
CL > CJ 


UJ 


o 

CL 


K 
-J z 

Ll o o 

*— > u 


CD 

<x 

cl 

£ 


0 D © S) 0 

« * * • * o 

(J 1 0 © Q Q « 

O If If ll IJ 

a ~ cm ro 

!> H H t— H 


cn 

o 

CL 

> 


f— I — ZD 
11 11 n 3 
cm ro — - 
I- H H K 
z 
o 
u 


ro ro ro ro 
\ \ \ \ 
— • r\j ro ^ 
h* ' 

It 


§1 

fs 


I! II 

rvj ro 

I- K 


II 

H cn 
n o 

^ 0 - 
H X 


to 

© 

Cl 

a 


in 

© 

© 

© 


ro 


ro 


CJ CJ CJ CJ 


CJ CJ CJ 0 


CJ (J CJ O u o 


CJ CJ CJ 


LU 

UJ 

a 

a. 


tn©- 4 Cvjro^intof s w©j>©»-*ojfo 5 y 2 , P!ir a:j ^CD— *fMro^in i J?rN-rocri©^r\jro , Q‘inuDo-©cri© 
MDrwr'-r^r^r^rvrvr-rvr^oocDaocoQDCDCocooDGDcricocncoo'ieno'iCoc'iO')©©©©©©©©©©-* 

©©©CD©©©©©©©©©©©©©©©©©©©©©©©©©©©©©©©©©©©©®® 


i 


PAGE 0006 DC 12 11:02 PM MON., 3 MAP. 










PAGE 0009 DC 12 11:02 PM MON., 3 MAR., 19S0 


-167- 


CL 

s 


3 

CL 


r j 

CD 

g 

u. 


cm 

CD 


3 

O' 

X 


s 


r 


E 

3 

Pi 






(_) 






U 

E 









3 






3 

X 









-1 






-J 

CL 















% 

© 







r“. 


CD 


s-\ 




CD 

Cl 

/*“ s 






L_' 


« 






* 

CL 

y/ 






CO 


H 


CO 




h- 

CL 

co 






►-I 


D- 


*— t 




CU 


»—i 






a 


LJ 

E 

Pi 




CJ 

CD 

p 





z 

►— t 


ID 

X 

►—1 




3 

Z 






Cl 

* 


-J 

CL 

% 




-J 

»- 1 

•> 





cu 



* 

CD 

\/ 




* 

h- 

s - 





CD 

CL 


X 

O 

rv 




X 

H 

CL 





O 

H 


X 

CL 

F 




5c 

O 

F 




^7 

or 

cn 


X 

Q. 

cn 




X 

P 

CO 




•—* 

CL 

»— • 


rvj 


*-i 




CM 

CL 

H 1 






* 


« 

CD 

% 




% 


* 

y~ 



O 

CD 




X 






z 


5 



H 

Z 



CM 

l-H 





CM 

o 






•— » 

CD 


CD 

t 

CD 


E 


© 


N- 

F 

T“ 

CO 

o 


•-i 


g 



cn 


CD 

F 


3 

CO 

l: 

CD 

H 


o 

-J 


x 


E 

3 

_J 

© 

X 

u 


►— < 

Ll 


Ll 

-J 

D- 


CL 


Ll 

© 

CL 


X 

X 


LL 

X 


X 

CL 

X 


X 


X 

*— < 

X 

CL 

X 

p 

** 


CD 

CJ 

CD ^ 


CD 

E 

CD 

CJ 

CD <~s 

£ 

CD 

H 

cu 

1 

Z 

* 

p 

* *3* 

E 

* 

x 

* 

p 

* ^ 

y 

cn 

«. 



<x 

cn 

ro — 

X 

cn 

C_' 

cn 

ro ~ 

to 

CTi 

*— i 

cn 

e 

• 

DC 


a: 

CM 1 

CL 

w 

CD 

w 

X 

CM 1 

*— 1 

w- 

P 

w* 

CO 

a 

CD 


CO 

W 

CD 


e 


CO 

w 

p 





UJ 

P 

CD 

u 

CL 

o 

P 

LJ 

a 

CL 

LJ 


(J 

X 

u 

p 

• 

fe J 

p 


h 

CO 

w 

✓ N. 

to 

M 

s 

P 

S! 

H 


f- 

X 

Ll LU 
X 3 

LU 

in 

*H 

LU 


in 

•— « 

LU 

X 

UJ 


LU 

x 

LU 

p 

LU Z 

CO 

z 

X 


LU 


X 




UJ 


p 


»— i 

X 

w 


-J 

CD 

p 


-J 

K 

-J 

CD 

-J 


3 

f— 

3 P 

LU 


z: 

1 

CL 

P 

z 

P 

O 

-J 

CL 

-J 

z 

P 

CD 

-J X 

P 

ll 

3 

<r 

X 

X 

X 

X 

-J 

<x 

3 

X 

3 

X 

3 

x a 

LU 

*— i 

CL 

CJ 

CL 

U 

CL 

u 

CL 

a 

CL 

LJ 

CL 

CJ 

X 

u u 

CL 















\r 



• 


• 


• 


• 


• 


• 


• 


• 


• 


• 


• 


• 


• 


• 


• 


• 


• 


• 


• 


• 


* 


• • ■ 


# 


• 

o a 

u a 


u u u u 

u 

c 

c 

c 


cj a u u 

u 

c 

c 

c 

U U 

u u 

u 

p u a 

u 

u u 


rs.oocr»©-*CMro^incor^oDcn©-«(Mro^incor^CD 

rororoM'<tf<3‘^M‘*TM ,, T^T , 'J , tnininLnLninininin 


PJ \l V V \ 'J V M m mujujuju;ui 

loronrororororororororororororororo 

CD CD CD CD CD CD CD S3 CD CD CD CD CD CD CD CD CD CD 



3 



to 
( T ) 

O 

CD 

fl 







0383 SUBROUTINE COL (CF1,CF2,XN,FR) 

0384 FR=CF 1 CXN**2 . ) +CF2s?<N 

0335 RETURN 

0386 END 



PAGE 0012 FTN. 11:02 PM MON., 3 MAR. , 1930 



- 171 - 


GD 

CO 

C 7 ) 


ro 


Z 

Q. 

LD 

ro 
• • 
o 


* 

Li. 


© 

CD 

© 


L 3 

<r 

CL 


CM LL 

21 2 

a: 

cd or co 
* ld « 
OQC 
CO CL LU 
CO CL CL 


<X 

z 

LU 


CD 
Z 
«— * 

1 = 

o 


CL 

z £ 

G 

UI Z 

CD tO 

t —4 

LU A 
Z 

o z 
o 

CL 

o <x 

IL F 

g 
LU Ci 

z 

*-« _j 

_J <x 
u 

a F 

U LU 
F CL 
F CD 
UJ LU 

LU 

P 5 g 

cr 
a 

o 

■ CO 
. CL 
CO LU LU 
F Z m 
z — £ 

*—* CL HD 
QUIZ 

°-§ch 

CD LU Z 
CD UJ O 
fs- tO X — 

in £ £ - 

tsg 

LL 3 tO 



in 

















in 






ft 











w 






/-ft 











CD 






CL 









/s 


ft ft 






LL 









3 








LU 









3 


uo cd 






CL 









CL 


* 3 






> 









a 


in 






ft 











w 






/-* 








ft 



F w 






CM 









/—ft 


Z Z 






fO 








CD 

ro 


LU CL 





ft 

w 








t— 

w 


)» O 


•, 



/•N 

UJ 









to 


►- 3 


r-\ 



LD 

_J 


% 




ft 


/-> 

r 


»\ ft 


F 


ft 

Z 

F 


/N 




/—> 



a 


r\ 




/■ft 




$—* 




GD 



a. 


cd in 


LL 




F 


z 





/—» 

X 

\_/ 


ft 


•— * 


z 

X 





ft 

■ 

CD 

DC 

ft—' 

a 


f\J w 




2 - 

ft 

■w 


LU 


,—ft 

/—> 

ft 



/—. 


w •— « 


/“v 



/-ft 

ft 


* 


CD 

Q 

/-ft. 

Z 

/-* 

iL 


CO p 


N 


X 

UO 





F 

F 

C\» 

CL 

— • 

LO 


LU F 


»—• 


ft 

CM 

Cl 


CD 


ft 


LO 

ft 

Ll 



z l “i 


w 


/-* 

ft-’ 

ft 


ro 


/—>. 

/—ft 

ft— * 

-ft 

ft 

A 


»-< x 


LU 


CM 

LU 

/-ft 


>-♦ 


^r 

CO 

LU 

in 

.—s 

*— i 


-J ft ft 


cJ 


(M 

pJ 

O 


Ul 



*v 

_ I 

uo 

ro 

ft 


ft sr\ — < 


F 


w 

F 

ro 


_i 


V.' 

w 

F 

w 

s—/ 

/—ft 


Z- in r 


*-+ 


LU 


w 


F 


L'! 

LU 

Z 

LU 

X 

ru 


X - ft 


F 


-J 

F 

UJ 




-J 

-J 

F 

-J 

w 

ft/ 


3 w ^ 


♦— i 


F 

>-H 

-j 


F 


F 

F 

»—» 

F 

ft 

to 


X P CD 


w 


»■— < 

■w* 

F 


*— S 


►—i 

l-M 

w 

*— * 

/—ft 

T~ 


* F CM 


ft 


F 

ft 



ft-' 


F 

F 

ft 

i— 

G 

G 

* 

^ 5 '-; 


/-* 

ft 

*— < 

ft 

F 

ft 

ft 


*— » 

•— i 

/-V 

»— » 

ft 



— X CL 


Z 

/->» 

V-' 

-J 


/■*. 

/—ft 


*w 

*w/ 

_J 

>— 

ft 

w 

CD 

ft ft X 


c 

_J 

ft 

£ 

>— ■ 

in 

T~ 


ft 

ft 

LL 

ft 

OJ 

ft 

CD 

GD ^ UJ 


CJ 

LL 

/-V 

X 


CD 

-J z 

.—ft 

/—ft 

UJ 


w 

/-ft. 


*— 1 (S 3 <«, 


LO 

LU 

— * 

3 


a. 

> X 

G 

c_« 

CL 

Ln 

>c 

\/ 

w 

% O ~ 

UD 

*— •« 

CL 

G 

~ 

F 

X 

LU 

LU 

F 

F 

> 

\ ■' 

w 

CL 

> 

~ r^ in 

in 

ft 

ft 




t, 

* 

■. 


ft 

ft 


ft 

F 

F 

in w ft 



.'-v 


/-ft 





/— > 


/*. 

— V 

/—ft 

CO 


OFF 


z: 

o 

CL 


CD 

UD 


CD 

in 


© 


in jo 
in cd 


0 O N O ro OJ N 

(M N ro ID in UJ N 

1 I I I I « I 

*-« «-4 <~l If) 

— cm ro ^ in to n. 



F F F ^ ^ 

^ — in 

H H H ^ ^ 

~ ~ x X 

S— W W *w v 




LU 

© 


.?"« tM ro ^ 


— cm ro ^ in ud rv 



uuuuuuuuuuuu 


GDOOOOO©DO^^^-<^^^*-^^(%jfvjr\!(\iCVJOJ(MCMDJ 
OC(SS)(S)(DQQCDOOOOO©( 2 CCOCDS)©D<SQS “ — 
OCDQQOO0OS)S)0OOCD0ISGCD0G)S)QQS)O 


Ln 

UD 

N. 

CD 

CO 

CD 


tM 

ro 


in 

UD 

rv 

CD 

cn 

G> 


CM 

CM 

CM 

CM 

OJ 

CM 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

sr 



GD 

CD 

CD 

CD 

GD 

GD 

CD 

CD 

CD 

GD 

GD 

GD 

GD 

CD 

GD 

GD 

G 

CD 

GD 

GD 

CSD 

GD 

GD 

GD 

GD 

G 

CD 

CD 

CD 

GD 

CD 

CD 

CD 

GD 

GD 

CD 


} 


■ 172 - 






/-N 















4 

© 



© 

£ 









<-s 












\r 

v—S 
















0 D 

4 


4 







0 D 






© 

■ 






• 

co 














h- 


• 

CD 


CO 

















00 

• 














4 


CM 

• 


• 







CJ 






CM 

4 





U. 

CM 
















^4 

CM 


CM 













— < 

© 





* 





x 










• 


Ll 










Lu 





CO 

X 

X 





in 

LL 




© 










co 


4 

ll 


Ll 







4 





• 

4 

p 






4 




to 










Ll 


in 

4 









in 





r.J 

LO 

4 





P m 




X 










4 


CD 

tn 


in 







CD 






X 

CD 





*, 

’XI 




At 










in 


r- ^ 


y 


4 





H 





LL 

P 





s 

p 
















4 

H 


P 


© 


4 




4 


4 


* 

4 

CM 





cn 

4 

/*» 



Cl 










h- 


£ 



4 


CO 


© 

4 


£ 

G 


© 


LO 

£ 






h 

£ 

CO 



Ll) 










4 


1 

£ 


S 


p 



© 


1 

CD 


X 


1 n 

1 

Ll 






1 

• 

£ 


LO 










£ 


+ 





4 


P 

CD 


-J- 

H 


H 


P 

4. 

4 





CD 

+ 

CM P 


O 


CD 





4^. 





C 

It 


II 


CO 


4 

H 


E 

4 




4 

=’ 

Ln 






£ 

•-4 

IZ 


© 







UD 



II 


4 





• 


© 

4 


4 

<sT 




■ £ 

4 

'4 





C> 

». 

LU 



CO 


4 





4—4 


n 



© 



or 


oi 



© 


o 





1 

© 

l— 






© 

* 

© 


~7~ 


in 





* 


rf 

© 


l£i 

CD 


o 


— 'i 


CM 


UD 

CM 


nj 



CD 

4 





P 

CO 

*0 

0 . 


© 


At* 





m 


• 

♦ — < 


1 “ 

1 - 


H 


LL 



CM 


P 





C 

P 

S 





Z 





CL 


P 





V 


OJ 

CO 


4 

l 


CJ 


4 


IL 



4 

X 


Ll 


4 

4 






co c: 

4 

P 

z 


P 


4 





H* 


4-4 

cn 



O 


rr 


IjO 



Ll 


ro 

4 


4 


© 

© 

II 





^ LJ 

r^» 

% 

u 


1 


• £ 





4 

r*. 

LU 

© 


• 

*> 


LL 


v 7 


in 

4 



— M 


— 


X 






* »-« 


‘£ 

© 


X 


© z 

© 




£ 

rr 

4 

© 


C>J 





H 


V 

Ln 


CJ 

rr 




K 

w 




Cl 


in ll 

CO 


XI 


© 

■£ 

LL O 

►— * 




CO 

• 

uo 

© 





V 


4 


P 




K 


H 


4 


CM 



CO 


\T ^ 

ll 

1 ! 



to 

tu 

*. *— 1 

4 

£ 



X 

fM 




LL 

X 


h- 


r 


4 

P 


X 

4 


4 


© 

X 

At 



CTi 


P z 

4 


P 


s * 

V 

ru 

Ln or 

in 

LU 



LU 

—4 




4 

H 


4—1 




£ 

4 


4 

= 


£ 


4 

it 





i CD 

in 

• 

c 


S—' 

© 

LU 


CL 


r-\ 

> 

LU 




— • 

»— • 


cn 


II 



r 


in 





CM 

in 






& *— * 


G 




Z P P P 

© 


CO 

~ 

4 

= 




CO 






H 



"v 

i; 


It 



T? 

B 



4 


CO 

t — 

CO 

00 


or 


4 4—1 


*7- 


4—1 

*— 4 

in 


o 


(- 

T~ 


LU 





II 


H 





X 

p 

X 



• 


CO 

•> 


© 


P 

LU 

= © 

'£ 

cn 


4 



II 



4 

LU 


t 


nj 





4 

■> 


/—v 


4 

4 

H 



x 


Z LL 

£ 

X 

cn 


r* 

CD 

x 

or 

UJ 


in 

or 

K- 

o: 



S 

F— 




X 


rv 

r- , 


= 

U 


X 


in 

s 

© 



© 


o o 


o 

-j 


X 

© 

ll 

LU 




o 

4 

UJ 

© 



z 


►-4 

CO 

* 


-u 

M 



X 


X 


M* 


X 



z 


1—4 

II 




LU 

or 

CO X 

H 



H 

ll 

= 


4—1 


II 

»— » 




w 


z 

X 


II 

cn 


cn 


K 

It 



/— N 



l- X 


Z 

¥—4 


> 

o 

z u 

L'J 

-J 

= 

4 



LU 






LU 

r.j 



w 




\ 


js 


4 


© 


E 



a l*j 

N 

© 

CO 

*— > 


p 

o z 



X 

£ 

or 

(1 

z 

h- 

r-'. 


t 




Z 

>-4 


^ 1 


/“\ 



w 


S 


w 


V- 

ro 


or dd 


CO 

LU 

£ 

cm 

CO 

*— LU 

© 

—I 

4— • 

LL 

LU 

• 

© 


© 

nj 

X 

0 D 

4—1 

LL 

X 

© 

z © 


© 

X 

'w' 

© 


© 


X 


© 

*— • 



X Z 

Z 


X 

T 

\ 


1 — CD 

cr 

© 

or 

O 

H* 

& 

o: 

LU 

• 

X 

■N - 


© 

4 


J, • 



s_ ‘ 




1 ! 

X 

T* 


X 



P © 

w 

© 


© 

-4 

u 

© a; 

© 

u 



UJ 

© 

© 

cr 

00 

*~ 

LU 

CM 


tn 

X 

r-j 

3 nj 


r: 


> 

nj 


r^j 


o 

CD 

CM 

p 

4 


~ Z 


© 

LL 

»— 4 



© LU 

© 

* — i 

© 

CO 

Z 

c_ 

CL 

o 

—4 

N— ■■ 



• 

CD 



~r~ — < 

“ 

— * 

X 

1 — 


P 


X 

K 



© 

• 



z © 

o 

CD 


or 

© > 


h- Z 

f— 

© 

as 


-j 

Ll 


£ 

Ll 

X 

f— 

z 

Ll 

LL 

X 


cr 

*— r 

X 

< — < 

X 

X 


P 

X 

z 

Z 


Ll Q 

re 



o 

—4 

l — 

© z 

or 

LU 



or 


P 


4 

X 

z 

4 

G 

4 

© 

4 

Z 4 



ZD 

X 

4 

X 

« 

DD 


X 

4 


g 


O X 

-i 

*— 1 

or 

u 

+ 

LU 

> o 

o 

or 

Ct 

CD 

© 

ll 

n 

LL 

in 

u 

© 

in 

Ll * 

£ 


LO 

o tn 


LO 

L_ 

o 

in 

o 

LO 


X 

X 

m 

V 


r-> f— 

Z 

CO 

p 

© 

z 

i— 

LU CJ 

or 

o 

|i J 

4-4 

CL 

o 

or 

© 

UD 

7~ 

Cl 

UD 

X 

1 

M 

UD 

— UD 

© 

© 

© 

X 

X 

X 

P 

15 

rv 

co 

X P 



U) Dr a: 


LU P 

*-> 

W 

Z, 

1 

or 

LU 




© 


p 

UJ 

G 

K H* 

+ 

P 

1- 

co H 


H 

X 

X 

f- 

X 

i— 

T~ 

G 


P 




* us Z 

P 

or 

CJ 


1 

>‘ 

• or 

LU 

X 

or 


or 

• 

© 

o 

4 

© 

G 

4 


£ 

z 

4 


CO 

4 

X 

e= 


4 

> 

4 


© 

cn 

4 

p 

z 


2EP 

CD 

i 

LU 

x 

P 

CO 

(J LU 

i 

r— UJ 

• 

o 

y 

a 

•— 4 

£ 

G 

v -- 

= 

• 

4 

LU 

£ 

i r 

co 

= 


£ 


E 

z 

CO 

cn 

£ 

z 

X 

• 

*— < 

a 


© P 


Z Ll 

or 

J 

-> 

as 

or 

© 

or h- 

1 

LU 

CJ 

f 

X 

o 

cr 

1 

© 1 

© 

1 

z 

0. 

1 

X 

1 

L-J 

X 

O 

i 

X 


cm x g cn 

a 

co 

| 


D 2 g 5 

✓s X ' 

LU > 

*7* 

*— < 

or 

o 


© 

+ 

X 

© 

a. 

© 

CD 

o 

+ 

o + 

© 

+ 

X 

X 

+ 

X 

4- 

© 

X 

X 

+ 

X 

X 

CM 

Z LJ 

*> 

co 

Ll 

X 

4—4 

r— 

4—1 

co 

LU 

zr 


X 

= 

u_ 

CD 

= 


K 



CJ = 

as 


I — X 1 


X 

S 

Z 

X 

X 

£ 

►— ■* 

ro 

ll 

a 





































• •! 

in 

— 4 r = 

sr 

= 

£ 

s 

s 


s c 

- 

£ 

s 

= 

c 

s 

s 

c 


s 

r 


£ 


c 


c 

c 


c 

c 


£ 


s 

£ 

£ 


c 

CD 

w 

WWW 

w 

w 

w 

w 

w 


w w 

W 

w 

w 

w 

w 

w 

w 

w 


w 

w 


W 


w 


w 

w 


w 

w 


w 


w 


W 


w 


p h h 1— 

P 

P 

P 

p 

P 


H P 

1 

P 

f- 

P 

p p p p 


P 

K 


K 


K 



H 


P 

H 


H 


H 

H 

P 


p 


§ 

© ©a 

© 

© 

© 

© 

© 



© 

© 

© 


© © © 
z z z 

£ 


© 

© 


© 


© 


© 

© 


© 

© 


© 


© 

a 

© 


© 

CVJ 

z z z 

Z 

zz 

z z z 


r r 

zz 

2Z 

r 

z 



z 

Z 


Z 


r 


z 

z 


z 



z 


Z 

Z 

4— 


z 

•»* 

nr 

or or or 

or 

or 

or 

or 

or 


or or 

or 

or 

or 

or 

or 

or 

or 

or 


or 

or 


X 


X 


X 

X 


X 

X 


cr 


CL 

X 

cl 


X 

P 

a 

o o o 

o 

o 

o 

o 

o 


o o 

o 

o 

o 

o 

a 

o 

o 

o 


o 

o 


o 


CD 


o 

o 


CD 

o 


o 


CD 

o 

o 


o 

Ul 

u_ 

Ll ll X 

ll 

Ll 

Ll 

ll 

LL 


Ll Ll 

LL 

LL 

Ll 

Ll 

Ll 

u. 

Ll 

Ll 


Ll 

ll 


Ll 


LL 


LL 

X 


X 

X 


X 


X 

X 

X 


X 









^4 















•—4 








•— 1 


•— i 







ro M* in x 

IV 

ro 

© 

© 

—4 


cm ro ^ in 

UD 

iv 

© 

cn © 



CM 

ro 




in 


CD 

V 


CO 

© 


© 


«— * 

CM 

ro 


tn 

rvj 

© 

© 

pn 






















M* 


*sr 





M* 

\r 


in 


in 

in 

LO 


in 

Ul 

X 

(J3 

ro ^ in x 

X 

© 

cn 

© 


cm ro 

in 

UD 

n- 

CO 

cn 

© 


CM 

ro 

in 

CD 

IV 

© 

cn 

© 

— © 

ro 

M* 

in 

CD 

IV 

© 

© 

© 


OJ 

ro 


a 

^ ^ 


*sT 

M* 

in 

in 

in 

in in do 

in 

in 

in 

in 

CD 

KD 

CD 

CO 

'JD 

CD 

*JD 

•JD 

( JD 

CD 

IV 

rv 

IV 

V 

IV 

fs- fv- 

rv 

V 

© 

© 

© 

0D 

© 

CL 

© 

© © © 

© 

© 

© 

© 

© 

© 

© © 

© 

© 

© 

© 

© 

© 

© 

© 

© © 

© 

© 

© 

© 

© 

© 

© © 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 


© © © © 

© © 

© 

© 

© 

© © © 

© © © 

© © 

© © 

© © © 

© 

© 

© © 

© 

© 

© © 

© 

© 

© 

© 

© 

© 

© 

© 

o 

© 

o 

© 


\ % c> 

O \ I- 

* 04 s 

Ou 

f~ cr 
in d. ^ 
c_ o v: 

i*i * 


•j * », 

CO * • 

-* u n ld 

• 

n. 

nj 

G-k \ \ 


* ro 

fv CT) 
IT' O 

co ro 

* © »Xt 

• • in 

— • co r- 

r \ ? X 


co 

* a * 

\ y: \ \ 

*o \ — • r.i 
\ _ o o 

cx £ d \ 
• — • - c * 

Cn G * £J 

O ~ 

o 

GO 

GD 


« n 
uj g 
cj 3 

Z3 X 
Cj = 

LU * /-» 
CL 01 = 


h- H H 

x s: x 

O' 0' 01 

o o o 

LtU.ll 
CD fv CD 

LH in in 


co r- co 
c j ^ w 
CD CL 01 
\ 01 01 
LU LU 
O 

u 

\ \ 
to z ^ 

eir =? 5? 


l-hh3 
<x cr cr 
2 : x £ <r 

CL or LU b- 

o o o <r 

LL LL Li- Cj 

<D «-* CVJ 
U) CD UD 


DC w CTi 
CL DC £V — 
CL UJ CL 


LU * \ 
» \ 

LU LO cr 
* .J \ H- 
«. DJ CD -J 
in h z ^ uj 

U CL C Z Q 


ff ? ■ g £ £ £ £ 

<rcrcr<x<r<rcr<x 

ccjcqocjcjcj 





GD 






* 


~ r j 






DJ * 




CL 


* ll 

X 



h“ 


X CL 

o 



cr. 


CL H 

CO 



►— • 


1— co 






CO •— 

p 






l- 



in 


*. d 




•k 


© a 

01 



LU 


O 'C 

LU 



-J 


v — 

f 



H 


*— * « 





. u- 

CT 



H 


✓ v * 

U- 


V/ 

*— * 


LU UJ 



CO 

* 


* ^ 

X 


►—1 


CO 

S' 

Cl 


« 

CO 

*—* 

CO co 

cl 



•-H 

o 

►— < 

U. 

LO 

r 

A 


cj a 


X 

o 

H 

£ 


<s 

cr 

Cl 

% 

o 

.* *, 

H 

CL 

LL 

•— « 

CL 


cr 

<-> 


'w' 

LL 


CJ 

01 

UJ 




<x 

Ff 

CJ 

<X 

o o 

C L 

CL 

W 

K 

LU UJ 

UJ 

n 

•*-« 

X 

cr 

X DC 

tu 

CL 

H- 

UJ 

C» 

LU LU 

CO 





X 

-J 

» 

-J 

CJ 

U -J 

cr 

-J 

cr 

-J 

cr 

-J _J 

£ 

cr 

LU 

cr 

LU 

d <Z 

o 

01 

u 

01 

U CJ 


o u o o uuuu u u u u u 


Ln«X)ivoooocD — curo^Lnuo^cDcriCD — iMro^rmuDiN-oocncD— 'C\iro^LnuDN.rocna^ 

CO CO GO 03 CD O'! 0"i CT) CD CO Ut CTc CO CTi CO © © © Q GO © © <55 © © — 1 —* — • — -* — * — » — CM OJ 

C50 Q O O Q O CD CD Q Q O (S) Q © G ^ ^ ^ >M 

©OOCQO©SOQS)QGlGQ0GlCC2O©C©D©C©©O©OG>(rim/nm/n 


5122 C 

5123 C 

1124 C 

1125 C.. MAX. ON LEFT PEAK OF EXP. INTEN' 

1126 C 







- 175 - 


© 

co 

0*1 


DC 


ro 


UJ 

rz 


3 

o 

DC 

CO 

ZD 

co 


Z 

O' 

a: 

CL 


CL 

o 


z 

ft 

E £ 

•• LLI 

g —* 
h * 
-J DC 

LU LU 
Cl tc 

A CC 

CO 

0 _ * 
LU 1C 
a DC 

(S CD 

*— 3 

CO ■> 


Ll 


u 

g 


Cl 

-J 

cu 

-I 

cc 

ro 

ri 

>— * 

co 

L‘J 

CJ 

K 

O 

li- 


ce 

3 


o 

-J 

CL 

_l 

G 

3 

Cl 

CO 

UJ 

CL 

or 

o 


Z 

CL 

li 

x 

E 

CL 


U.i 

P\ 


UJ 

CD 

g 

CL 

F 

CO 

o 


i 

CO 

3 

© 


o 

CO 


A 

*-4 



o 



LL 




ci 


Z 

LL 



X 


U- 






w 










£ 

L'J 



* 



*— 4 


o: 





*— < 


C 




cr. 








*“> 

3 










CL 



LL 



CO 


o 

• 




cn 


(J 

• 



Z 


X 




•. 


F 

F 










a 


CO 

*, 



LU 


CO F 




UJ 


CO l- 



0J 





*-4 

■— » 


-> 

p 









z 

3 


LU 

a 



DJ 


*— « 

LD 


r\ 


CL 


•—4 

CD 


r% 


O' 



z 

ii 


x 


L*J 







X 

CL 

o 


CL 

co 





% 

■ 


r*>. 




« 



/-> 

• 


F 




3 


ii 

II 

3 



CD 





CO 


w 

CO 



LL 


CM 





Ll 


CM 



3 

F 


G 






r% 




/-s 





•— 1 

CD 







O 


It 


U! 



a 


ll 


LU 

w 

LD 


z 



11 

CD 

— < 

3 


*— • 

LU 







ro 

F 


O 






»—» 

»— i 

3 

LL 




HH 

1—4 

3 

LL 

• 





•— * 

CM 

+ 

* 

A 

F 

3 

— » 





• • 

CO 


CO 




• 



w 

Z 

w 


• 



w 

Z 

<w* 



F 






►-* 

3 

re 

T~ 

A 

A 




g 

© 

G 


to 




X 


CO 

LL 


CO 


z 


CD 

Ll 


CO 

~7~ 


(J 



IV 

O 

1! 

w 


O 

*— • 

MD 

ID 




z 


LU 








UD 

w 

F CO 




V 

w 

F © 



3 



CM 

A 

\/ 

3 

© 

CJ 

1 

w 

W 






— 1 


f\J 




1- 





G 


z 




z 

G 

z 


e 








z 






F 










o 

LL 

O 

It 




O 

LL 

o 

II 

DC 




a 






o 

LU 

LU 




or 

(M 

-J 


-1 






Cl 


u x 


Cl 


A 


CJ 

T~ 

w 


co 



Cs 






U F 

F 




2 




-J 




z 





G 


z 





1—4 



z 


© 








*— 4 

•— i 




F 

G 


G 




►— « 


H 



Z 


M 

II 




Z LL 


O 


ii 








s 

o: 




Z 

LL 

CJ 


CJ 

- 



ll 


3 


© 

or 


Ll 




© 

Qi 



u 








© 

rv 







in 

• 






• 




CD 



• 




rv 




• 








CM 

CM 






• 

© 

• 






• 







• 








• 















• 

© 

• 






• 







• 








• 















• 

© 

CJ 

CJ 



o 

CJ 

CJ 

CJ 





uuuu 







c 

c 

c 

u 









CJ 

u 

a 

c 

c 

c 

LU 

CD 

© 

© 


cm ro in cd rv © 

CO 

© 


cm ro 

in 

CD 

IV 

© CO © 


cm ro 


in 

CD 

IV 

© 

© 

© 


OJ 

ro 

\r 

in 

•D 

IV 

CD © 

G 

CD 


IV 

Tn- 

IV 

N. 


r^ 

NNN 

© 

© © © co © © © 

©©©©©©©© 

© 

© 

CTi 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

a. 






















«-4 






CM 

CM 

CM 

CM 

CM 

CM 

0 J 

OJ 

CM 

PJ 


© 

© 

© 

© 

© 

© 

© © © © 

© 

© 

© © © 

0 CD S) 0 GD © 0 0 

© © 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 






PAGE 000? FT12 10:35 PM MON. , 3 MAP 


r 


- 177 - 


r* 

c-:> 

c 



f. 


\ 

c_ 

r* 
c •. 


CO 

2 

o 

H 

cr 


2 

i 


L«- 


\ 

S 


CO 

CL* 

p 

£ 

tr 

CL 

<r 

CL 

LU 

cn 

L'J 


CL 

CD 

ll 

UJ 

2 


P P 

s : S 
u 


U ! <D 

2 — 
* 

P *‘X> 

P w 

2 

C? LU 
CJ P 




£ 

CO 

Pi 

CE 

_J 

ZD 

CJ 

-J 

<r 

CJ 


CD 

II 

2 *: 


ro 

ro 


o 


CL 

o 

CO UJ Cm CVJ ^ 
II ZD — • UD 
Z * 

-- co »jd ud 

*w w w w 

LO!Z 

n CL a LU L'J LU 
L L' U h h h 


L»J 


CD r*. /■> /“v 

^ in co cn 
>. ^ ^ •. 
CD LD CD UD 


~ CD 
CSJ — 
% «. 
LD CD 


- /*n n ^ n co 


* * 
•JD CD 


V ~ ~ ~ ~ 


CD 


DT Cl CL CL 

3 3 3 3 


P P 

»— < *— i 

CL DL 
3 3 


LU 


p; ^ 


or or 

3 3 


P- 1 — 

* ■ * ►— < ►— < ►M 

or d; c: or 


lu 


% *, 

CD CD 


P- H 

CL CL 


* « 

CD IJD 


K P- 

c. or 


in 

o» 


oj 


ro 

ro 


Ll 

* 

<2 

* 

2 

2 


p 


> 








X 

P 


UJ 


w 





3 


p 


o 


►— < 

^-s ^ z-' v 




<E 

2 

V . 




p 

— oj ro ^ in 




CJ 

<Z 

H-t 


ci 


P 

W w W W W 



CL 


U 

II 

P 

CL 


p 

\/ S.,' s/ s/ s/ 


CJ 

*n LU 1 - — 

P 

CO 


P 

a 


X 

CL CL CL CL CL CD 


<2 

1 — Ll CL CD CL 

LU 


p 

p 

Cl 


w 

O O O O D CO 



p 2 LU Q LU 

CL 

* 

A I" 

x 

2 

2 

P 

3 3 3 3 3 CO 

LL 

X 

X •— < *-i «— • i— . 

O 

> 


P 


LU 

f 

CL 

O 

Ll 

CO 

P 

2 

UJ 

cj 

►— ( 

Ll 

Ll 

LU 


u u u u 


u u u u u u 


X 

W y-K 

o 


F 

in ~ 

u 



CL — 0J 



UJ 

O * * 

UJ 


F 

Ul CD IJ3 

2 


g: 

"T W w 

*— • 

cn 

-j 

P 

2 

cn 

3 

LU LU 

CL 

G 

CJ 

3 

_J P P 

Pi 

G 

X 

<X Or G ' 

UJ 

If 

CJ 

LJ 3 3 

& 

UJ 

• 


• 


• 


• 


• 


• * 


cj a 

c 

c 

c 

c 

c 

c 



ro\rinLDis-CDChcs^cvro^tnLDN-CQcriCD--'CjPo^rinLDN.CDCOGD-'Ajro^‘inLD>-CDcncD-'^'M 
in in Ln lo in in in ud U3 U3 ud co ud u? ud co cn r^ iv rv n. rw>- rv rv co a» od eo 03 co co cd co qj cn <ri en w, 
ojr^jcuo*j(\ir^ojcuivjcvjrvjcvjc\iCM(\j(\i(\i(\ic\JC\jrvjruruc\j£\Jcvjf\Jc\ioJrMcvjcvjcuc\joj oj oj im c\j n r\i oj 

OQQG 0 QOQQQOOQQOQOOOOOOOQSOOSQQOOOOQQOGGQOS) 


ro ^ 
cn cn 




-179- 



rvj 






w 




CD 


CL 




3 


CL 




+ 


LU 

* 

F 


t 

LU 


a 

F 

to 

z 

•"■4 



Z 

3 

M 

X 


—1 

*— • 

+ 

<X 

As 


u 

cr 

LU 

F 

CL 


C\ 

F 

0 

CL 

ll 


X 

CL 


LU 

LU 


-1 

UJ 

X 

(_> 

or 


CJ 

CJ 


3 

> 

i 



x 

DJ 

3 


cn 

3 

- - 



CD 



CD ?? 

cn 

CM 

x 

CO 

3 as 

F 

\ 

CO 

Cn 

F 

+ ro 

LU w 


A 

CTi 

V—* 

*— 


AS 

*-* 

F 


. ck 

e-3 

CD 


O 

C* 

3 CL 


w 

* 

w 


* LU 


*— t 

• 

+ 


r. -J- 

<r 

F 

or 


a 

-J C J 

LU 

_J 

<£ 

if 

LU 

LL rr 

CL 

LU 

3 

X 

Cl 

lu *::• 




rs 


CL a. 

CJ 

+ 

ro 

AS 

CJ 

1 OJ 

CJ 

A 

•— « 

U 

A, — ‘ 


in 


w 


ro w 

F 

W 


or 

F 

w or 

X 

•-1 


A\AsA\A.ASAsAS®"“» 

— ruro-tfincDcn — — 

W W W W W W WWW 

CL a CL CL £V- c_ CL or CL 

CL GD .CL' CL' Q_ CL CL CL QT 

LU UJ L ’ LU LU LU LU LU LU 


C 3 —■* 

fO — 

3 CD CD h ^ 

• W «r -1 Q r>«* 

E h H h w a: 

x « o' a uj 

3 C F 

X as CO a K a r\ 

\ ru \ CM □ CM «- 

F ^ h «. CO «. v 
* ‘X * CO II CD CO 

• W • W /-\ w w 

(M (M «— 1 

ll LU II LU LU LU 

3 F LD F ~ F F 

^ M *4 M M M 

F P' F DC CL CL OT 

« 3 O 3 LU 3 3 


J U CJ 

0 3 3 3 3 

CO x » V »III 
CO 3 * — © 3 3 3 * 

aia:u.«HixiF 



U- 




Ll 

ro 


A 

£ 


3 

w 

-J 

in 

3 

w 


CD 

O 


CD 

O 

X 

[u 


3 

X 

w 


O 

LL 

or 

u 

> 

As 

CL 

u 

CL 

As 

CD 

As 

LL 

£ 

3 

* 

LL 

£ 

CL- 

UJ 

* 

LU 

W 

\ 

* 

£ 

s 

n 

A 

CD 

w 


3 

II 

AS 

CD 

w 

Z 

o 

LO 

CD 

w 


II 

ro 


*— * 

II 

s 


3 

CO 


CJ 

o 

-J 

Ll 

w 

V 

O 

o 

CL 

Ll 

w 

g 

or 

LU 

II 

c? 

LU 

F 

-J 

LU 

CL 


-J 

LU 

or 


F 

CD 


LU 

CL 

CL 

DC 

UJ 

pr 

or 

or 

LU 

0*) 

or 

> 

> 

LU 

3 

> 

> 

UJ 

3 

A 

or 

3 


CD as — • c^J TO ’sT LH CD CO o 
CD CD CD CD CO CD CD CD CO ■CD CO 


3 LULULULJLUUJLiJUJLULULlJ 

3 or cl cl or or or or cl qt or cl 
O 3 3333333333 


CJ U U CJ CJ CJ 


O CJ O CJ CJ 


CJ CJ CJ CJ CJ CJ CJ CJ Cj CJ 


ZT Sr^^^^^T^TLnininLnLninininLnLncDCDCDCDCocDcocDCDiDrs-r^r^N-ivN-rN-iviN- 
{2KKKSKKK!2!2!2K!2EJ2!2!2!2 WMrOMror ^^^^^wMr?iro^MMnnnp5roroPi 
6)00080036505)000000000051080008006)0 00810(000000 







182 


<r 

© 


L»l 

Z 

z 


CD 

CO 

© 


r 

o 

cu 


DC 

cr 

Z 


ro 


r 

0 - 

© 

ro 


co 

tn 


© 

z 

3 

CD 

or 

CO 

z 

o 

cr 

CD 

CD 
Z 
»— • 

DC 

a 

u 

UD 


cr:* 

CD 

CD 


z x 

<Z I 
CJ '-N 
CO t— 

1—4 w ’ 

Jt 


© 
© 

UJ UJ © 

U ZD — 


LU © 

)— 


ft— 4 |r*f 


© 

CJ 

cvj 

*-• I - 

\ 

Z 

<r 


H Z 

w 

P t— 

p 

o 

© 

> a 

LU u 

A** 

f <■ 
(1 

z ~ 

co 



© 

o OJ 

© 



#<■* 

CJ 3 




CD 

in 





iv 

m 


CM 



CO 

U! 

z 


r 

© 

►“4 

Ci 

o 

CO 


a 

u 


LU 

cc 

C* 

LU 


w \ 

*— • ,*— > 

H *""• 
u w 

UJ IL 
G\ 3 
+ $? 

*— i *— * 

LL •-« 

CD ii n 

», /->, s 

II •, ^ UJ 

^-*OJ3 

w w z 

CD cn CD »- 
CM L'J UJ H 

z z z 

O M H Q 

© _! _J LJ 


CD 

CM 


CO 


© 

CD 

UJ 

cc 


UJ 

© 

CD 

a 


UJ 

CD 


tn 

<*A AA 

uj •-» ro nr 

I w w OJ 

h LO ^ lil r 

~ CD Cl P_ DC CD LL © — 

h C. O O O to 2 3 LL 

X 3 2 3 LO - I a H 

II II II II It II II II H <*s AA ,-A -v /-S -N -A A> /“> 

r\ a rv ^ a 0 •-< CJ TO nj IjO © V © CO 

— CMfonrLn©ivcQ©-‘ 


© © 


^ y cd 
© 7T Z X 


co -* 

II II 


U- 
© LU 

„ _ • DC 

© Z © H © > 

II II p II tl II 


© 
CM 

W W v» > 



LU 

LU UJ 

LU 

LU 

UJ 

LU 

LU 

LU 

© 

Ul 

LL' 

UJ 

© 

© 

UJ 

LU 

© 

© 

a: 

© 

.J _l 

3 

-1 

© 

© 

© 

_l 


3 

3 

3 

_J 

i 

i 

3 

© 

© 

QC 

© 

© © 

© 

00 

© 

© 

CO 

© 

© 

© 

© 

CO 

© 

© 

© 

© 

© 

© 

DC 

<Z 

<x © 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

<Z 

a 

cr 

r- h- ,H H H K H 

1- 

h- 

j_ h* 

H 

l- 


H 

H 

1- 

f- 

h- 



u u 

CJ 

u 





cj 

o u 

CJ 



CJ CJ u u 




















© 


CM 

ro 

nT 

tn 

© 

n- 

© 

cn 

© — 

cm ro <3* tn 

© IV CO © © -< 

cm ro nr m 

© 

(V 

co 

© 

© 


CM 

ro 

ny 

in 

© 


CD 

© 

© 

© 

© 

© 

© 

© 

© 

© 

c?» 

© 

© 

r^- n. 

0- 

fv 

iv 

IV fv JV CD © 

© 

© 

© 

CO 

co 

© 

CO 

CO 

© 

© 

© 

© 

© 

© 

C*Y 

© 

© 

© 

© 

nr 

nj 

nT 

nr 

nT 

nr 

nr 

m 

ny 

V 

nr nr 

nr 

nr 

nr nr 

nr nr nr nr nr nr 

n* 

nr 

nr 

nr 

nr 

nr 

nr 

nr 

n- 

nr 

ny 

nr 

nr 

nr 

nr 

nr 

nr 

nr 

in 

© 

© 

© 

© 

© 

© © © 

© 

© 

© © 

© 

© 

© © 

©©©©©© 

© 

© 

© 

© © © 

© 

© 

© 

© 

© 

S' 

© 

© 

© 

© 

© 

© 

© 





PAGE 0015 FTN. 10:36PM MON.. 3 MAR.. 1930 





PAGE 0016 FTN. 10:36 PM TON.. 3 1930 




PAGE 001? FTN. 10:36 PM MON., 3 MAS. , !?00 




PAGE 0018 FTN. 10:36 PM MON., 3 MAR., 198 




PAGE 0019 FTN. !0:3E PM 





0038 C 

0039 C 

0040 C... VERSflPLOT-07 PLOTTING ROUTINE 

0041 C 

0042 C 




PAGE 0003 PL 12 10:27 PM MON.. 3 MAR., 1930 


1 


- 192 - 


o 

u 

0- 

e3 


O 

U 

CL 

A 



U 

cn 

> 

<z 











> 




LU 


P 





/■A 







F 




A 

A 

u 

z 




> 










K 



LU 

=> 




F 







LU 

h- 


O 

Z 

z 

A 





u 

»- 






A 

o 


U 

NH 

o 

% 




UJ 






% 

u 


CL 

Z 

Gi 

z 

u 




A 

o 






Z 

CL 



> 

Ll. 

*-A 

LU 




a 

u 






M 




* 


r 

CL 




Z 

CL 






z 

A 



X 

tn 

> 

CO 




M 







> 

Z 


*i 

P 

*- 

a 




z 

A 






% 

© 


CD 

u 

Z 

X 

h- 




> 

Z 






>c 



t-H 

LU 

•-4 

F 

Ll 




% 

c 






F 

Li 


A 

~ A 

o 

u 

LU 

rw 



X 







u 

© 

4-S 


r~\ **■« a 

CL 

\ LU 

U 

CO 



F 

l: 


ip 




t" 

LU 

^ ro 

z 

CO ^ z 


a 


A 



u 

e 

W-4 





A 

CL 

A CO 

o 

A *“4 

u 

CD % 

CL 

Li 



LU 

LU 

rv 





A 


fL A 

CL 

•.Li Z 
^ F >• 

© 

* Z 

O 

£ 



A 

CL 


s y 




z 

CJ 

* !U 

Ll 

u 

Li — < 

LL 



* 



CL 





(J 

F co 

£ 

CL h* •. 

F cn ro 

♦— » 
F 

F 5 

r— ✓> 

CO 

IS 



Z 

*-4 

U 

u 


H* 

cn 





H 

cn -< • 

tu 

CO % 

LU 

% 



z 

CO 






z 

CO •— 

Z 

*-• a 

Ci 

~ © 

z 




X 



% 




© 

CD 



* © i 

o 

a a 


ro 



A 

Ll 


© 





»—* 

^ © 

o 

o © * 

LU 

© CD 

u 

• 



© 

LU 

© 

© 




© 

CL 

© © 

CL 

B ^ — 


© * 


CO 

a 


/-> A 

U 

© 

ro 




A 


© CM 


v , * % 

V — 

z 

LU 

<— > 


/-» o 


rj 

« 


CM 


^-4 

A 

CM , • 

• 

< H Z 


* * 

© 

z 

— # 



A 

% 

K 


\ 



Z 

F 

F 

X > O 

Ll 

> z 

1—4 

•— « 

i 


1 *-4 

Z 


•s- 


X 


© 

© 

X > 

X 

LU LU CO 

O 

F u 

A 

u 

►H 


M % 

© 

LU 

LU 


© 


© 

© 

LU LU 

tu 

* % *— « 


• cn 

o 

A 

W 


w CM 

o 

* 



r: 



CL 

A A 

^ Li a 

h 

s/ —• 

cn 

Li 

X 


* a 

CM > 

Gi 

\ • 

w 


5~ 


A 

or 

Li iL 

u. 

CO CO 

LU 

CO • 


cn 

cm 


C 

CO 

CO 


X 


F 

© 

CO CO 

o 

*— * *— 4 

CO 

— > 

u 

»-« 

+ 


+ * 


►—« 


© 

/-s 


S’ 


M 


A A LU 


A F 

<r 

A 

•—4 


CM X 

CO 

A 

A 

2 

•-4 


LU 

cn 

A A 

H 

•-••-•.© * 

A 

•— • A 

u 

►— i 

CD w 


w w 

h- 

•»< 



<W' 


% 

h* 

*-*4 

UJ 

» a © X 

CL 

* X 


% 

• to 


CO 

Z 

A 

% 

* 

I— 


X 

z 

A A 

cn 

— — LO LU 

•— « 

— LU 

h 

*-* 

«— < LU 


LU LU 


•— t 



X 


LU 

«—< 

•»< 


w w I w 

z 

w w 

LU 

w 

n Z 

•— t 

Z Z LU 

o 

W 

W 

II 

LU 

LU 

w 

o 

W W 

A 

Z 

F 


CL 


»-n •— t 

w 

*•* «■* © 

CL 



♦—* 

II 

— > 


CL 


CL 

U U (t UJ 


U LU 

O 

U 

u x 

u u z 


U 

U 



Z 

LU 


u u 

•-« 

£ & tn 5 

LU 

pc 

££ 

LU 


—■ ii 
CM 

II 

II *-4 

A u h 

© 

© 

$ 

«. ^ 

CM 

CM 

w 

P 

Z 

© 

© 

s ^ 

F 

LU LU •-» U 

F 

LU U 

LU © 

•— 

CM CM U Z 


LU 

LU 

K 

z 

u 

—4 

LU LU 


n 




ll 

O w 

w 

^ © O 



O 

X 

o 




A 

-J u z u 

A 

u u 

A 

U <*N 

A X X > U U 

A 

U 

U 

A 

LU 

CJ u 

A 

u u 

a 

uuuu 

<r 

u u 

<r 

u -« 




© 

u 

u 




u 

© 

u u 

LU 

a a cn a 

LU 

<x <r 

LU 

<r w 




LU 

a 

a 




a 

LU 

© © 

04 

o o *-* u 

Di 

u u 

Q£ 

u > 




CL 

u 

u 




u 

CL 

u u 


CM 


CM 

rvj 


uuuu 


ouuu 


uuuu 


uuuu 


uuuu 


inuDN.ODO'iCD^cvjroTfinujrs-oomcD^rMroTtinujrs-CDcno-^fvjro^JtnuN.CBcrio — c\j ro ^ in co 
oocococoa> ova cr. cr. cncr»A a cncn© ©©©©©©©©©-* — — -^ — ~~~-<*- cm cm cm cm cm cmoo 

©C3©CDG)CD0DCD©GD<S)OC3CDC3~*~ MM " a 4**-" • ^ m «•« »> « ~ ~ ~ « 

© © © © ©©©Q©CD©©©Q©©© © (2(21SEG3CD!SO(S©®(S1S©Ci®0©SOOCjCDC3 


PAGE 0004 PL 12 10:27 PM MON.,. 3 MAR., 19S0 




PAGE 0005 PL 12 10:2? PM MON.. 3 MAP. , 193 




PAGE 0006 PL 12 10:2? PM MON.. 3 MAR.. 1930 






-196- 


3. Histogram and Distribution 

This step had to be divided into three programs due to memory limita- 
tions. 

FT16 - obtains the residuals from FT12 and constructs the histogram arrays. 
It also fits the histogram with a Gaussian and determines the distribu 
tion of data points within the histogram. 

PL16 - Plots the histogram and the Gaussian. 

PL17 - Plots the distribution. 
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Quantitative flow visualization in sodium vapor seeded 
hypersonic helium a) 

R. B Miles, E. Udd l b > and M. Zimmermann 

Department of Aerospace and Mechanical Sciences. Princeton University. Princeton. Mew Jersey 0H540 
(Received 29 August 1977, accepted for publication 4 January 1978) 

Sodium vapor seeded into a hypersonic helium flow through an orifice on a model is used at a tracer to 
determine velocity profiles A narrow linewidth tunable dye laser is scanned across the Doppler shifted 
absorption line of the sodium and the resonant fluorescence is observed Curves showing the velocity- 
shifted pressure- and temperature-broadened absorption lines at several points in the flow are presented. 
Fluorescence is strong enough to use this technique for quantitative flow visualization: at each laser 
frequency, a specific velocity component is highlighted and easily observed by eye 

FACS numbers 47 40 Ki. 32.70.Jz, 51 70. + f, 06 30 Gv 


We report here the development of a new technique 
for the measurement of high -speed flow parameters in- 
cluding velocity, pressure, and temperature, and the 
quantitative visualization of the flow field. Atomic 
sodium vapor seeded into the flow has a Doppler -shifted 
temperature - and pressure -broadened absorption pro- 
file at any point within the flow. An expanded narrow - 
linewidth laser beam will cause all regions of the seeded 
flow with a specific velocity component in the illumina- 
tion direction to fluoresce brightly enough to be easily 
seen with the eye. Tuning the laser frequency selects 
different velocity profiles. 

Fluorescence for flow visualization has been found 
useful for both qualitative and quantitative studies in 
incompressible fluid flows into which fluorescing dyes 
are seeded and observed under laser illumination . 1,2 
For the gas phase, researchers have relied on electron - 
beam fluorescence, 1 electric -discharge -induced fluo- 
rescence, ‘ and index -of -refraction techniques, r ’ such 
as Schliercn, shadowgraphs, and interferometry for 
flow visualization. Of these methods, electron -beam 
fluorescence has also given quantitative density, tem- 
perature, and velocity information in low -density 
gases,' and index -of -refraction techniques have given 
density information in high -density gases. The standard 
laser Doppler velocimeter is limited in low -density 
high-speed facilities by particle lag and low signal 
intensities. 7 The technique presented in this paper 
applies to both high- and low -density flows and has vir- 
tually no particle lag problem. 

In this experiment a concentrated high -temperature 
mixture of sodium and helium is injected through a 
model into a Mach 16 helium flow and illuminated by a 
narrow -linewidth tunable dye laser. As the laser fre- 
quency is scanned through the hyperfine lines of the D 2 
transition at 589 nm, the absorption profile of the 
sodium line is determined by observing the resonant 
fluorescence. The choices of a hypersonic helium flow 
and injection through the model were made for sim- 
plicity since helium does not quench sodium and the 
model-mounted injector is easy to fabricate. The tech- 
nique is amenable to high-speed dry nitrogen and air 
flows, and full flowfield seeding appears to be feasible. 
The use of tunable lasers to observe sodium lines is 
well developed, H 9 and the theory of this technique has 
been discussed previously. 10,11 


An argon -ion -pumped Spectra Physics 580 dye laser 
beam is expanded through a spatial filter and telescope 
and passed through the wind tunnel at 40 to the flow 
axis. The laser frequency is set first using a spectro- 
meter, then a hot sodium cell, and finally an atomic 
beam. An optical spectrum analyzer monitors the fre- 
quency stability of the laser. Fluorescence is observed 
when the laser frequency is tuned to a Doppler -shifted 
absorption line of the sodium. The flow velocity com- 
ponent in the illumination direction is found by a 
straightforward comparison of the fluorescence versus 
frequency from the seeded flow with that from the 
atomic beam. 


• 

l 


— 



• ' 




FIG. 1, ComiKisitc photograph of the seeded flow. From top to 
lx>ttom the frequency' is increased, thus highlighting sequen- 
tiolly higher-velocity components, I\>int A is the injection 
orifice: i*>ints B and C are the locations of the velocity sur- 
veys 5,3 mm and 18 mm downstream of the injection point. 
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A 10° half -angle cone model is mounted in a 1,0x10° 
N/m 2 room -temperature -plenum Mach 16 helium tunnel 
for this experiment. A heated sodium oven is located 
inside the model, High-purity helium flows at 0.01 
g/sec with a back pressure of 1,3X10° N/m 2 into the 
oven through a resistively heated stainless -steel tube 
and mixes with the sodium vapor* The sodium -helium 
mixture then passes through another resistively heated 
stainless -steel capillary and is injected into the flow 
through a 0.9-mm-diam orifice on the model, The 
sodium density reaches self -trapping levels at an oven 
temperatuie of 220 °C, so normal operation is in the 
region of 190 °C. The capillary is maintained at approxi- 
mately 500 °C to avoid clogging. 

With the flow, the hot sodium -helium jet is swept 
downstream, Fluorescence can be seen all the way 
along the cone and continuing down the tunnel. As the 
laser is tuned, different velocity components of the flow 
are highlighted. Figure 1 is a composite of four pictures 
of the seeded region illuminated at different wave- 
lengths, The corresponding highlighted velocity compo- 
nents vary from low velocity in the top picture to high 
velocity in the bottom picture. 

For velocity measurements, a cylindrical lens is 
placed in the beam and the laser focused to a line over 
the axis of the cone. A detection system is focused to 
intersect the line so a point in the flow directly over 
the cone is observed. The calibration cell scan of the 
D 2 hyperfine lines shown in Fig. 2 is compared with the 
scans from points 0,25 mm, 0.75 mm, 1,25 mm, 1,75 
mm, 2, 25 mm, and 2.75 mm above the cone 5. 3 mm 
and 18 mm downstream of the injection point. In Fig. 3 
the scans 18 mm downstream are shown normalized, 
vertically offset, and scaled so that 0 MHz corresponds 
to no Doppler shift. The overall line shape is a com- 
bination of the temperature- and pressure -broadened 
profiles of the six hyperfine transitions of the D z mani- 
fold. These transitions are used because the strong line 
(3 2 S 1/2 [f r =2] to 3 2 / > 3/2 [Ft=3]) seen in Fig. 2 provides a 
convenient reference wavelength. 

The axial velocity at a point may be determined by 
measuring the frequency shift of the peaks and using 
simple geometry. Figure 4 gives the axial velocities 



FIG, 2. Atomic beam calibration cell scan of the hyperfine D 2 
lines of sodium vapor. 
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FIG. 3. Fluorescence as a function of laser frequency off reso- 
nance from points 0,25 ± 0.25 mm, 0. 75 ± 0, 25 mm, 1. 25 
± 0.25 mm, 1.76 ±0, 25 mm, 2. 25 ±0.25 mm, and 2, 75 ±0. 25 
mm from the surface 18 mm behind the sodium-helium injec- 
tion point. The top curve is closest to the surface. The error 
bar on the top point of each curve applies to all points In tile 
curve. The curves are sketched through the points for clarity, 

at various heights above the surface 5.3 mm and 18 mm 
downstream of the injection point. The error bars are 
large because the detector aperture sampled a 0.25- 
mm-diam region and the data reduction was done by 
hand. The scattered signal was bright enough to require 
attenuation before entering the detecting photomulti- 
plier, so volumes several orders of magnitude smaller 
should be observable. 

The observed Voigt linewidth is a convolution of pres- 
sure and temperature broadening of the sodium vapor 
plus additional terms due to turbulence, the linewidth 
of the illuminating laser, and the natural linewidth of the 



FIG. 4. Measured velocities above a point 5.3 mm (circles) 
and 18 mm (dots) downstream of the injection orifice. 
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FIG. 5. Normalized computer-gener- 
ated fit to the experimental points taken 
18 mm downstream of the sodium-heli- 
um injection point and 1.25 mm above 
the surface. 


sodium. The calibration cell scan (Fig. 2) shows ttiat 
the combined laser and natural full widths at half- 
maximum are not broader than 25 MHz and therefore 
can be ignored relative to pressure and temperature 
broadening. Pressure broadening of sodium by helium 
is 12 

A/,=4.2X10 10 P/P 0 (T 0 /T) l/2 , 

where P 0 and T 0 are the pressure and temperature at 
STP and A/^ is the full linewidth at half-maximum. The 
full width at half -maximum of temperature -broadened 
sodium is 

A/^i.25xio°(r/r 0 ) 1/a . 

A normalized computer -generated fit to the experi- 
mental points at 1.25 mm above the surface and 18 mm 
downstream is shown in Fig. 5. The best fit yields a 
Voigt linewidth of 1396 MHz made up of a pressure - 
broadened (Lorentzian) component of 416 MHz and a 
temperature-broadened (Gaussian) component of 1157 
MHz. The corresponding temperature and pressure are 
233 °K±8%and 9 Torr ±12%. The sensitivities are de- 
termined from the variation of pressure and tempera- 
ture among acceptable fits and are strongly dependent 
on experimental noise. The temperature drops to 190 °K 
±15% 2 mm farther into the flow; however, the pres- 
sure remains approximately the same. Fits at other 
points show the temperature decreasing farther out into 
the flow; however, due to increased noise ir* our mea- 
surements, values of temperature and pressure could 
not be determined for these points . 

The measurement of velocity close to the cone sur- 
face under these circumstances is difficult by other 
techniques. Pitot surveys tend to interfere with the flow 
and require the temperature and static pressure in 
order to deduce the velocity. Electron -beam fluores- 
cence is hampered by scattering of electrons from the 
surface. 


This technique virtually eliminates the particle lag 
problems associated with standard LDV and opens 
exciting possibilities for measurement of multiple 
velocity components, turbulence, temperature, pres- 
sure, and density, Since seeding through the model 
perturbs the flow, noninteractive flowfield seeding 
techniques are currently under development in our 
laboratory and the expansion to dry nitrogen and air 
flows seems promising, 

The authors wish to thank G.K. Bienkowski, S.M. 
Bogdonoff, F.R. Hama, S.H. Lam, G.S. Settles, and 
particularly J. A. Smith for technical and editorial sug- 
gestions. W.C. Honaker, J.C, Hoppe, andW.W, Hun- 
ter at NASA Langley have interacted closely with this 
work, giving strong encouragement and suggesting ex- 
perimental configurations. 
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Resonant Doppler velocimeter 

Richard B. Miles 

Aerospace and Mechanical Sciences Department, Princeton University, Princeton, New Jersey 08540 
(Received 15 October 1974; final manuscript received 26 February J975) 

Narrow linewidth tunable lasers augur a new kind of laser Doppler velocimetry employing resonant 
absorption and fluorescence from trace atomic species rather than scattering from particles. This 
technique may provide better turbulence and small volume information than present velocimetry. 


A narrow bandwidth tunable dye laser incident in the 
direction opposed to a flow seeded with trace amounts 
of sodium vapor will be absorbed at the Doppler-shifted 
resonant frequency of the sodium and the light will be 
rapidly reradiated in all directions. The fluorescence 
intensity may be correlated with the absorption; the 
maximum fluorescence occurring when the tunable laser 
is matched to the Doppler -corrected resonant absorp- 
tion line center of the observed region. Measurement 
of the laser frequency, then, will yield the velocity of 
the gas in the observed region. Pressure and tempera- 
ture discontinuities will also be observable as intensity 
variations, rendering this tool a convenient device for 
qualitative flow visualization as well as quantitative 
measurements. Similar applications have been sug- 
gested for the study of plasmas, 1 and experiments using 
a tunable dye laser to study hyperfine structure in a 
sodium atomic beam have demonstrated the feasibility 
of the technique. 2 ' 3 

The Doppler shift for sodium is 1.7 MHz for each 
meter/sec of flow velocity. Since the natural linewidth 
is 9.8 MHz, the limiting resolution of the resonant 
Doppler technique is several meters/sec. Pressure 
and thermal broadening will substantially decrease the 
minimum velocity sensitivity. Several experimental 
configurations with these minimum velocity resolutions 
and other pertinent information are given in Table I for 
helium and nitrogen. 

Laser stability is not expected to be a limitation. 
Bjorkholm and Ashkin 4 achieved stable operation to 
within 20 MHz for times up to 30 sec and were able to 
repeatedly scan through the sodium resonance lines. 

Wu et al. 5 were able to operate their dye laser with a 
short term jitter of 200 kHz and observe the hyperfine 
structure of iodine to one part in 10 9 . Commercially 
available lasers now produce single mode electroni- 
cally scanned beams at the sodium resonance frequen- 
cies with scan rates up to 40 msec and resolutions of 
less than 10 MHz. 

The hyperfine structure of the sodium is on the order 
of the frequency shifts to be measured and therefore 
must be taken into account. Splitting of the 3 2 P 1/2 level 
is 192 MHz and, if unresolvable, will cause a 100 m/ 
sec indeterminacy. Splitting of the 3 8 Si/ 2 level is 1.77 
GHz and should pose no problem. Either splitting may 
be used to calibrate the instrument. A sodium atomic 
beam calibration cell may be used in which P 1/2 state 
splitting is easily observable. 2 * 3 The wind tunnel flow 
will be both pressure and temperature broadened so 
p Uz state splitting will generally not be observable; 
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therefore, the most desirable line to observe is the 
3 2 S x / 2 (P-l) to 3 z Pi/z (P = 2) line since its companion 
is a factor of five lower and should not cause difficulty 
in accuracy. 

Rows 1-4 of Table I list operating conditions in typi- 
cal wind tunnels of interest. The velocity sensitivity in 
Row 5 is determined from the pressure and thermally 
broadened linewidth. Full width at half maximum val- 
ues are 

/ fcrW 2 / T W 2 

" 757 W) — . <*> 

where $, = 0.89X10 1 ° for N 2 , 6 5, = 1.27X10 l ° for He, 7 
and 5 n = 9.76X 10 6 is the sodium natural linewidth. N 0 
= 2.69X10 19 molecules/cc and N is the static number 
density of nitrogen, / 0 = 5. 0B5X 10 14 Hz. A square root 
dependence on temperature is assumed, although at 
very low temperatures this may not be accurate. The 
velocity sensitivity is rather arbitrarily chosen to be 
± 10% of the full width at half maximum of the combined 
thermal and pressure (Voigt) linewidth profile plus the 
natural linewidth with which the fluorescence is being 
compared. 

• (3) 

Pressure effects cause a shift of the resonance line 
to the red in the case of nitrogen and to the blue for he- 
lium. 7 

% s =-0.3«xl0“i(X)‘" , (4) 

5/,..0.073><10» ] |(^)''\ (5) 

These shifts produce a velocity correction factor shown 
in Table I. 

Maximum collected signal intensities will depend on 
the incident laser power density, the sodium atom den- 
sity, static temperature and pressure, quenching, and 
on the collection optics., The incident laser power den- 
sity should be held below the saturation level of the 
sodium, particularly for pressure dominated broaden- 
ing since saturation further broadens the resonance. 
Assuming the 3 2 S 1/2 (P = 1) to 3 2 P 1/2 (P=2) line con- 
tributes 5/48 to the total 3S-3P oscillator strength / 

= 0.982,® the line center cross section for this naturally 
broadened line is cr= 1.77X 10' 10 cm 2 . The unquenched 
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TABLE I. Velocity sensitivity and signal level for various wind 
tunnel configurations. 


Helium 

j Nitrogen 

Much 

mmibar 

8 

10 

10 

20,7 

2,9 

to 

20 

Stagnation 
pressure (Atm) 

20 

108 

20 

102 

0,8 

50 

100 

Stagnation 

temperature 

(*K> 

200 

230 

300 

297 

291 

1000 

2000 

Free stream 

velocity 

(m/sec) 

1723 

1000 

1753 

1748 

010 

1407 

2027 

Velocity 

sensitivity 

±3% 

*7% 

±0.6% 

±0,0% 

±49% 

i 2,3% 

±1,1% 

Velocity 

correction 

(m/sec) 

+ 27 

+ 07 

+ 1.8 

+ 3,3 

-1128 

-9,3 

-0,2 

Fluorescence 
per atom 
(photons/soc^ 

0.33 

xiO* 

2.58 

X10 4 

3,08 
xi O'- 

3.30 

xl0 f ' 

2,45 

XIO 2 

1.08 

xlO 5 

1.99 
xlO ‘ 

Maximum 
sodium density 
(atoms/cc) 

4.09 

XIO 10 

1.15 

X10 n 

8,00 

X10 9 

0.01 

XIO 9 

2,91 

XIO 11 

2,05 
X 10 10 

1,47 

xlO 10 

Maximum 
detected signnl 
from 1 cc (W) 

2.5 

xl0* fl 

2,5 

X10' fi 

2.5 

Xl0‘ fi 

2.5 
x 10' e 

0,0 

,0’ 8 

1.9 

X10“ C 

2.5 

xlO" 11 

Mass flow 
ratio 

1.7 

X10 7 

2,3 

XlO 7 

1.3 

X10 T 

2,8 

XIO 7 

G, 1 
X1Q 7 

1, 1 
XlO 7 

5.1 

xio 5 


lifetime of sodium is 16.3 nsec so the saturation in- 
tensity is approximately 0,06 W/cm 2 . Quenching and 
broadening will increase the saturation level? however, 
for maximum signal calculations, an intensity of 0.01 
W/cm 2 was assumed for convenience. 


The number of photons scattered per atom, n , is 
determined by the cross section at the maximum of the 
Voigt broadened line, the incident flux density, and the 
quenching factor Q: 


(0.01 W/cm 2 ) 

h/ Q 


o(/=/o)<?. 


( 6 ) 


Q = [l + 72.1(7’/lOO) 1/2 7V/AT 0 ]" 1 for nitrogen 9 and <?= 1 for 
helium. The Voigt broadened maximum may be de- 
ter mined from the thermal linewidth A/* and the pres- 
sure broadened linewidth A f p by 


<r(/=/ 0 ) = 


5 fe\ (ln2) 1/2 
48 2m e 



4 —“SH 

x^-. r ,245*i&)<cm*> 


where /is the oscillator strength, e is the electron 
charge, 7 ] is the free space impedance =377, and m % is 
the electron mass. 

From the number of photons scattered/atom, the 
maximum signal may be determined for any experi- 
mental arrangement. Collection optics, of course, only 
subtend a small portion of the scattering angle. The 


// 5 optics assumed here reduce the number of photons 
seen by 400; however, the signal varies as the inverse 
square of the collect /number. The density is limited 
by the condition that the medium be optically thin. This 
is specified by limiting the absorption length to longer 
than the wind tunnel distances: 1/N v jj^L I has been 
chosen as 10 cm for sodium density and mass flow ratio 
calculations presented. The scattering volume ob- 
served is assumed to be 1 cm 3 , but the detected inten- 
sity simply scales with the volume. The large calcu- 
lated intensities suggest that very small volumes or low 
atom concentrations are observable. This is consistent 
with observations by Fairbank et al. f 10 who have mea- 
sured sodium densities as low as 10 5 atoms/cc with a 
signal -to -noise ratio of better than 1000. 

It is apparent from Table I that substantial inaccura- 
cies occur for low Mach number of high pressure con- 
ditions. This limitation is due to thermal and pressure 
broadened linewidths causing uncertainty in the mea- 
surement of the velocity associated frequency shift. 
Heavier vapors or gases may be used to overcome the 
thermal broadening; for example, iodine vapor will re- 
duce it by 70%. The lower cross sections of these : 
molecular gases may be compensated for by higher in- 
cident laser power densities. The major loss in per- 
formance will be the higher tracer gas densities re- 
quired to offset the thermally distributed ground state 
populations. To achieve similar signal levels in iodine 
vapor, the incident laser power density would have to be 
increased about 3000 times, and the iodine vapor den- 
sity would have to be 1000 to 4000 times greater than 
the sodium atom density. 11 Since iodine and chlorine 
do not react with air, either may be a reasonable can- 
didate for supersonic low pressure air wind tunnels. 

I would like to thank F. Hama, J. Smith, and I. Vas 
of the Gas Etynamics Laboratory, Princeton University 
for many helpful discussions. 
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SUMMARY 


RESONANT DOPPLER VELOCZMETER 


Richard B. Miles 
Vrincaton University 
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Resonance fluorescence is a potentially useful tool for studying hypersonic flow- 
ing gases* This paper discusses the measurement of velocities and demonstrates visual- 
isation techniques using sodium seeded into helium flows. A narrow linewidth dye laser 
is tuned onto the sodium resonance line* The doppler shift of the absorption frequency • 
yields the velocity profile, and the fluorscence provides direct flow visualization. 
Other flow properties such as the turbulence frequency, temperature, etc* may also be 
Measured. 

1* Introduction * 


Tracer particles in a flowing gas have extensively been used for laser doppler 
velocimetry in subsonic and transonic flow. In the hypersonic regime, however, these 
particles tend to lose their ability to follow the flow* This "particle lag" becomes 
acre serious as the flow velocities increase and makes standard laser doppler velocimetry 
a rather inaccurate tool for monitoring low density high velocity flows and high tur- 
bulence frequencies. To overcome this problem, smaller particles must be used. As the 
size of the particle decreases, however, its effectiveness as a visible light scatterer 
also decreases, thereby lowering the signal levels until finally with particles less 
than .1 micron the signal levels virtually disappear. Thus, the ideal tracer particle 
la useless as a scatterer. 

Work currently underway at Princeton University's Gas Dynamics Laboratory is in- 
vestigating the possibility of seeding hypersonic gas flows with atomic sodium and using 
resonance fluorescence rather than Mie scattering to generate high signal intensities. 

The strongly resonant sodium D lines have been studied since before the turn of the 
century. These lines are so strong that number densities as low as several thousand 
sodium atoms per cubic centimeter have been detected in the laboratory. 1 Sodium atoms 
introduced into foreign gases such as helium and nitrogen quickly equilibrate thermally 
with the foreign gas and tend to follow the flow closely. Currently available narrow 
linewidth dye lasers can be used to probe the sodium tracer atoms. Since the absorption 
frequency is doppler shifted in the direction of the incident light, monitoring the 
fluorescence intensity as a function of the incident laser wavelength gives the velocity 
component in the incident light direction. The wavelength of the fluorescence may also 
be monitored to give the component of the velocity in the observation direction. Other 
Characteristics of the sodium absorption and fluorescence profile give additional in- 
formation about the gas mixture* for example the fluorescence intensity may be related 
to the number density, the temperature and the pressure. Illumination of the entire flow 
timlA will produce fluorescence which may be useful for flow visualization. 

Zn the laboratory, (Figure 1) a beam from a narrow linewidth tunable dye laser 
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is split into two components . The first is directed axially down the wind tunnel and 
fluorescence from that beam is monitored in a transverse direction. The second beam is 
directed into a sodium atomic, beam calibration cell. The fluorescence from the sodium 
atoms mixed with the helium or nitrogen gas flow within the wind tunnel is compared to 
the fluorescence of the pure sodium atomic beam in the calibration cell. By sweeping 
the frequency of the dye laser, the frequency offset of the fluorescence maximum from 
the wind tunnel will give the velocity of the gis mixture within the tunnel. The ob- 
served linewidth difference will give information about the pressure and the temperature 
Of the mixture within the tunnel. If the dye laser is tuned onto the maximum of the 
fluorescence profile from the tunnel, fluorescence intensity fluctuations can be corre- 
lated with turbulence. Various single and two photon processes may be conceived to 
measure other parameters of the flow. 


XI. 


Discussion 


Flows seeded with particles depend, of course, on the Mie scattering cross section. 
This cross section is strongly dependent on the size of the particle, and for particles 
less than .5 \i in diameter, falls off at approximately the fourth power of the diameter if 
ths flow is illuminated with visible light. 2 on the other hand, the ability of a seed 
particle to follow the flow is strongly enhanced as the size is decreased. For example, 
for water droplets in air, the relative mean square error at 10 KHz drops from 2.7% for 
a l.Op diameter droplet to .02% for a .3p diameter droplet. 3 In hypersonic flows, tur- 
bulence frequencies up to a megahertz or more must be measured. Even at 100 KHz, the 
• 3 v particle cannot follow the flow, 4 and at 1 MHz the true frequency spectra cannot be 
measured. To overcome particle lag and low scattering cross sections, very high densi- 
ties of < .3y particles must be used. Such high densities are prone to particle agglom- 
eration and may also affect the character of the flow. 

An atomic species seeded into the flow will equilibrate in several collisions, 
therefore there is virtually no particle lag. Although the resonant scattering cross 
section is low (a » 10~ 12 cm 2 ) , the mass of the atoms is small so large numbers (>10 10 /cc) 
of tracer particles may be used without affecting the flow characteristics. 

Due to the nature of atomic fluorescence, the gas mixture, and the illuminating 
laser, various broadening processes are observed which limit the accuracy of the measure- 
Mint. lii - ut; LV-l;U il* OF m 

A. Lifetime broadening. ORIGlr.A — v ,yj ^ ^ 00R 

The larger the absorption cross section, the shorter the lifetime in the 
excited state. Thus an atom will re-emit the radiation within several tens of nano- 
seconds. This means that the "observation time" cannot be longer than several tens of 
nanoseconds . 

B- Collision broadening. 

In a gas mixture, a collision with another particle dephases the seed atom and 
causes the "observation time" to be reduced to the time between collisions. Even if the 
atom remains in the excited state, these collisions cause the energy level to be more 
uncertain and therefore broader. 

C. Thermal broadening. 

At any instant of time, each atom has some instantaneous velocity. The dis- 
tribution of instantaneous velocities is centered around the average velocity; the width 
Of the distribution is determined by the temperature and mass of the atoms. Since the 
absorption line center for each atom is doppler shifted by the instantaneous velocity, 
ths absorption profile for the mixture is broadened. 

D. Laser broadening. 

Our ability to measure the absorption profile is limited by the bandwidth and 
stability of the laser. Current tunable dye laser systems have output linewidths in the 
10-20 MHz range. Some dye lasers have been stabilized to 250 KHz 
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E. Velocity fluctuation broadening. 

Fluctuations in the average velocity will be caused by turbulent eddys. 

These fluctuations look similar to the thermal broadening and will cause the absorption 
profile to broaden yet further. T?v'sy differ, however, in that the thermal broadening 
occurs in time comparable with the atom-atom collision frequency, and velocity broadening 
occurs in time comparable with the turbulence frequency. 

Zn general all these ambiguities are mixed together to produce some measured ab- 
sorption linewidth with its peak shifted by the average velocity. Although they may make 
the accurate determination of the velocity difficult, thermal, collision, and velocity 
broadening yield additional information about the flow. 

XXI. Example - Sodium in Helium 

The strong sodium D lines are split into ten hyper fine lines as shown in figure 2.* 



Figure 2 Sodium Resonant Line Strengths 
and Energy Levels 


Theoretical line strengths and spacings are shown in the figure. With tunable dye lasers 
and atomic beam cells, these lines have been measured down to 10% of their natural width 
of 9.6 MHz. ® In the flow, these lines will be substantially broadened by collisions, 
thermal broadening, and velocity fluctuations. A calibration cell is necessary to pro- 
vide a comparison between the unperturbed line and the wind tunnel. 

Although any of the lines will suffice, two are of particular interest. The 

2 2 

S 1/2^ Fss ^ to P 1/2^ F=2 ^ line is paired with a line with an intensity lower by a factor 
Of five and is removed by at least 1.77 GHz from other strong lines. Thus this line pro- 
vides an easy unambiguous reference. A problem occurs, however, if the atom is excited 
numerous times within the observation region. Since the selection rules permit AF=0, ± 1, 
after several excitations, the probability is high that the atom will relax to the 
S 1/2^ P ~ 2 ^ state and effectively be lost. This phenomena is called resonant pumping, and 
to avoid it the 2s x/2^ F=2 ^ to 2p 3/2^ F ~ 3 * line may be chosen. Any atom reaching the F * 3 
State can only relax back to the F = 2 state from which it came. 

Assuming the 2 Sjy 2 (F=l) to 2 Pjy 2 (F=2) line is observed, the line center value of 
the absorption cross section is* 


o(f.f o > 


2.55 x 10 
Af, 


-3 


[exp (0.833 jjS-r) 


833Af 

c 


x (1-erf 


where o is in cm 2 , Af c and are the collision and thermal linewidths in Hertz. The 
flow is assumed to be non turbulent. The error function arises from the combination of 
the thermal and collision broadening into a Voigt profile. Using the S i ^ 2 (F*2) to 
**3/2 line wil1 increase this cross section by a factor of 2.8. 

Absorption is detected by observing the fluorescence from the atoms relaxing back 
to the ground level. In the absence of quenching all the photons absorbed are reradiated 
into 4 it steradians. It has been shown that helium does not quench sodium, so the de- 
tected signal can be straightforwardly calculated from the incident photon flux, the 
cross section, and the detection aperature. 

It is apparent that the scattering intensity increases as the number of sodium 
atoms increases and as the incident photon flux increases. As the density of atoms be- 
comes large, the fluoresced light is itself reabsorbed and subsequently reradiated. This 
is the self trapping or "optically thick" regime. Since the region to be observed must 
fluoresce and these photons must be able to reach the detector, the sodium number density 
must be maintained low enough to permit a photon to pass from the observation region to 
the window without being absorbed. In general, sodium densities in excess of 10 10 /cc or 
so cause self trapping over distances of 10 cm. 

Table 1 shows calculations for several helium tunnels. 6 The velocity correction 
is due to a slight shift in the absorption peak caused by the helium ambient pressure. 

The velocity sensitivity is chosen to be one tenth the sum of the full width half maxi- 
mum linewidths of the calibration cell fluorescence and fluorescence from the tunnel. 




TABLE I 

Velocity Sensitivity and Sodium Density 
•for Various Wind Tunnel Configurations 




Mach 8 


in Helium 
Mach 10 


Mach 16 


Mach 20.7 


Stagnation 
pressure (Atm) 

20 

108 

20 

102 

Stagnation 
temperature ( K) 

300 

258 

300 

297 

Free Stream 
velocity (m/sec) 

1725 

1606 

1753 

1748 

Velocity 

sensitivity 

± 3% 

+7% 

: 

+0.6* 

+0.6* 

Velocity 

correction (m/ sec) 

+27 

+67 

+1.8 

+3.3 

Maximum sodium 
density (atoms/cc) 

4.69 

xlO 10 

MS 

xlO 11 

8.06 

xlO 9 

9.01 
xlO? . 

Mass flow 
ratio 

1.7 

xlO? 

2.3, 

xlO 7 

1.3 

xlO 7 

2.8 

xlO 7 


As the laser intensity is increased into the saturation region, the collision 
broadened linewidth begins to increase as/1 + I/I gat and the peak absorption decreases 
as 1/(1 + 1/1 Thus as the intensity is pushed farther into the saturation regime, 

the sensitivity drops due to the linewidth increase. 

A saturation intensity of .06 watts/cm is calculated from the naturally broadened 
sodium atom cross section. The temperature and pressure broadening substantially de- 
crease the on-line center cross section and therefore increase the saturation intensity. 
For each configuration a saturation intensity must be calculated, and these are given 
for several helium wind tunnels in row 1 of Table II. If the laser intensity is in- 
creased to the saturated intensity, the absorption cross section drops to 1/2 its 
unsaturated value and the linewidth increases by /?. This may be considered the maximum 
signal level per unit volume before sensitivity is lost. Row 2 of Table II gives the 
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. TABLE II 

Saturation Intensities and Maximum Signal 
From the Minimum Volume 



Mach 8 

Mach 10 

Mach 16 

Mach 20.7 


5.1 

1. 23x10* 

8.5X10” 1 

9.6xl0 _1 

• (CTO 2 ) 

-12 

2.1x10 

8.6xl0“ 13 

1.24x10” 

11 l.lxlO -11 

Max Signal 
from 1 cc 
•t I-I, at (w> 

«.3xl0~ 4 

1.6xl0 -3 

l.lxlO” 4 

1.2xl0” 4 

Beam Area 
from 10 mvr 
laser at . 

1 “ x aat <«■ > 

2.0xl0~ 3 

8.1xl0” 4 

1.2xi0” 2 

l.OxlO" 2 

KKPROi I'; 
ORIGINAL i 

Smallest 
Sample 
.Volume (cc) 

l.OxlO” 4 

2.6xl0 -5 

1.4xl0“ 3 

1.2xl0” 3 

Maximum 

6.3x10”® 

4.0xl0” 8 

1.5xl0” 7 

1.4xl0” 7 


Signal (w) 


OF THIS 

P V 


on line center value of saturated cross sections for each configuration. If the intensity 
is chosen at this level, the maximum detected signal from one cubic centimeter is shown 
in row 3, assuming the number density of sodium atoms is still confined to the optically 
thin region. 

The output from a commercial tunable dye laser is approximately 30 mW. If 10 mW 
passes through the sample volume, the beam area needed to reach saturation is shown in 
row 4. Assuming the detection optics is arranged to look at a cylindrical volume with 
a length equal to the diameter, the detected signal is given in row 6, calculated from 
the sample volume given in row 5, and assuming f/5 collection optics. Thus the volumes 
shown in row 5 represent minimum observation volumes before saturation begins to dominate 
with a 10 mW laser; of course smaller volumes may be observed by simply throwing away 
some of the scattered signal. It is important to note that the detected light corres- 
ponds to a rather large photon flux and can be enhanced by a factor of 25 by going from 
f/5 to f/1 optics. 


IV. Turbulence Measurements 

From Table II it is apparent that rather small volumes may be observed without 
substantially reducing the collected signal. Sample regions are on the order of a milli- 
meter or less. 

If the incident laser frequency and amplitude are held constant, fluctuations in 
the output fluorescence intensity will be caused by temperature, pressure, density, and 
vslocity fluctuations within the flow, assuming the sodium is well mixed. The frequency 
spectrum may then be read directly from an electronic spectrum analyzer connected to the 
detector. 

Since the detected signals of 10~ 7 and 10“ 8 watts correspond to between 3 x 10*° 
and 3 x 10** photons/sec, even fluctuations in the megahertz region have sufficient photon 
flux to eliminate statistical noise. All the minimum sample volumes shown have more than 
10® sodium particles, so again scattering is not masked by low particle density statistics. 

V. Seeding 

The choice of which atomic or molecular vapor to use as a seed is presently limited 
by the availability of tunable narrow linewidth lasers. By far the most successful and 
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Figure 4 Sodium- He Hum Jet from the Figure 5 Sodium-Helium Jet from the 

Injector Oven into Atmospheric Injector Oven into 4 Torr 

Helium Illuminated by Sodium Helium Illuminated by Laser 

Vapor Lamp. on Resonance. 


of helium and sodium emerging from the oven into one atmosphere of helium. This plume 
t iJ transparent under standard illumination, but appears opaque when illuminated with a 

sodium vapor lamp. The density of sodium atoms is high enough to cause the plume to 
become optically thick and appear opaque. 

In Figure 5 the oven is placed in a low pressure (4 torr) chamber and observed 
through windows. The dye laser is tuned onto resonance and strong fluorescence from the 
I free jet is visible. Again the jet is optically thick so the laser beam cannot pass 

through. The multiple beams come from reflections in the dye laser and beam steering 
optics. 

VIII. Conclusion 

The possibility of using an atomic gas as a seed for doppler velocimetry seems 
particularly attractive now that narrow linewidth tunable dye lasers are commercially 
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dependable laser system is the argon lassr driven cw dye laser. With the several dyes 
available, the operable range is from 5400 A to 6400 A. Within this range lie the 
resonant lines of the iodine molecule extending over a broad range around 5535 A and 
the sodium atom at 5890 A. Most other resonances lie farther toward the ultraviolet and 
cannot be reached by this laser system. Tunable lasers in the ultraviolet are pulsed and 
therefore have inherent linewidth and jitter problems. 

The iodine molecular spectrum contains numerous rotational vibrational lines, many 
Of which overlap when broadened in a mixture. Due to the nature of molecular resonances, 
the absorption cross section for each line is much lower than the atomic cross section, 8 
and the absorption is further reduced by the large number of ground states over which the 
Molecule is distributed. Even with these problems,* iodine may prove desirable if sodium 
seeding is too complicated or if the tunnel uses air* 

The atomic resonance of sodium yields strong absorption and fluorescence and un- 
ambiguous lines. The problem with sodium, of course, is its reactive nature and the high 
temperatures necessary to achieve high densities. A simple oven in the stagnation cham- 
ber does not work because the sodium vapor condenses or is lost in other ways before it 
retches the throat. A more sophisticated injection system is required which either 
introduces the sodium downstream of the throat or injects it into the flow just before 
the throat. The resident time of the sodium within the tunnel must be short compared 
with the condensation time or any reaction rate with impurities. At velocities in the 
thousand meter per second range, the sodium must remain in atomic form for the order of 
a millisecond. 

X£ the gas flow is pure helium, then the loss mechanism is dominated by the forma- 
tion of dimers assuming the sodium is fully vaporized in the injection oven. Dimer forma- 
tion# however, requires a three body collision, two sodium atoms and one helium or sodium 
atom to take away the excess energy. At the sodium densities of interest, this three 
body reaction occurs slowly. 

Sodium particles of even submicron size may not be seen, but do provide agglomera- 
tion and condensation centers in the flow. Particle concentrations depend on the oven 
temperature and ambient gas pressure, but it seems certain that particles are almost al- 
ways present. Duthler et al. conjecture that particles in the range of .005p exist in 
the hot mixture and grow depending on the cooling rate of the mixture. 9 These particles 
tend to keep the vapor pressure of the sodium in equilibrium with the ambient temperature. 
At the low static temperatures in the observation region (generally below 10°K) , equili- 
brium vapor pressures of sodium are negligible. Thus the equilibrium rate must be slow 
compared to 1 millisecond for the seeding to work. 

Sodium also will react with oxygen or water impurities in the flow. Generally 
impurity levels must be kept lower than 50 ppm to avoid losing the sodium between the 
injection point and the observation region. 

VI* Experimental Configuration 

Figure 1 shows a sketch of the experimental configuration. The laser is mounted 
on a vibration isolation table to facilitate frequency stability in the wind tunnel 
laboratory. The laser fequency is swept through the sodium resonance, and the signal 
from the calibration cell is compared with that from the wind tunnel. 

. t Currently work is concentrating on the sodium seeder oven. The newest configura- 
tion is shown in figure 3. Pure helium enters at the bottom and passes through a heating 
core. The crucible is heated by the helium, and the hot sodium-helium mixture passes 
through a double walled insulated injector needle. This oven is designed to be placed 
in the stagnation chamber with the needle reaching just below the throat. 

VXX. Flow Visualization 

The usefulness of the sodium vapor as a flow visualization technique is demon- 
strated in figures 4 and 5. The fluorescent plume apparent in figure 4 shows a mixture 
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available. If seeding difficulties can be overcome, this device may be used to monitor 
not only the flow velocity, but turbulence frequencies as well. 

With* sodium vapor in the tunnel, numerous spectroscopic techniques become avail- 
able to further characterize the flow. Two photon absorption, for example, may be used 
to eliminate thermal and velocity fluctuation broadening. If the incident power density 
is increased far above saturation, the fluorescence intensity depends only on the number 
density of sodium atoms and not on the pressure, temperature, or velocity. Specific 
velocity components may be selected by two step absorption using two lasers. Other non- 
linear spectroscopic tricks may be used to look oniy at the pressure, temperature, and 
velocity. 
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